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I ' Abstract 
Oh! 

• A series of informal seminars at graduate-student level on the subject of coupling 

dependence in quantum field theory, with an elementary introduction to the notion 
^ , of resurgent function that forms the appropriate framework for the coupling depen- 

^ I dence of strictly renormalizable theories. While most of the discussion is pedagogical, 

there are also a few things for the expert: we demonstrate, by studying a model spec- 
tral integral, that an amplitude may possess both the t'Hooft set of singularities in the 
coupling-constant plane, and Borel-plane singularities of the infrared-renormalon type in 
its perturbative part, and yet be uniquely reconstructible from its resurgent symbol, the 
appropriate generalization of the semiconvergent perturbation expansion. In the same 
model we demonstrate the virtues of a quasi-perturbative expansion, obtained by resum- 
mation of the resurgent symbol in its nonperturbative direction, and which in contrast 
to the perturbative one is Borel summable. On the basis of this expansion, we discuss 
a systematic approximation method for the reconstruction of correlation functions with 
the resurgent coupling dependence typical of strictly renormalizable and asymptotically 
free theories. 
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Generalities on Coupling Dependence 



We know everything that is knowable about a quantum field theory (QFT) once we 
know the full set of its Euclidean correlation functions, 

Gn{xi...xn) = {0\ip{xi)ip{x2) ..■(p{xn)\0) , (0.1) 

for all configurations of elementary fields, or equivalently their Fourier transforms 
defined by 

" xi...d^ XN e' • ^1 + - + ''^ ■ Gn (xi ...xn) 

(0.2) 

= {27r)^5^{h + k2 + ... + kN)GNiki...kN) 

at Euclidean momenta ki. Here we denote by D the Euclidean spacetime dimension, 
which we will allow to deviate from the realistic value of D = 4, both for the purpose 
of looking at lower-dimensional model theories {D = 3, 2, 1) and for dealing with ul- 
traviolet divergence in D = 4 through dimensional regularization {D = 4 — 2e). For 
simplicity we mostly consider, instead of ( p.2| ), the momentum-space, Euclidean, proper 
vertex functions F^v, that is, the connected, amputated, and one-particle irreducible 
functions, since these are the simplest building blocks from which the full Gj\f can be 
pieced together purely algebraically. Their full set of functional dependences is 



N 



Tn (ki, k2, ... kN; / in), m (//), ^ (//) ... ; , (0.3) 



where the momenta obey the constraint T,ki = of total-momentum conservation, and 
n denotes the arbitrary renormalization mass scale introduced by renormalization in 
a scheme R not usually indicated explicitly - here, a dimensional-regularization-plus- 
minimal-subtraction scheme. Moreover, g'^ (p), m (p), ^ (p) ... are the renormalized 
coupling constant, renormalized elementary masses, renormalized gauge-fixing parame- 



ter ... of the theory, defined at the "sliding scale" fi |WEI96| and therefore "running" 
with that scale, since upon changing the arbitrary fi they all need to be readjusted to 
keep the observable quantities constant. 

In these seminars our focus will not, as usual, be on the momentum dependence 
of the functions (p. 3D, but rather on their dependence on the coupling constant g^ (n). 
(Renormalization-Group (RG) analysis shows that these two are in fact related, and we 
are going to use this relation in places, but the link is really tight and restrictive only 
for two-point functions depending on a single squared Euclidean momentum k"^, since 
for higher correlation functions the RG governs only the behavior under uniform scaling 
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of all momentum arguments but says little about individual momentum dependences.) 
We are less interested here in the dependence on other parameters such as m (n), ^ {/j,), 
and therefore will mostly get rid of those by considering massless theories (m = 0), 
adopting special choices of gauge fixing (such as ^ = 0), etc. This, of course, is not to 
say that the other dependences may not be of paramount interest in dealing with specific 
problems, but the interaction strength (fi), which makes for a nontrivial theory in the 
first place, does play a distinguished role, of which these talks are intended to provide a 
glimpse. 

To anyone with a basic knowledge of QFT, the full (/x) dependence of amplitudes 
generated by a realistic, interacting, renormalized QFT would at first sight appear to be 
of unfathomable complexity, nothing short of a horrible mess. It is therefore surprising 
that for a class of interacting theories with mild ultraviolet problems, the so-called 
superrenormalizable theories (essentially, theories in which the mass dimension [g'^] of 
g'^ (fi) is positive, and which therefore mostly live in spacetime dimensions < 3) 
this full coupling dependence can in principle be recovered from (though it is not given 
explicitly by) the most straightforward and unsophisticated ansatz conceivable, that of 
a "perturbative" power-series expansion about the free-field situation g'^ = 0: 

oo 2 / "v 

rr {{h}; a (/.); /.) = ^ a^rS^^P^^* {{h}; /z) , « = ^ , (0.4) 

p = 

where the zeroth-order terms, the "b are vertices" T^"^^ , are just the ordinary Feynman 
rules to be read off from the action of the theory. This situation will be briefly sketched 
and commented upon in chapter 1 below. We recall that an expansion of type (|0.4| ) is 
straightforward to generate formally, either by expanding the path-integral representa- 
tion of correlation functions (see ( 1.25| ) below) in powers of the interaction term of the 



Lagrangian ("integral" approach), or alternatively by iterating the coupled equations of 
motion of the vertex functions around their zeroth-order terms in a, which again are the 
-p(o)pert ^£ ( ^Q ^ ("differential" approach). Either route leads to the familiar representa- 
tion of the p-th order term '^'^^ as a sum of Feynman diagrams with the external-legs 
configuration N and p successive momentum integrations ("loops"). The main practical 
obstacle at present to actually evaluating terms of higher perturbative order p{p > 3, 
say) is the explosive growth with p of the number of diagrams contributing that will 
soon come to haunt us in the form of eq. ( |0.5| ) below, but for the superrenormalizable 
theories, which enjoy the property that overall convergence of loop diagrams in the ul- 
traviolet improves with increasing p, the time is probably not far away when computer 
programs, following one of the above two routes, will furnish us with r^^''°'^*'s for all 
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configurations of practical interest at the orders p (in the range of 10^ to 10^) that 
are nowadays routinely computed in series expansions for statistical-mechanics models. 



That expansion ( p.^ ) is at best a first step in analyzing coupling dependence has 
to do with the fact that, in the superrenormalizable and in more complex theories alike, 



the series turns out to be badly divergent. Dyson | DYS52 | as early as 1952 was the 
first to show, by an ingenious physical argument, that in the case of quantum electro- 
dynamics (QED), the correlation functions must possess some kind of non-analyticity 
at real negative a, and that consequently the radius of convergence in the complex a 
plane for expansion (|0.4| ) is zero. Later work by many authors, work which received a 



particular boost from the instanton results of Lipatov [ LIP72 |, established techniques 
for calculating the leading large-p behavior of the coefficients rj^^^*^"^*, with the results 
always fitting into the generic form 



l + O [- 



(p ^ oo) , (0.5) 



where c carries the momentum dependence, and A, b are real constants. (The large-p 
calculations form an interesting subject in themselves, into which these seminars cannot 
really enter, so I must refer you to the editor 's introduction and the reprinted primary 



sources in the edited volume ||LGZ9C|] if you wish to pursue this special topic.) In par- 
ticular, the instanton contribution to large-p behavior, which is calculated along the 
"integral" route mentioned above by expanding the path-integral representation semi- 
classically around a finite-action solution to the Euclidean classical equations of motion 
("instanton"), has the constant A equal to that finite action. That this instanton satu- 
ration does catch the full large-p behavior is, as far as I know, currently a hope. Indeed 
in the strictly renormalizable theories additional contributions have been identified in 
perturbation theory, which we shall have occasion to mention in sect. 2.4; they also fol- 
low the pattern ( p.5| ) but with different constants A, b, c. In the absence of a complete 
classification of possible contributions to large-p behavior, the mathematically minded 
person would not regard the divergence of the a expansion as proven, but since the 
pattern of eq. ( p. 51) has resurfaced in all large-p studies performed so far, we shall as- 
sume here pragmatically that up to linear combination with possibly different sets of 
constants, but with the p! factor common to all terms, this relation does give the large-p 
behavior. The factorial growth, which in the superrenormalizable theories simply re- 
flects the fast increase in the number of combinatorial possibilities of drawing a p-loop 
Feynman diagram, precludes convergence of the series (p.5D at any a. 



It is then obvious that for coupling dependence, expansion ( p.5| ) is not simply "the" 
answer. To make some use of the information contained in it, one essentially has three 
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options. If a regime exists where the running couphng a (/Li) is small, one may exploit 



property of semiconvergence (section 1.1) to halt the 6xpctnsion set a finite p — Pmax 

and 

use the resulting polynomial approximant directly as an estimate for the function r^r - 
an estimate whose accuracy, while limited in principle, may in practice be quite satisfac- 
tory. This is what Poincare, in chapter VIII of his celebrated Nouvelles Methodes de la 



series, referred to as la regie des astronomes - and "les astronomes" were, of course, 
the perturbative physicists of his day, the practitioners of celestial mechanics calculating 
planetary-orbit perturbations. This is also what is routinely done in the electroweak 
theory, where the gauge couplings are indeed small in the currently accessible range of 
scales, /u ~ f 250 GeV (the electroweak scale). 

In regimes where a (fi) is not small enough to provide good semiconvergence, one 
may try to adduce physical arguments to justify, for particular kinematical situations, 
an estimation of F^v by "partial resummation" , i.e., the exact summation to all orders 
of certain subseries of (|0.4D whose terms one has reason to believe will dominate the 
particular phenomenon. In most cases the series summed are subseries of nonzero con- 
vergence radius - geometric, exponential, or Bessel-function series - which stay clear of 
the p! growth problem. Examples are the "leading-logarithms" geometric resummations 
one performs (as discussed in sect. 2.1 below) in order to trade the coupling a (^) for a 
coupling a (Q^) running with some momentum variable (one of the kf or ki ■ kj or a 
combination of these) of r^r, or calculations of bound states like positronium by solving 
a Bethe-Salpeter equation with low-order irreducible kernel, which amounts to summing 
an integral-kernel-valued geometric series in a. Such summations are useful and (as in 
the bound-state problem) sometimes necessary to extend the range of applicability of 
low-order perturbative information, and to the extent that the resummed function often 
has properties qualitatively different from those of the finite polynomial approximants, 
they may be called " nonperturbative" . One should however keep in mind, particularly 
when looking beyond the superrenormalizable theories, that they will not recover terms 
that are invisible in principle in the semi-convergent perturbative expansion, and in the 
context of these seminars I will therefore refer to them as "weakly" nonperturbative. 

On the other hand if one is interested in the more theoretical question of whether, 
or to what extent, the series determines the "true", exact function F^r, one must look 
for uniqueness properties justifying the application of what mathematicians call limiting 
procedures, or summability methods - the fine art of sensibly assigning a value to a 
divergent series. Such uniqueness results exist for the superrenormalizable theories and 




character notwithstanding, possesses the 



Mecanique Celeste [P0893| where he introduced the modern concept of semiconvergent 
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will be sketched in sect. 1.2; sect. 1.3 will then go on to discuss the most frequently 
used technique, the Borel-Laplace transform, of actually performing the reconstruction 
of r^r. In this context, sect. 1.3 will also take a look at the analytic-continuation step 
required for that transform - apparently a purely technical problem, but one that will 
guide us in more central, analogous problems arising further down the road in sect. 3.1. 
- All of the material of chapter 1 is " old" , and most of it in fact was old already when 
physicists became interested in it in the late 1960 's, dating back as it does to the work 
of mathematicians of the late 19th and early 20th centuries, but I believe it is of timeless 
interest and value, and deserves to be retaught on occasion. 

The situation is markedly different in the realistic field theories in D = 4 - in par- 
ticular, the asymptotically free gauge theories central to present-day particle physics - 
that are "strictly" or "marginally" renormalizable but not superrenormalizable (chapter 
2). There, the very definition of the theory involves a crucial new element, a divergent 
coupling renormalization. (More precisely, the tri- and quadrilinear interaction mono- 
mials in the Lagrangian of the theory, which in the superrenormalizable cases underwent 
at most finite rescalings, now demand their own ultraviolet-divergent renormalizations, 
but here we will continue to somewhat loosely refer to just the coupling renormalization 
that is at the core of these). If go (^) denotes the bare coupling originally appearing in 
the interaction terms of the Lagrangian, with mass dimension [^q] = e = ^ (4 — D) 
m. D = 4 — 2e dimensions, then its relation to the dimensionless and renormalized 
coupling g (^u) that parametrizes the functions ( p.3p must in these theories be written 
as 

gl{e) = g\^i)^?' Z^{g\^,),e) , (0.6) 

where the coupling-renormalization constant, Zq,, is ultraviolet divergent as e — > at 
any finite order in g"^ (n) except for the zeroth order, (0, e) = 1. This is in contrast 
to the superrenormalizable theories where is ultraviolet finite term by term and can 
always be transformed to 1 by a finite rescaling of the elementary fields (p (x) (and, 
possibly, of m (/i)). You have become accustomed (or so I presume) to the idea that the 
renormalization process is every bit as important in determining the dynamics of a QFT 
as the essentially classical Lagrangian, and you will therefore be prepared to hear that 
this new ingredient draws a fine but fateful line between the superrenormalizable and the 
strictly renormalizable field theories, a line signalling in fact a deep and fundamental 
divide. (I confess that doubts sometimes beset me as to whether the authors of the 
myriad scholarly papers studying superrenormalizable models in order to learn something 
for the realistic theory have always fully realized what they were up against.) As we 
shall recall in sect. 2.1, this new feature directly translates into the emergence of an 
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" RG-invariant mass scale", A, that is non-analytic at a = in such a way as to permit 
no power-series expansion in a at all, not even in the semiconvergent sense []GNV |. In 
such a theory, the very act of writing expansion (p.4| ) already implies that one has lost 
an important class of terms in the correlation or vertex functions, and although certain 
traces of the loss remain, it is essentially irretrievable. (The terms lost will become 
important as soon as one or more of the Lorentz-invariant momentum arguments ki ■ kj 
in our functions are not much larger than A^.) To specify the minimal information 
required for a reconstruction of such a function one then needs, in place of ( p.4D , a 
more detailed formal construct called a general resurgent symbol, or a resurgent double 
expansion, which still has a " perturbative direction" proceeding in powers of a, but now 
also features a "nonperturbative direction" providing for a systematic approximation of 
the new, strongly nonanalytic a dependence. As discussed in sect. 2.3, the reconstruction 
(through an extension of the Borel-Laplace process) then takes place in the framework 
of what might be called, in relation to the summable semi-convergent series, the next 
larger, mathematically well-delineated class of functions, known to mathematicians as 
the resurgent functions. In sect. 2.2 we sketch how an old acquaintance, the operator- 
product expansion (OPE) of Wilson |WIL72| in the case of an asymptotically free theory, 
can be recognized as a resurgent double expansion, and provides hints as to the special 
form this expansion will take in an asymptotically free situation. Section 2.4, in which 
I hope you will find something new, examines a simple model of a coupling-resurgent 
function, close enough in form to simple field-theory amplitudes to claim relevance for 
our subject, but still fully tractable by analytic means. This example will (I hope) 
convince you of the importance of treating a resurgent function as a whole: you will 
see how its analytic properties seem to sabotage the Laplace reconstruction from the 
perturbative expansion (|0.4D , but how in fact the reconstruction is entirely well defined 
when the strongly nonperturbative parts, whose separation is structurally averse to the 
function's resurgent character, are fully taken into account. 

The final chapter 3 is more special. It attempts to explain the basic ideas (while 
glossing over much technical detail) of a particular strategy for dealing with the one 
question to which the present mathematics of resurgent functions does not seem to con- 
tain a ready-made answer, and on which the physicist, for the time being, is therefore 
left to his own heuristic devices. This is the question of how to obtain some systematic 
formal representation of a resurgent amplitude (from which resurgent analysis proper 
may start) when the double formal series is not satisfactory for this purpose but needs a 
certain amount of analytic continuation to begin with. This need arises from the simple 
fact that the basic equations of motion generating the resurgent amplitudes of QFT are 
not, as in almost all cases studied in mathematical theory, ordinary differential, integro- 
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differential, or difference equations with respect to the coupling, but a rather complex, 
hierarchically coupled, and rcnormalizably divergent system of integral equations with 
respect to the momentum (or coordinate) arguments of correlation functions, equations 
into which the coupling enters only as a parameter. Actually, the emergence of functions 
resurgent with respect to a parameter - rather than the variable acted upon by a differ- 
ential operator - from differential equations has been studied in mathematics under the 
name of quantum resurgence, and properly speaking it is this kind of resurgence that 
will concern us below. But there is an extra difficulty here, arising from the fact that the 
nonperturbative coupling dependence enters as a dimensionful parameter, the mass scale 
A, and is therefore inextricably linked with the dependence on the momentum variables 
kf. in the integral equations, the loop integrals extend down to the smallest momenta 
\ki\, and in order to get started at all, we need to have something to insert for the inte- 
grands that is valid also at |/cj| ^ A, whereas resurgent symbols of the OPE variety only 
furnish asymptotic representations for 3> A. We are therefore going to need, from 
the very start, at least some partial continuation (namely, with respect to the A variable) 
of the symbols. Section 3.1 outlines one practical form of such a partial continuation, 
drawing on what we may have learned in sect. 1.3, which can be cast in the suggestive 
form of generalized Feynman rules, objects of zeroth order {p = 0) in the "perturbative 
direction" but providing a systematic continuation-through-resummation in the "non- 
perturbative direction". Section 3.2 sketches the special self-consistency mechanism by 
which the generalized Feynman rules can establish themselves in the hierarchical vertex 
equations, and mentions some nontrivial consequences of the fact that this mechanism 
is tied to the renormalizablc loop divergences of those equations. Section 3.3 is a brief 
outlook on some of the features of the coupling-resurgent iterative solution that can then 
in principle be generated from the generalized Feynman rules - an outlook certain to be 
incomplete, as many other such features probably remain to be explored. My goal here 
will have been attained if you come away from this with a realization that the two-step 
nature of this reconstruction process is not due to the shortcomings of the physicist 's 
heuristic patchwork, but is firmly rooted in the mathematical nature of resurgent func- 
tions. 

Despite their title, these talks are not lectures on mathematics. Although we are 
going to quote a few mathematical theorems without proofs and indulge in a certain 
amount of name-dropping, we will mostly proceed, as physicists do more often than not, 
by plausibility, formal calculation, and simple examples. The reason is lack of expertise, 
certainly not lack of enthusiasm, on the part of your speaker. I do feel obliged, how- 
ever, to draw your attention to the fact that there already exists a substantial body of 
mathematical research on the fascinating objects called resurgent functions. (I learned 
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what little I hope to have understood from selected chapters of the books CNP93 | and 
| STS96[| .) A large part of this research builds on the fundamental contribution of one 
person, the French mathematician Jean Ecalle, who in Les fonctions resurgentes (1981 
- 85), a three- volume work totalling almost 600 pages [ECA81|, first singled out resur- 
gent functions as a remarkable and unexplored phenomenon within classical complex 
analysis, studied comprehensively their structural aspects ~ the algebras and subalge- 
bras, auto- and homomorphisms, dual-space and fiber-bundle structures engendered by 
them - that are of premier interest to the mathematician, introduced, in the form of 
the so-called alien derivatives {derivations etrangeres) and of the associated alien dif- 
ferential equations, a whole array of powerful new tools and concepts for characterizing 
the peculiar continuation behavior of resurgent functions around singularities of their 
Borel transforms, and developed on this basis a new classification of ordinary differential 
equations of truly remarkable scope and depth, where large classes of equations related 
by quite general analytic mappings can now be characterized by a single " complete set of 
holomorphic invariants" . In these seminars, because of the limitations imposed by time 
and by my insufficient mathematical preparation, I am going to use only a few of the 
very simplest and most qualitative concepts from this beautiful edifice, but even these 
few will, I am confident, demonstrate to you one hallmark property of Ecalle 's theory: 
its power to structure, to clarify, and to provide larger perspective on, what previously 
may have looked like heaps of isolated and poorly understood results. 
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1 Superrenormalizable Theories: Summability 



1.1 Semiconvergent or Asymptotic Series 

We are accustomed to dealing with functions / (z) of a complex variable z analytic 
around z = and therefore having a power-series expansion in z within a circle of 
convergence around the origin. Semiconvergent, or asymptotic, expansions at the origin 
usually arise in situations where z = is not only a point of non-analyticity of / (z), but 
is in fact a member of a whole set of non-analyticity points - an endpoint of a branch cut, 
say, or an accumulation point of a sequence of poles, etc. - that reach up to the origin 
from somewhere in the complex plane, and the analyticity situation, not only globally 
but already locally around z = 0, is therefore "anisotropic" to begin with: the residual 
domain of analyticity, starting at z = 0, is usually concentrated along some direction, 
or sector of directions, in the z plane. In our context, where z is the coupling or a, the 
positive real axis is distinguished among possible directions of this kind: it represents 
the physical values of the coupling, and if a field-theory amplitude were singular or two- 
valued at points on this axis, the underlying QFT would quite simply be nonexistent as 
a physically acceptable theory. We will therefore always assume that our functions / (z) 
are analytic in at least some narrow region TZ along the positive real axis (mathematical 
theory, to be sure, deals with much more general situations). This, for example, may 
simply be a strip of constant width along the axis (Fig. 1.1 (a)), or a sector of opening 
angle 7 with apex at the origin (Fig. 1.1 (b)), or, to mention a seemingly exotic case 
that is nevertheless of relevance for QFT, a region bounded by circular arcs tangent to 
the positive real axis at 2; = (Fig. 1.1 (c)). The "sectorial" analyticity of Fig. 1.1 
(b) is the standard situation postulated in practically all mathematical treatments, and 
indeed often made part of the very definition of asymptotic series, but we should keep 
in mind that physics may force other kinds of domains on us. However, they all have in 
common that they do not allow for a finite circle of analyticity around the point z = 0. 

A formal power series in z, i. e., an object of the form 



with real or complex coefficients a„ and with nothing being said as yet about its possible 
value, is called a semiconvergent or asymptotic expansion for our function / (z) if for all 
z in the analyticity region TZ the error estimate 



00 




(1.1) 



n = 




(1.2) 
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Fig. 1.1 Examples of regions R of analyticity 
including the positive real axis 



is valid for all nonnegative integers N, with constants Cn + i such that the error bound 
eventually diverges for large everywhere in TZ: 

Civ + i|z|^ + i — > oo {N ^ oo) . (1.3) 

One expresses this situation by writing 

oo 

/(z)~^a„z", (1.4) 

n = 

with the ~ sign being read as " has the asymptotic expansion" . In contrast to the case 
of a convergent series, ( |1.2D does not guarantee that the error on its l.h.s. can be made 
arbitrarily small at all points in TZ by going to sufficiently large A^. The expansion may 
nevertheless by useful for obtaining estimates of / (z) at small z in TZ: for \z\ <C 1, the 
exponential decrease of with may, for the first few A^'s, win over the growth of 
the Ctv + i, with the divergence ( |1.3[ ) making itself felt only at larger A^'s, and if one 
halts the expansion " in time" , one may come away with an error that is tolerable for the 
purpose at hand. Such an error pattern, illustrated in Fig. 1.2, is behind all "pragmatic" 
applications of asymptotic expansions, where one strictly avoids asking questions about 



a possible meaning or value of the series (LI) to all orders. 
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This pragmatic handling of asymptotic series is facihtated by some simple formal 
properties: two formal series semiconvergent to functions /i (z) and /2 (z) in the same 
region TZ can be linearly combined, with complex coefficients Ci and and multiplied 
to give formal series that are semiconvergent to the corresponding linear-combination 
and product functions, (i fi + (2/2 and /i • /2, respectively. A series ( pTiP asymptotic 
to / (z) may be integrated termwise from z = out to some z G TZ to give a series 
asymptotic to a primitive of /, 

/ f{w)dw ~ Yl (1-5) 
« = i 

On the other hand the series produced by termwise differentiation, 

00 

(n + l)a„ + iz" , (1.6) 

n = 
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can be trusted to be semiconvergent to the derivative function /' (z) in TZ only if one 
knows in advance that this derivative exists, is continuous, and does possess an asymp- 
totic expansion. 

Can a function have more than one asymptotic expansion? It is easy to see that the 
answer is no |WIG77|: divide the estimate ( |1.2| ) by \z\^ so that one factor of |2:| remains 
on the r. h. s., and let z ^ from within TZ or, as we will put it in view of the special 
nature of our regions TZ, " from the right" . You get 



N 

f{z) - Z ^nZ^ 
n = 

7n 



as z 



+ 



(ah iV) , 



(1.7) 



a property often used as an alternative definition of semiconvergence. Spelling this out 
for successive A^'s one has 



fiz) - ao 
fiz) 



z~*0+ 



ao 



ai 







2-^0+ 



ao = f (0+) , 
«i = /' (0+) , 



ai 



02 



0-2 



f (0+) 
2! ' 



z-»0+ 



etc. That is, if a function f {z), analytic in a region TZ as specified above, has an 
asymptotic expansion at all, then that expansion is identical with its Taylor series "from 
the right" at z = 0: 



1 fd^f 



[n 



0,1, 2, .... 



(1.9) 



Moreover, by writing the Taylor expansion of / {z) with remainder term at z ^ TZ, one 
obtains the integral representation 



N 



/(-)-E 



a„, z 



n = 



N f(N+l) 



{tz) 



z 



(1.10) 



of the error term. If the derivatives l/*-^"^^-*] are all bounded throughout TZ by A^- 
dependent constants, then so is the bracketed integral, and the bound gives back the 
constants Ctv + i of (|l.2|). This point is important - for the bound of ( |l.2| ) to hold, it is 
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not sufficient to know the large-n behavior of the Taylor coefficients at z = 0+ (such 
as relation (|0.5|) ). We also need some boundedness properties of the function / and its 
derivatives away from zero, which the formal series alone does not supply. 

Importantly, the statement converse to this result is not true: a given formal series 
(1.1) can be a semiconvergent expansion for more than one function f (z) analytic in 
TZ, and in fact to infinitely many such functions. Indeed if /i (z) and /2 (z) differ by a 
function g (z) analytic in TZ and having all its derivatives vanishing as z ^ 0+, 

/2 (z) - fi {z) =g{z); (^) = (all n) , (1.11) 

then both have the same asymptotic expansion. Such a g (z) would vanish identically if 
it were to be analytic in a region around z = 0, but with our much weaker condition of 
"anisotropic" analyticity in a region TZ, it is quite possible for g (z) to be nonzero. From 
( |1.2| ) we can immediately deduce an alternative characterization of such functions, 

|^/(z)| < Z)jv + i|^r + ' (z G 7^, all iV) , (1.12) 

with constants D^. That is, they "vanish faster than any power of z" as z ^ 0+ in 
TZ. Of course, the two functions /i and /2 differing by g have different constants Cat 
(namely, C^^_^^ = C^^^-^ + + i) in their error estimates (|l.2D , so these constants are 
characteristic of the individual functions, not of their asymptotic expansion. 

In order to pin down more precisely this class of functions which "fail to show up" 
in an asymptotic expansion, and also to estimate the minimal error observed in Fig. 1.2, 
we need to adopt more specific hypotheses about the constants Cat in (|1.2|) . In view 
of what we said in the introduction about the large-p pattern (^]^) in QFT, it seems 
natural to assume from now on that 

C7^^) = c^j) \A\-^ N\ {N > 1) , (1.13) 

(although as already emphasized this does not follow from ( p.5| ) alone). Here only the 
positive constant C(j) depends on the individual function / (z), while B and A are to be 
the same for all functions having the same asymptotic expansion (this is not essential 
but simplifies calculations). This form clearly has the property (|1.3| ). It would obtain, 
for example, if in ( |1.1[)| ) we had the estimate 

l/^"^) (^)l < C(^) lAl-"" (Nlf , zGTZ. (1.14) 

The case B = is known as functions of Gevrey class 1 in the mathematical literature. 
With this assumption, a function g (z) as in (1.11) satisfies 



N 



\g{z)\ < C(,)iV^ ^ Nl (iV > 1) . (1.15) 
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We can find the lowest among this infinite sequence of bounds by minimizing the r. h. s. 
with respect to N at fixed \z\. For this purpose, we temporarily treat as a continuous 
variable, writing A^! = NT{N). The condition for ( 1.15 ) to be minimal at A'^ = 
then turns out to be 



In 



B + 1 







(1.16) 



where (x) is the logarithmic derivative of the gamma function. There is no closed-form 
solution to this equation, but we are primarily interested in small-2 behavior, \z\ <C \A\, 
where we expect N„i ^ 1, so we may use i/j (x) = In x + O (x~^) jGRY94 1 , the Stirling 
expansion applied to t/', to write 



In A^, 



+ 



B + 1 



The leading-order solution obviously is A'^^ 
expansion in -|- and obtain 



O 



-I 

2 I 



1 



= . (1.17) 
so we can write as a Laurent 



N„ 



{B + 1) + O 



(1.18) 



Plugging this into the r. h. s. of ( |1.15| ) while using for A^! once more the Stirling formula, 
we obtain the lowest bound on g (z) in the form of a small-z expansion. 



\9iz)\ < 



\9) 



'2tx\A 



B 



-B 



e-M[l + 0(|z|)] 



(1.19) 



where B' = B -\- h, and where the last factor in square brackets is a power series in 



l^l arising from the higher corrections in both eq. ( 1.18 ) and the Stirling formula. (Of 
course, unless the Nm of ( 1.181) happens to be an integer, the actual minimum will occur 
at [A^m] or [Nm] + 1, but in the case A',^ S> 1 of primary interest here, this leads only 
to a minor readjustment of the constants in ( 1.19| ) [^TS96| ). We will refer to functions 
satisfying a bound of this type as exponentially suppressed (remember that even [fB > 
0, the function exp -j^^ vanishes so fast as |2;| — > as to suppress arbitrary negative 
powers). Thus if a function g{z), analytic in TZ, has vanishing asymptotic expansion 



and in addition satisfies a bound of the special form ( 1.15 ), then our earlier statement 
"vanishes faster than any power" can be replaced by the more precise statement of 
exponential suppression. It is clear that the same result ( |1.19| ), with C(g) replaced by 
C(j) , also holds as a minimum-error estimate for the semiconvergent approximation of a 
function / (z) having a nonvanishing semiconvergent expansion, and in addition obeying 



the special condition (1.13) on its error-bounding constants Cn- 
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On the other hand if g (z) is known to be of the form 



g (z) = d^g) e-- [l + O {z)] {X > 0) (1.20) 

with real constants j4', B' , d(^g^, then all its derivatives contain the factor exp(— A' / z) 
and therefore vanish as 2: — > from the right, no matter what B' and the bracketed 
power series are, and so it at least has vanishing asymptotic expansion. That in addition 
it also falls into the class of eq. ( |1.15| ) can be proved generally only for the "sectorial" 
analyticity situation of Fig. 1.1 (b) with opening angle 7 < vr, where one may write 

_M , /I \ 

< e l-l with |A| = ^ cos ( - 7] , (1.21) 

but it may be true in individual other cases, and in any case the weaker statement is 
quite sufficient for our purposes. 

A very interesting observation in this context is that the square-bracketed power 
series in ( [1.19[) , arising as it does from the semiconvergent Stirling approximations for 
tjj (N) and T (iV), is itself semiconvergent, and that semiconvergent series are also admis- 
sible for the bracketed series in ( 1.20| ), without changing our conclusion. Thus within the 



class of functions characterized by (|1.13| ), the subclass of "escaping" functions that are 
left undetermined by a semiconvergent series are typically of the form of an exponential- 
suppression factor (possibly preceded by an arbitrary power z~^ that does no harm) 
times another power series which may itself be semiconvergent. This observation is one 
of the possible heuristic starting points for the more general theory discussed in chapter 
2. 

Of course, life would be much simpler if we could get rid of the undetermined part 
altogether when trying to reconstruct a function from its asymptotic expansion. In such 
cases the series is referred to as summable: there is then a unique function / (z) having 
this series as its small-z expansion. It is to such a situation that we turn next. 



1.2 Nevanlinna-Sokal Analyticity and Application to QFT 

A sufficient condition that has been widely invoked for unique reconstruction of a 
function from its asymptotic expansion, in both mathematics and physics contexts, is 
based on a theorem found by Phragmen and Lindelof already at the end of the 19th 
century. It again refers to the "sectorial" analyticity situation of Fig. 1.1 (b), but this 



15 



time with an opening angle 7 larger than tt (see Fig. 1.3): 

Phragmen-Lindelof Theorem: A function g (z) analytic 
in a sector S of radius R and angle 7 > tt with apex at the 
origin, and fulfilling the condition of exponential suppression 
(1.19) throughout S, must vanish identically. 



Im z 




I refer you to, for example, ref. STS96| ] for a proof of this theorem, which is based 



on an application of the maximum principle of complex-function theory. Note that a 
finite radius R of the sector is quite sufficient for this theorem; the essential point is the 
large angle of the analyticity sector. 

An immediate application of this theorem is to a function / (z) analytic in such 
a sector of wide opening and having an asymptotic expansion (|1.2[ ) there, with error- 
bounding constants Cn given by ( 1.13| ). We have seen in the last section that for such 



a function the "escaping" functions g{z), by which it differs from all other functions 
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having the same asymptotic expansion in the same analyticity domain, obey the bound 



( 1.19 ). We therefore have the 



Uniqueness Theorem T. If a function f (z) is analytic in 
a sector S of opening angle 7 > vr, and possesses in this sec- 
tor an asymptotic expansion (1.1) with error-hounding con- 
stants given by l \l.lS ), then it is the only function with this 
asymptotic expansion in S . 



(1.23) 



This apphcation was made in 1912 by G. N. Watson |WAT12], who combined it with 
a proof that one particular method of reconstructing the unique / (z) from the series 
- the Borel-Laplace process we are going to discuss in the next section - works under 
the same hypothesis. His result, known as Watson 's theorem., became the favorite 
summability provider for more than half a century, right through to the first summability 
proofs for (superrenormalizable) QFT 's in the mid-1970 's. It is a quirk of mathematical 
history that nobody seems to have been irritated by the fact that the Laplace-transform 
representation established by Watson - our eq. ( 1.451) below - did not recover / {z) in 
its full Phragmen-Lindelof analyticity domain of Fig. 1.3 used as input to the proof, but 
only in a much smaller subdomain - the hatched circle of Fig. 1.4. Didn't this look 
suspiciously as if analyticity in that smaller domain should have been sufficient in the 
first place? 



Sokal [ 5OK80 | in 1979 proved precisely this, and also redressed an injustice of 
history by pointing out that the result had in fact been given as early as 1918 by the 
Finnish mathematician F. Nevanlinna |NEV18|, though in a journal perhaps not readily 
accessible to many researchers. We may state a first segment of this result as our 

Uniqueness Theorem II: If a function f (z) is (i) analytic 
in the circular region Cr given by 



1 1 

Re- > - 

z R 



1 



R 



(R real and positive), and (ii) possesses in this region the 
semiconvergent expansion (1.1) with "factorial" error bound 
I{1.1^) , then it is the only function with this asymptotic ex- 
pansion in Cr. In particular, with R 00, analyticity in 
the half-plane \ arg z\ < ^ is a sufficient form of (i). 



(1.24) 
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Fig. 1.4 The circular region C/j in theorem ( 1.24| ). 



The second, more technical segment of the Nevanlinna-Sokal result, which again 



concerns the Borel-Laplace reconstruction, will be stated and discussed as theorem (1.41) 
below. 



Our two examples of sufficient summability conditions, theorems ( |1.23 ) and ( 1.24| ), 
create (when applicable) a situation which, if you think about it, is quite remarkable: 
under their hypotheses a divergent power series, initially a purely formal object, is never- 
theless found to be associated with, and thus in a sense to determine, a unique analytic 
function, in much the same way as a convergent expansion does. They are therefore 



sometimes [ 3TS96 referred to as quasi- analyticity principles. What makes them work 
is the potent mixture of analyticity and strong error bounding they postulate for the 
associated function in suitable domains - I only need to drop the catchword of Liou- 
ville's theorem in complex analysis in order to remind you of how powerful that kind of 
combination can be. Yet it is prudent to remember that both theorems spell out sujfi- 
cient, not necessary, conditions for unique reconstruction. Uniqueness may still hold in 
special cases not fitting into their assumptions, and section 2.4 will indeed treat you to 
an entirely nontrivial example of this. 

For the physicist interested in QFT, the most important applications of these the- 
orems undoubtedly have been the summability results for the coupling dependence of 
correlation functions in superrenormalizahle theories, beginning in 1975 with a proof 
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| EMS75 | for the cf)^ theory - a neutral scalar field with quartic self- interaction - in 
D = 2 Euclidean dimensions, and later extended along similar lines to ■0 (Yukawa) 
theory in D = 2 and, most significantly, to the superrenormalizable cj)^ theory in D = 3, 



the workhorse of the modern theory of phase transitions and critical exponents [ MAG77 1 . 
(It is widely believed that analogous summability properties also hold for other super- 
renormalizable theories, and in particular for abelian and non-abelian gauge theories in 
D = 2 and D = 3, though for some reason these do not seem to have been analyzed 
with the same degree of mathematical rigor.) These proofs involve a good deal of heavy 
mathematical technicality, and I cannot even begin to do them justice here, but let me 
try to give you at least a superficial glance at their main ideas. The starting point is 
always the Euclidean path-integral representation of correlation functions: taking the 
(f)'^ theories as an example, we have 

r, t 2^ /PA^[0]0(:ri)...<A(xjv)e-^n0] 



JV^l[ct>]e-SI\ 



where 

^2 



Si[4>] = Ij y" d'^xct^ixf [D = 2, 3) , (1.26) 

and where the Gaussian functional measure, P /i [</>], incorporates the bilinear part of the 
Euclidean action. If the integrals were ordinary multi-dimensional integrals, it would be 
straightforward to conclude they are analytic functions of g^ in at least the half-plane 
|arg(5f^)| < ^, i. e. Re{g'^) > 0, since the integrands are obviously analytic and 
the integrations strongly convergent in that region. In fact the integrations are infinite- 
dimensional and, for the axiomatic field theorist, therefore call for an exact definition 
of the Gaussian measures T) ^[(j)], and since ultraviolet renormalization in these theories 
can be dealt with at the level of one- and two-point functions, that definition will also 
have to account for renormalization. Conveniently, in the years before 1975, the mea- 
sures D fi [(p] for the prototype superrenormalizable theories had already been studied 
exhaustively, and it had been found that the (small finite number of) renormalization 
counterterms suggested by perturbation theory were in fact sufficient for such a construc- 
tion. Therefore these measures were already under good enough control mathematically 
so the analyticity argument could be carried over without difficulty. 



To establish the second pillar of summability, the factorial error bound of eq. ( 1.13| ) 



in the analyticity domain Re {g"^) > 0, the successive derivatives with respect to g^ 
of expression ( 1.25| ) were analyzed with the help of the cluster-expansion technique 



of statistical mechanics, which in particular allows division of the denominators into 
numerators - the "removal of vacuum bubbles" in physics language - to be made explicit. 
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and were estimated using suitable correlation-function inequalities also established by 
earlier constructive work. The result typically was an estimate of the form of eq. ( Ll^ ) 
above, 

gn 



GN{...;g') 



<Ci{C2r{nir {Reg'>0), (1.27) 



which when used in the remainder term ( 1.10| ) of the Taylor expansion around = 
yields the desired bound. 

An inconspicuous but nevertheless crucial precondition for the success of these 
proofs is that in eq. (1.26) we did not need to distinguish between a bare coupling 



go that appears in the interaction Lagrangian, and a renormalized coupling g that serves 
as the perturbative expansion parameter. Indeed as remarked above, in the superrenor- 
malizable theories the two differ at most by a finite rescaling, so one may expand directly 



in the g^ of eq. ( |1.26D without loss of generality, since a finite rescaling merely changes 
the meaning of the constant | A| in eq. ( 1.13| ). In the strictly renormalizable theories this 



simple connection is lost: between the Lagrangian go and the expansion parameter of the 
perturbation series, there is now the abyss of an ultraviolet-divergent renormalization 
( |0.6| ), which can no more be absorbed into a rigorous construction of a Gaussian measure 
P /i [(/>] . That this simplicity is lost is presumably the main reason for the continuing 
lack, after more than 25 years, of comparable results for realistic theories. 

The super renormalizable QFT 's are as close as one can come to a perturbation 
theorist 's paradise: in these theories the perturbation series, while divergent, never- 
theless "contains the whole truth" in that it has associated with it unique correlation 
functions without escaping remainders. Conversely, in such theories there is in principle 
nothing that cannot be represented by, and recovered from, a semiconvergent expansion 
in the coupling. The term " nonperturbative" in these models can at best have a rather 
restricted meaning: it can refer only to what in the introduction we dubbed weakly 
nonperturbative phenomena, whose calculation requires some all-orders resummation of 
their perturbation series. This makes both for the special charm of these theories and 
for their limited usefulness when trying to learn something for the realistic, marginallly 
renormalizable cases. 

Before moving on to the technicalities of the most widely used reconstruction proce- 
dure, it may be useful to emphasize that once you have verified some sufficient condition 
for uniqueness of / (z), it is in principle a secondary matter what method you use to 
get hold of that function. You may simply guess it if you are clever. Among more 
systematic methods, the choice depends largely on the pattern of divergence in your 
given series, but even for the factorially divergent series of type (^]^) one usually faces in 
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perturbative QFT, the Borel-Laplace procedure discussed below has no monopoly. To 
mention just one alternative, many functions having asymptotic expansions with this 



divergence pattern can be reconstructed by the method of Pade approximants [BAK70|: 
one approximates / (z) by a rational function with numerator and denominator degrees 
and respectively, 



f{z) - fM,N{z) = i + ,^, + ,^,2 + ...,^,^ ' (1-28) 



and adjusts the M + N + 1 coefficients of this function so as to make its first M + 
+ 1 Taylor coefficients at z = agree with the coefficients ao ...um + n of the given 
divergent series (|l.l| ). Of course, this method is easiest to apply if you have some advance 
information about the large-z behaviour of / (z) to guide you in the choice of the relative 
degree M — N. For example, if f (z) is known to fit into the general form 

oo 

f{z) = I dx^^ , a{x) > 0, (1.29) 

J z -\- X 



i. e. to admit a spectral representation with a cut on the negative real z axis (so that in 
fact we have the luxury of Phragmen-Lindelof analyticity in a " sector" with R = oo and 
angle just infinitesimally smaller than 2 7r), and if the coefficients in its formal small-z 
expansion 

oo 

an = {-rjdx^^, n = 0, 1, 2, ... (1.30) 



while forming a factorially divergent sequence, are individually finite, then the "diagonal" 
Pade approximants with = M, f(N,N){z), are known JBAK70 | to converge to / 



uniformly on compact domains in the cut plane. (If (|1.30) exists also for n = —1, so that 



/ (z) — > — a„i / z for |z| — > oo outside the cut, then the A^ = M + 1 approximants, 
of course, will converge better). Such functions are known as Stieltjes functions, and 
their Pade reconstruction has been used successfully to sum perturbation series in the 
(j)^ theory at D = 1 (the theory of the anharmonic oscillator) [|LOE69 |. However, even 



in the absence of such detailed information, the method is often applied heuristically 
because of its simplicity, generality, and relative robustness. It does have its drawbacks 
- when applied directly to a factorially divergent series, its convergence in low orders 
M + N may be erratic, and for this reason it is more frequently used as an analytic- 
continuation engine within the Borel-Laplace method, to which we turn in our next 
section. However, in principle it can stand on its own feet. 
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1.3 The Borel-Laplace Summation Process 



In the case of a summable formal power series 

oo 

ff{z) =^anz^ , (1.31) 



n = 



the most widely used method of recovering from the series the function / (z) uniquely 
associated with it is the venerable Borel-Laplace (BL) process [BOR28]. It requires some 



mild additional restriction on that function, but it is capable of handling the "factorial" 
coefficient growth 

an — > cA-T{n + B + I) (n ^ oo) , (1.32) 
which is essentially the same as that of ( p. 5]) since r(n + -B + l) (n + l)-^r(n + l) 



for n — > oo |GRY94|. It is therefore suitable for summing the perturbative expansions 



of superrenormalizable QFT. 

The first step is to map / (z) by an invertible integral transformation B onto a 
function b (s) of another complex variable s, the Borel transform of /: 

p+ioo 

b{s) =B[f-s):= [ dCe^'^fiC), (1.33) 

p—ioo 

where 

/(C) = J / (^) • (1-34) 

(This peculiar-looking prescription is adopted here to conform with the convention used 
below in eq. ( 1.45| ) for the Laplace transform.) The notation B [/; s) is used for a 



quantity that is both a functional of the function / and an ordinary function of the 
variable s. The path of integration in the C, plane, ReC, = p, is chosen such that all 
singularities of / (C) remain to the left of it. 
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Fig. 1.5 (a) Region of analyticity under Phragmen-Lindelof conditions 
and integration contour in the Borel transform. 



(b) Deformation of integration contour into r_ 



In the situation where uniqueness of / is guaranteed by its fulfrhing the conditions of 
theorem ( 1.23| ) , the domain of analyticity of / (C) is the image of the " wide-angle-sector" 
of Fig. 1.3 under the inversion C = 1 / z, and is therefore the shaded region in Fig. 1.5 
(a), while under the hypothesis of uniqueness theorem II, ( |1.24 ), it is simply the half- 
plane ReC > 



so in both cases we should choose 



— < p < oo . 



(1.35) 



From eq. ( 1.33| ), it is clear that the growth restriction we must impose on our function 
f (z) is modest: f (z) should be such that its associated f{C), eq. ( |1.34| ), possesses a 
Fourier transform with respect to the imaginary part a = ImC, of its argument. 



The transformation can next be applied, termwise, to the series ( 1.31 ), giving a series 
called the formal Borel transform, 



hf (s) = B iff] s) = an < 



n = 



1 



s{p+ioo) 



2Tri 



(1.36) 



s{p—ioo) 
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where rj = s p, and if the integrals in both ( |1.330 and ( 1.36 ) exist, this series is in turn 



semiconvergent to b(s). The integration path is Rerj = sp for real and positive s and 
is "tilted" for complex s, but since the only singularity of the integrand is now a pole at 
r/ = 0, it can always be deformed, by bending its ends over to the left, into the contour 
r_ of Fig. 1.5 (b). The resulting integral is Hankel's integral representation |GRY94| 
of the inverse gamma function, 

-L f dT^e^rj-^ = (1.37) 

at integer a = n + 1. (If this looks unfamiliar to you, take a integer and work out the 
residue of the integrand at r/ = 0.) Therefore 

oo 

a 



n 

n = 



We could have simply defined 6j by this series, but since we are going to need the concept 
of Borel transform for objects more general than the series ( |1.31j ), we have chosen to 
derive it here from the more general integral transform ( p.. 33 ). 



Thus we have manufactured a series with much faster-decreasing coefficients than 
( |1.31 ), and if the a„ coefficients grow no faster than ( |1.32| ) at large n, then 6/ is actu- 



ally something much better than a mere asymptotic expansion of 6(s), namely a truly 
convergent expansion: it is given by the generalized geometric series 

n = ^ ' 

plus a faster-convergent remainder series having a least the same convergence radius, 
and is therefore convergent in the circle 

\s\ < \A\ , (1.40) 

where it defines b (s) constructively. This is not yet all of b (s); we shall need to perform 
an analytic- continuation process A to continue b from the circle ( 1.4(]| ) into as much as 



possible of the s plane, and in particular along the real positive s axis if we are to recover 
the original / (z) uniquely at the real positive values of z in which we are primarily 
interested. It is no surprise that the conditions of theorem ( |l.24|) also guarantee that 
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this can be done - a statement known as the |S0K8C] 



Nevanlinna-Sokal theorem: If / (z) fulfills the conditions 
of theorem (1.24), then its Borel transform b{s) is analytic 
be yond the circle ( |1.40D in the striplike region \s — < 
\A\ on both sides of the real s axis R"*", as illustrated in 
Fig. 1.6. Moreover, b{s) is exponentially bounded there by 
const, exp^^^ , and / (z) is uniquely given by the Laplace 
integral ( |1.45| ) below for all z in its region of analyticity Re 



(1.41) 



z ^ R- 

We can see from ( |1.39| ) that (at least for B > —1) the sum of the bf series is singular 
at s = for ^ > 0, and at s = — 1^| for ^ < 0. Thus under the conditions of 

the Nevanlinna-Sokal theorem, the A in ( 1.32 ) can only be negative: the large-n terms 
of ( [1.31 ) must have alternating signs. We emphasize that ( |1.41 ) again gives sufficient 
conditions; the desirable behavior of b (s) on the positive s axis may still obtain in other 
cases with suitable modified hypotheses. 



Im s 




Fig. 1.6 Region of analyticity in the s plane of the Borel 
transform b (s) in the Nevanlinna-Sokal theorem 



With such behavior of b [s) secured, we can finally reconstruct from it the function 



/ {z) in its region of analyticity by inverting the transformation B of (1.33): the inversion 
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theorem in the theory of Laplace transforms [ GRY94 ] states that 



/(C) = / dse''^'b{s) . (1.42) 







[If you feel you need proof of this, write eq. (|1.33| ), with ( parameterized as p + ia, in 
the form of a Fourier transform, 



hoo 



e-''b{s) = ^ I dae'^'^fip + ia) , (1.43) 

— oo 

and evaluate the Fourier transform of the product 

e{s)[e-P^b{s)] , (1.44) 

by the Fourier convolution theorem plus contour closing, using [i (u; — ie)]~^ as the 
Fourier transform of the unit-step function ©(s)]. Therefore we recover our function 
/ (z) through what for our purposes we shall call the Laplace transform, 

oo 

f{z)=C[b;z):=-[dse-H{s), (1.45) 





for values of z such that the integral converges. Under the conditions of the Nevanlinna- 
Sokal theorem ( |1.41| ) above, this is obviously the case for Re ^ > i. e. in the entire 
original region of analyticity of /. Under the stronger Phragmen-Lindelof or Watson 
conditions, that region does not become any larger, but b (s) then has in addition a 
"sectorial" region of analyticity in the s plane, with opening 7 — vr, that may be useful 
for some purposes. For example, the integration path in the variable i = f , namely, the 
ray 

t = - ■ (0 00) , (1.46) 

can in this case be rotated back to the real axis without changing the value of the integral 
provided |argz| < ^ (7 — vr), and (|1.45|) is then equivalent to 



f{z) = J dte-'b{tz) , (1.47) 


which otherwise can be justified only for real, positive z. This form lends itself more 



readily to generalization to more than one variable | STS96 |, and incidentally is best 
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suited to performing the final check on our reconstruction procedure: from ( 1.47 ), the 
Taylor coefficients of / from the right follow by differentiations under the integral sign 
and use of (|L3|), 



1 



dt e 



nl \ ds ^ ^ 



s=0 



^ dte 



-t 



and since 



dt e'H^ = r(/3 + 1) 



(1.48) 



(1.49) 



this becomes equal to an, so ( 1.45| ) indeed gives the function / (z) that has the series 
(1.31) as its asymptotic expansion. 



A remark on conventions: our definition of the Laplace transform, eq. ( 1.45| ), is 
adapted to the study of semiconvergent expansions at small z, and has been chosen so 
as to agree with the convention used in the majority of the field-theory textbooks that 
discuss the subject [ WEI96[| | ZNJ89 ITZ80 |. It is also the convention used by Watson 
| WAT12| ] and by Sokal |SOK80| ]. Most of the mathematical literature on resurgent 
analysis (such as references | ECA81 |, [CNP93|, and parts of | STS96| ) studies expansions 



at a large argument C, and therefore needs to be translated to our context through z = ^ ; 
moreover it uses a different definition, namely eq. ( |l.42| ), for the Laplace transform. This 
makes for simpler appearance of a number of formulas (for example, the convolution 
theorem, line 10 in Table 1 below, gets rid of the extra z factor on its 1. h. s.) but on the 
other hand requires the separate introduction of a unit element, the delta distribution 
5 (s), for the algebra of formal Borel transforms hf (s), since ( 1.42 ) fails to map h = 1 
onto / = 1. To further complicate matters, a recent QFT textbook |SIE99| also uses this 



"mathematical" definition of the Laplace transform, while reference | STS96(| occasionally 
switches between small-argument and large-argument considerations. So be prepared to 
spend some time climbing around between the various conventions. 

To summarize, the BL reconstruction proceeds in three steps: 

(1) " Borelize" the given asymptotic series: 

B 



ff{z) (eq.(1.31)) 



hf{s) (eq.dOsI)) . 



(2) Continue 6/ to 6 analytically along the positive real s axis: 

A 



hf {s) 



h{s) 



(1.50) 



(1.51) 
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(3) Obtain / (z) from b (s) through Laplace transform: 

b{s) A f{z). (1.52) 

The importance of the second step, the analytic continuation A, can hardly be overem- 
phasized: without it we would make no progress at all. Indeed if we forget that step and 
simply plug expansion ( 1.3^ ) into the integral of ( [1.47 ) - i.e. if we use ( |1.38 ) beyond 



its radius of convergence - then obviously we do nothing but repeat the calculation of 
eqs. (|1.48|) and (|1.4g| ), and are led back to our original semiconvergent series (Il.31| ) 
without having come any closer to recovering the "true" function f (z). (Actually, one 
standard way in which semiconvergence arises is as a punishment for your having used 
a series in an integration outside its circle of convergence.) Needless to say the simplest 



method of continuing ( |1.38| ) analytically, namely to sum the series in closed form as we 



did in ( 1.39| ), is hardly ever applicable in realistic cases, where the a„ coefficients result 



from complicated perturbative calculations, and often are only available numerically for 
a limited number of orders n. The practical methods that have been used - particularly 
in studies of critical phenomena based on the superrenormalizable (j)^ theory in D = 3 



— are surveyed in chapter 37 of |ZNJ89|. Here I would like to mention just one method 
which, while no more the most refined one for the specific purpose of calculating universal 
quantities in phase-transition theory, is still widely used in statistical mechanics because 
of its simplicity and flexibility. This is again the method of Pade approximants, whose 



direct application to the original asymptotic series ( 1.31 ) we already mentioned in sect. 



1.1, but which is usually applied more profitably to the better-converging "borelized" 



series (1.38), where its convergence in the lower orders as a rule is less erratic. Thus one 



approximates 

b{s) « V,L)(«) = (1.53) 

with polynomials and Ql of orders K and L respectively, and with Ql (0) normalized 
to 1 without loss of generality. The K + L + 1 coefficients of this fraction are then 
determined either by matching the coefficients of ■s'^, s^, s^, ...s^~^^ in its expansion 
around s = to those of bf (s) (Pade-I method), or by matching the entire fraction at 
K + L + 1 points Si within the circle of convergence of 6/ to the numerically summed 
bf (Pade-II method). The L roots si of the denominator polynomial, 

L 

Ql (s) = l[{s- si), (1.54) 
1 = 1 

then approximate the singularities of b (s) - it is an advantage of the method over, 
e. g., conformal-mapping methods that one does not have to know the positions of these 
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singularities in advance. Nor, in principle, does one need to know much about the large- 
\s\ behavior of b (s), provided one is willing to work through the entire "Pade table" , i. e. 
the entire set of possible {K, L) pairs at given K + L + 1. U the true 6(s) bears some 
resemblance to ( 1.39| ) at large s, then the approximants with K < L will converge best, 



and one may write a partial-fraction decomposition 

L 



b(K,L) (s) = > , ri = (1.55) 

which allows the Laplace transform ( |1.45| ) to be carried out analytically according to 
line 11 of table 1 below: 

f{z) ^ f^K,L) = -- E rie-TEi (^) . (1.56) 

1 = 



If / is Borel summable and // has real coefficients, the si must all be off the positive 
real s axis, and thus come either at negative real s or in complex-conjugate pairs, so all 
the Ei functions in ( 1.56| ) are unambiguous. 



For the last step (1.52) of the reconstruction scheme, it is often convenient to have 



a little toolkit of basic BL transforms available, so Table 1 offers such a list. All entries 
in this table are straightforward to verify by elementary calculation. (Many lines carry 
some clauses or restrictions that are usually obvious and not mentioned explicitly.) The 
first line, for u = n an integer, repeats what we have already used to derive (1.38), but 



it also shows that it is just as easy to "Borelize" more general formal series with some 
noninteger power in front: 



oo 

Z' 



n = / n = 

The fourth line is interesting because it states that B and C commute with the infinites- 
imal generator of dilatations, so that / (A z) and b (A s) form a transform pair if / (z) 
and b (s) do. The tenth line expresses the well-known convolution theorem of Laplace- 
transform theory, with the convolution bi * 62 of two functions being defined by the r. 
h. s. 

Lines 11 and 12 feature a simple higher transcendental function studied already by 
Euler in the 18th century, the exponential-integral function Ei (z) defined [ |GRY94 ] by 
the integral 

z 

Eiiz) = du— (z / 0) . (1.58) 
./ u 
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Table 1: Some Borel-Laplace transform pairs 





f{z) = jr[b; z) 


b{s) = B[f; s) 


1. 


z"" {y > -1) 


Tiu + l) 


2. 


^•/(^) 


s 

[ dt b{t) 
J 




3. 






4. 


df{z) 

d - 


db{s) 

r/.s 


5. 


df{z) 
dz 




6. 




s ■ b{s) 


7. 


z 

^ J dw f{w) 




s 

\ j dt b{t) 




8. 


/(^) -/(o) 

z 


db(s) 
ds 


9. 


1 — UJ Z \1 — CO z J 


e'^'bis) 


10. 


z ■ fi{z) ■ h{z) 


s 

(bi * 62) (s) = J dtbi (t) 62 (s - t) 



11. 


Z \ Z ) 


1 

S - So 


12. 


z \ z J 


^ e (s - so) 


13. 




9(s - so)6(s - So) 
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The equivalent spectral representation (w = z — u) 



Ei{z) = e"" / dw (1.59) 

J z - w 



shows that the function is analytic in the z plane except for a cut along the positive real 
axis, where its discontinuity is 

Ei{x + iO) - Ei{x - iO) = -27rie{x) (x real) . (1.60) 

The cut is thus logarithmic, as shown explicitly by the alternative representation 

z 

Ei{z) = 7i? + In(-z) + / du "^ ~ ^ , (1.61) 

J u 



in which the integral is an entire function. Formal expansion of the denominator in 
(1.59) and integration beyond the convergence circle results in the divergent series 

e-" °° 1 

EH^) = -T.P'-^^ (1.62) 

p = 



a _a „ (a 
- e 

z \z 



so that 

oo 

Ei (") = ^ (a-Pp!)zP (1.63) 

p = 

is a simple example of a factorially divergent series of type (p.5|); line 11 of Table 1 just 
describes the BL reconstruction from this series. Line 12, on the other hand, follows 
from the application of line 3 to the Laplace-transform representation 



Er (-^ 









^ dse 2 







e(s-a)ln(^-) , (1.64) 

b 

where a is real and positive. This representation embodies the modern-day significance 
of the Ei function as the simplest among a family of elementary building blocks for 
resurgent functions, the so-called resurgence monomials | ECA81 |. 

The last two lines of Table 1 actually go beyond the framework we discussed so far: by 
considering Borel-transform functions with jumps, they remind us of the fact that for 
the Laplace integral ( |1.45| ) to be well-defined, it is quite sufficient that b (s) should only 
be piecewise continuous on the real s axis, rather than fully analytic in a region along 
it. In this they point ahead to chapter 2, where we will exploit this greater generality. 
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As a matter of etiquette, we must acknowledge that the three-step scheme of (|1.5C| ) 
to ( |1.52|) should actually be preceded by a zeroth step: 



(0) Check conditions for summability, i.e., uniqueness of f (z) . (1.65) 



But you know that real life isn't like that - outside the arcadian realm of the 
superrenormalizable theories, physicists hardly ever have enough information to do such 
a check. So as a rule they will proceed heuristically, apply the BL machinery to the 
limited information they have, and expect to be punished by some singularities on the 
real Borel axis should they have committed some invisible offense. If you feel bad about 
this state of things, you may find solace in the words of the great Poincare, who clearly 
had a heart for the poor physicists when (on the limits of " pragmatism" in dealing with 
asymptotic series) he wrote | P0893(| : "Les astronomes ne les connaissent pas toujours 
d'une facon bien precise, mais ils les franchissent rarement, ... d'ailleurs leur instinct 
les guide et, s'il les trompait, le controle de T observation les avertirait promptement de 
leur erreur" . 
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2 Strictly Renormalizable Theories: Resurgence 



2.1 Coupling Renormalization and Lambda Scale 

When passing from the superrenormahzable theories discussed in sect. 1.2 to a "re- 
ahstic", strictly renormaUzable theory, such as {cI)^)d = 4:, ( Yukawa) £> = 4 , or a nonabehan 
gauge theory in D = 4, we encounter no changes in the form of the corresponding 
Lagrangians, or in the symmetries displayed by these Lagrangians. Also, the renormal- 
izations related to the bilinear parts of the Lagrangians - field renormalizations and, if 
bare-mass terms are present, mass renormalizations - remain ultraviolet divergent, as 
they were already in the superrenormahzable cases, the only change being that they are 
now divergent in all orders of perturbation theory, rather than in a few low orders only. 
The only fundamentally new element in the definition of the theory is the divergent cou- 
pling renormalization (^]^). It becomes necessary because the three-point and four-point 
vertex functions that evolve through quantum effects from the trilinear and quadrilin- 
ear Lagrangian interaction terms now become superficially divergent, which amounts to 
saying that the bare coupling becomes dimensionless as -D — > 4. 

The information about this new constitutive element, which will turn out to have a 
strong effect on all the correlation functions of the theory, is originally contained in the 
coupling-renormalization constant Za of ( |0.6| ), which perturbatively takes the form of a 
Laurent expansion in the dimensional regulator e = (4 — D) / 2, 

{g^,e) = 1 + Ci (5') I + (9') ^ + - • (2-1) 

Each (^n is in turn a (semiconvergent) series in starting at order {g'^)^, since ultraviolet 
divergence of type e""" only comes from diagrams with at least n loops: 

/ n \2n / g \2n+2 

Cn{g^) = Cn,n[j^) +Cn,n + l(^) ••. (n > 1) . (2.2) 

The set of all coefficients Cn,fc (iT' > 1, k > n) perturbatively represents the complete 
information about the coupling renormalization. However there is a much more compact 
and convenient way of encoding the same information, suggested by renormalization- 
group (RG) analysis, namely, in terms of the RG beta function, 

/3(5(/^)) =M%^. (2.3) 
a/i 

This function is UV finite (the renormalized coupling g (p) being UV finite by definition), 
and its dependence on the regulator is therefore much simpler. Since this is an important 
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point in all that follows, let me briefly repeat how these two encodings are connected 



| TH073| . We can in principle determine /3 (g) from Za by applying the operation / 
on both sides of (|0.6|) and using the fact that ^q, a bare quantity, knows nothing about 
the sliding scale ji, so d^g I = 0- We find 



-£5 



l+g^^ (In Z«) 



(2.4) 



Upon inserting expansion ( p.l[ ) this becomes 

^3 r 



/?(5) = -eg + /Cl + Y 



C2-CiCi-5'(Ci)^ 



+ o 



(2.5) 



Here the primes denote differentiation by g'^ . Since [3 {g) must be finite as e — > we 
conclude, first, that all terms with powers of ^ on the r. h. s. must vanish separately. 
This gives a sequence of conditions on the Cn which generally determine C,^ for n > 2 in 
terms of lower C,'s. For example, the ^ term displayed in ( |2.5D gives 



(2.6) 



These relations are interesting in themselves; they may be exploited e. g. to derive a 
succession of partially resummed forms for the perturbative expansion of Z^, but since 
such forms are more easily derived from the closed-form result of eq. ( |2.12| ) below, we do 
not dwell on them here, noting only that all the with n > 2 are in principle determined 
by the function (i {g'^) alone, the residue of Za at e = 0. Thus (^) is actually a highly 
redundant form of conveying the information on coupling renormalization. 



(2.7) 



Second, comparison of the finite parts in eq. ( |2.5D now gives 



P{9, e) 



-eg + g 



dg^ 



and therefore a (semiconvergent) perturbative expansion 



|9(9, e) 



+ 



(2. 



= -(fc + l)Ci,fc + i (A: = 0, 1, 2, ...) . (2.9) 

You will remember the particular significance of (5q, the leading beta- function coefficient 
at e = - the theory is asymptotically free, i. e. its coupling decreases like an inverse 
logarithm at large scales /i, if and only if (3q > 0. Also, let me recall without comment 
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that the first two coefficients, /3o and Pi, are independent of the renormahzation scheme 
i? in a large class of schemes that includes the dimensional schemes | WEI96 |. 

Having expressed, in eq. (^), the beta function in terms of the residue at e = 
of Za, we may now invert the procedure and express in terms of the beta function. 
From (2.4), we have 



d 



In Zn 



1 



where we abbreviated (g / 47r)2 = K, and where 

x{k) -- 



(4vr)2c;(52 



/3o + /3l K + /?2 ^2 + ... 



(2.10) 



(2.11) 



contains the information of the beta function at e = 0. Integrating with respect to 
under the initial condition Za (0, e) = 1 then gives the integral representation [ rH073 | 



(3/4^)2 



Za {g^, e) = exp 



(2.12) 



This exact expression shows the well-known double face of renormahzation constants: 
when expanding the integrand " perturbatively" , i. e. in powers of k with /5o/e as 
a leading term, and integrating termwise, you get back the perturbation expansion of 



( |2.1| ) and (2.2) in which every term is ultraviolet divergent as e — > - though with all 
coefficients now expressed in terms of those of (3 (g) or (5^), and with the relations for 
the higher Cn 's, like ( |2.6| ), therefore built in: 



Za (g, e) = 1 - Po 



_9_\^ 1 
4ir 



(2.13) 



On the other hand if you ask about the behavior as e of the full, resummed 
expression ( 2.12| ), you find exp (—00) = ! To see exactly how the function vanishes 
requires a minor detour in order to get rid of the lower-limit divergence of the k integral 
at e = 0: write 



1 



d 



exIk) 1 + ex (k) dK 



— ln[K + ex{K)] 



(2.14) 



where the prime is now a differentiation with respect to k, and perform a partial inte- 
gration. You get 

Za {9, e) = F X 
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X exp < e 



1 + ex' (z^) 



1 + ex' (0) 



(2.15) 



[l + ex {^')? 



\n[K + exiK )] 



The remaining integral is finite at e 
only in that limit; so 



0; the exponent has terms of order e and e In e 



[l + 0(e, eln e) {g\ii))] 



(2.16) 



is an exact statement. (The writing emphasizes that the additional terms still depend on 
(n) in a complicated way). This will not surprise lattice theorists, since when inserted 
into (^]^) it says that in the removal-of-regulator limit of a nonperturbative calculation, 
the bare coupling should be tuned to zero according to 



9o (^) 



-^^.[i + 0{e,eln e) {g' (fx))] , 



(2.17) 



4vr J l3o 

and for them this is entirely familiar if they replace e with something proportional to 
their regulator a, the lattice constant. 

There is yet a third way of encoding the information on coupling renormalization, 
which in a sense is the most physical since it relates to the representation of dimensionful 
observables. In the strictly renormalizable theories the coupling becomes dimensionless 
as -D — > 4, and if the theory has no explicit mass terms in its Lagrangian (or, as in QCD 
with only light quarks, masses that are orders of magnitude too small to explain anything 
observed), then there seems to be at first sight no dimensionful quantity available to set 
the scale for dimensionful observables. For example, a pure Yang-Mills theory can still 
generate nonzero glueball masses, or nonzero scattering lengths (low-energy scattering 
amplitudes) between its gauge bosons. The solution to the paradox [GNV| is that 
the theory creates its own mass scale "spontaneously" through quantum effects - and 
essentially the same quantum effects (i. e., loops) that are responsible for the necessity of 
coupling renormalization. Since observable quantities are to be expressed as multiples of 
powers of this scale, it cannot simply be the arbitrary sliding scale ^, nor can it depend 
on arbitrary gauge-fixing parameters. In an otherwise massless theory, this scale, called 
A, must therefore depend on and on g (/x) in such a way as to be an RG invariant, i. 
e. it must be a solution of 



A(g (/x), fi) 



A(<?, /u) 







(2.18) 
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the simplest RG equation conceivable. Mass dimension requires A to be of the form 
■ Se {g (/u)), with the dimensionless function iS^ (g) fulfilling 



1 + Pig, e)^ 
dg 



£e (g) = 0, 



(2.19) 



so separation of variables gives the representation 



^eigifJ-), ^) = At exp < 



dg 



(2.20) 
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with an integration constant gi, a trivial RG invariant. A simpler but somewhat more 
formal way of seeing the existence of A ist to write eq. (|2.3| ) as 



1 



dg 



fJ- ' Pig, e) 

and integrate between any two values /^i, ij,2 of the sliding scale, which gives 



me 



-^{g (Mi),e) 



(2.21) 



(2.22) 



with ^' (g, e) denoting a primitive of 1 / /? {g, e); this says that the combination ( p. 20 ) is 
/i- independent. Thus a divergent coupling renormalization, as expressed in the RG beta 
function (2/7), directly translates into the "spontaneous" emergence of a A scale. 

The integration constant gi, like the higher beta- function coefficients (^]^) with 
k > 2, depends on the renormalization scheme adopted, and in fact is part of its 
specification. Changing gi to g'l while keeping the higher /3fc's unchanged (a special 
change of scheme) gives another scale A' differing from A by a finite dimensionless 
factor: 

h{g'i) - ^(51) • 



A. 



exp 



(2.23) 



Thus the A scale, while an invariant within one renormalization scheme, is still different 
for different schemes, with the A's for two 
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schemes differing by a pure-number factor. In that sense it is not yet on the same level of 
theoretical nobility as a physical mass, which must in addition be scheme-independent, 
but it is a necessary and important intermediary between the arbitrary and the mea- 
surable masses and lengths. 

From the point of view of coupling dependence, the essential point about the A scale 
emerges in the limit e — > 0: by expanding perturbatively. 

Pig, 0) Po9^ [ \(3oJ Uvr 

where the terms O (y^) again contain the scheme-dependent higher beta coefficients, we 
have 



(2.24) 



An 



fie 



9 

47r 



/3l 
2^ 



1 



Pi - P0P2 



a 

At: 



(2.25) 



2/3q / x,„. 

where the (semiconvergent) series in curly brackets depends on the renormalization 
scheme i?, and where for the class of schemes we use we have chosen gi such that, 
in the notation of ( 2.22| ), 



*(5i,0) 



In 



Po 



(4vr)2_ 



(2.26) 



Thus A at I? = 4 is the perfect example on an "escaping" function of the form of 
( 1.20 ), a quantity that will never show up in a perturbation expansion in (7^. In making 



this identification, it is clearly essential that Pq > 0, i. e. that the theory should he 
asymptotically free. In asymptotically non-free theories the scale Aq, rather than being 
exponentially suppressed, would be exploding as g"^ 0-|-, and none of the things we 
are going to discuss in the following would be applicable. 

Moreover, one should not lose sight of the fact that when talking about "behavior 
as g^ — > 0-|-" here, we are actually talking about a rather fictitious, unphysical limit - 
a limit in which we imagine 5^ being tuned to zero while ji is being kept constant. Only 
in this purely formal limit can we speak of Aq being "exponentially suppressed" - in the 
physical world there is no way of getting Aq to vanish, since the only parameter there 
which we can "tune" at will is /i, and Aq by construction just does not change with fi. 

Now eq. (|2.25| ) would not create an essentially new situation if the only way Aq 
occurred in correlation functions were the trivial one: one may of course invert eq. ( 2.2C| ) 
to express g in terms of A^ / /x, obtaining at e = the well-known expression ( [[WEIOq] ) 



47r 



Po Ao 



2Pi 
Pi 



In Ac 



+ O 



1 In Ao 



^0 



X2 
-^0 



(2.27) 
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where 



An 



(2.28) 



One may then parameterize the perturbation series in terms of the inverse of this loga- 
rithm rather than in terms of a, but this wih not introduce any new information beyond 
the perturbative one, and this remains true even if we go one step further and invoke 
the RG scahng properties of vertex functions to rewrite the expansion in terms of a 
couphng a that "runs" with Q^, some representative combination of the function's mo- 
menta, rather than with /i^ [ WEI96[| . (This would be the same function as ( 2.27| ) with 



simply replaced by Q"^ .) Since 



a(g2) 



a 



1 + 



a 



power series in - — (3q In — ^ 
47r V /i^ 



(2.29) 



(as follows by combining eqs. ( 2.2C| ) at scales fx and \/Q^), this amounts to a resumma- 
tion of the perturbation series that is useful for many purposes, but again creates no new 
content beyond that series. But then if your QFT text mentions no other application of 
the lambda scale than eq. ( 2.27| ), you should protest - they are hiding the best from you! 
For example, an s-wave scattering amplitude in a massless theory at low center-of-mass 
energies ^/s would look like 



his) 



1 



doAo + + + 



(2.30) 



with do Aq the inverse scattering length, and with the dimensionless coefficients di pos- 
sibly having some perturbative, power-series dependence on the coupling ( |2.2S| ). Or a 
bound-state pole in a four-point amplitude T of a massless theory would appear as 



T{p , P, p; a) 



a^{p /Aq; a) ■ a$(p/Ao; a) 
P2 + cAg 



-|- [regular terms] , 



(2.31) 



(where the notation indicates that the factorizing residue, in contrast to the purely non- 
perturbative denominator, is a mixed perturbative-nonperturbative quantitiy beginning 
with at least an overall factor of a^). In other words, the truly interesting quantities 
are powers and inverse powers, polynomials and inverse polynomials, series and inverse 



series in Aq, and in these, unlike in (2.27), Aq appears in its full, non-analytic nastiness, 
eq. ( ^.25 ). We will, for the purpose of these seminars, refer to such quantities as strongly 
nonanalytic at a = or strongly nonperturbative, as opposed to mere resummations of 
a perturbation series, where Aq only appears in inverse logarithms. 
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2.2 Operator- Product Expansion and Lambda 



Our digression into the Operator-Product Expansion (OPE) will be brief and highly 
selective, serving only the narrow purpose of reminding you of a long established way 



I WIL72 | of asymptotically representing the strongly nonperturbative coupling depen- 



dence. For a systematic development of the OPE, refer to | WEI96 |; this reference also 



contains, in the special context of renormalons, a discussion equivalent to the one below. 



For a terse account of Zimmermann 's perturbative renormalization |ZIM70| of the OPE, 



you may consult | KUG97 ]. 



Again we only look at one simple example that will exhibit all features of interest 
— the Euclidean vacuum expectation of a single product of two generic elementary fields 
with naive mass dimension one at variable separation x, 

D{x) = (OlvP Qx) ^ {-]^x^ |0) , (2.32) 

in the context of an asymptotically free theory {(3q > 0) without Lagrangian masses 
{m = mo = 0). (General OPE treatments would consider correlation functions with 
this pair of operators inserted into an arbitrary product if {yi) if {y2) ...if{yn) of fields 
at fixed, distinct, nonzero points yi). Loosely stated, the OPE is the postulate that the 
field product permits a small-x expansion of the form 

/I \ / 1 \ 

^ ( 2^) ("2^) = 5Z 5Z Wn,i{x)O2n,i{0) , (2.33) 
\ ^ ^ / n = i = l 

where the 02n,i (z) are the full set of local operators of mass dimension 2n, elementary 
or composite, that bear the conserved quantum numbers of the field product. Thus at 
the core of the OPE is a completeness hypothesis which, to my knowledge, has nowhere 
been proven with full mathematical rigor but is nevertheless highly plausible: the full 
set of local operators at x forms a basis for the expansion of operator products clustering 
around x. 

These operators are enumerated here according to their naive mass dimensions 
2n = 0, 2, 4, ... and to a counting label i ranging over the finite number / (n) of such 
operators with given n, the enumeration beginning at n = with / (0) = 1 and the 
dimensionless unit operator, 

Oo,i = ll. (2.34) 
To make up for the increasing mass dimensions, the c-number functions Wn,i{x), the 
"Wilson coefficients", consist of powers (x^)"~^ times a perturbatively calculable series 
in the coupling that contains powers of In x^. 
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More precisely, what is meant |WIL72] by ( |2.33| ) is that this expansion should hold 
in the sense of "weak operator convergence" , i. e., not with respect to any operator norm, 
but only when both sides are sandwiched between normalizable states and (^'2] of 
the quantized-field system. The two-point function ( |2.32| ) is then the simplest sandwich 
of this kind, where = [^'2) = |0), the physical vacuum state. In a renormalizable 
theory, the resulting matrix elements, the "vacuum condensates" 

(0|O2n,i(0)|0) (n = 0, 1, 2, ... , i = 1, 2, .../(n)) , (2.35) 

should finally appear in a renormalized form, referring to a definite renormalization 
scheme R and to a scale fi within that scheme, but since the detailed way the renormaliza- 



tion works in an expansion like ( 2.33 ) has some aspects of interest for our later discussion, 
we will begin by considering a not-yet-completely-renormalized version, parameterized 
by Q (fi) but still at a positive value of the dimensional regulator e = (4 — D) / 2. When 
taking Fourier transforms, such as 

Wn,i{x) = I d'''^'xWn,i{k';a;^)e"' '' , (2.36) 

the above mentioned structure of the Wilson coefficients is mirrored in the momentum- 
space structure 

~ ( ^ j ^ perturbatively calculable series in a} , (2.37) 

with the coefficients of the series containing, in general, both ^ ultraviolet divergences 
and the typical logarithmic modifications In (fc^ / ^u^). Therefore when passing to the 
two-point vertex function r2 = —D~^ the expansion will look like 

- r2 (fc ; a; e) = A; Vn,i (a; In ; - 1 r^- , (2.38) 

n=Oi=l 

with perturbatively calculable coefficients, 

(^2 -j^ \ °° / ^2 \ 

a; In — ; - ) = ^ / In — ; - ) a (/i)^ . (2.39) 

As for the vacuum condensates, which at this stage have their own ultraviolet diver- 
gences at e — > 0, let us start by considering the expansion purely within the context 
of perturbation theory. There, we have no known way of " renormalizing" composite 
operators; i. e. redefining them so as to have finite matrix elements, other than defining 
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them as free-field normal products :02n,i-, but since those all have zero expectation 
values in the perturbative vacuum, all condensate terms with n > 1 in ( ^.38 ) would 
then " renormalize" to zero. Thus we conclude that the n = term, involving the trivial 
operator ( p. 34 ), is the perturbative portion of D~^, 



rf a; e) = k^Voi { a; In 



(2.40) 



and, since the zeroth-order perturbative part is simply k'^ 



Vol a; In 



fc2 1 



1 + O (a) . 



(2.41) 



The situation would be different, of course, for the more general correlation-with- 
insertion functions on which most discussions of the OPE focus: there one encounters 
matrix elements of the free-field normal product between non-vacuum states, which in 
general have nonvanishing perturbative limits. Pure vacuum expectations like ( 2.32| ), 
which in a sense are more vulnerable to perturbative mutilation, are particularly suited 
to our purpose because they compel us to make a crucial point: as originally envisaged 
by its inventor |WIL72], expansion ( 2.3^ ) is to be understood here in a more general, 
nonperturbative sense, where vacuum condensates for 2n > 2 can be nonzero even 
though the OPE by itself offers no immediate recipe for calculating them in a given 
renormalization scheme. This indeed is the standpoint taken in the successful semi- 



empirical applications of the OPE through QCD sum rules |S VZ79| , another wide and 



interesting field on which these seminars will be unable to enter [NAR89|. 



Even without calculating the condensates we may, to some extent, discuss their 
renormalization formally to see more precisely what is "incalculable" about them. Under 
renormalization, composite operators of the same mass dimension are known to mix 
linearly: 

(0|O2n;.(0)|0)„,. = («' -J {02n,,)R,, • (2.42) 

The renormalization matrices Z^'^ ") are of the general form 



1-3 



(2n) 



a. 



6^J +0{a) , 



(2.43) 



with all terms O (a) ultraviolet divergent, while the renormalized condensates {02n,j) R,^ 
are by definition ultraviolet finite but in general dependent on jjL and on the the scheme 
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R. The right-hand side of (|2.38| ) becomes 

fc2 



Vo,i ( a; In ^ 



1 \ oo l{n) ( l{n) 1 



For such an expansion to remain within the realm of renormahzable QFT in the presently 
accepted sense of the term, one of the minimal prerequisites is that nonlocal ultraviolet 
divergences of the form 

which would require counterterms of the same form in the Lagrangian, should be absent. 
Therefore the curly-bracketed quantities in ( |2.44| ) must all be finite for at least n > 2; the 
Wilson coefficients must be rendered finite by the composite-operator renormalizations 
separately for each n: 

J2 ^"'^ ("5 e) ^S"^ I) = ("5 1^ ^) ' fi^it^ as e ^ , (2.46) 

where, because of ( |2.43| ), 

= K,,(a = 0) + 0(a) . (2.47) 

This "internal" renormalization at each n > 2 is indeed what happens in the purely 
perturbative calculations alluded to before, for quantities where our condensates are 
replaced by more general matrix elements having nonzero perturbative limits. It is also 
implicit (though, curiously, never discussed explicitly) in all cases where one- or two-loop 
Wilson coefficients have been calculated under nonperturbative assumptions allowing for 
nonvanishing vacuum condensates jNARSgl] |RRY85| ] iLAV92|| . This is analo gous to the 



well-known result of renormalization theory that nonlocal divergences of type 

W (>n (2.48) 

cancel separately in each order of the perturbative expansion. Here, since we lack the 
time to delve into explicit examples, we will simply accept this general statement as 
plausible. (Incidentally, one of the uses of the techniques described in chapter 3 may be 
to provide a more precise framework for such demonstrations.) 

In the n = 1 term, the counterterm required for an ultraviolet divergence would 
be constant in momentum, of the form of a mass renormalization, and therefore not 
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excluded by renormalizability, but it turns out to be really necessary only in the case of 
m 7^ 0, a nonzero Lagrangian mass; for a massless theory, where this term represents a 
contribution to dynamical-mass formation not visible in the Lagrangian, the finiteness 



of eq. (2.46) will persist at n = 1. 



On the other hand for the perturbative term with n = 0, associated with the triv- 
ial operator ( |2.34| ), there is no internal renormalization to bring ultraviolet finiteness 
- the unit operator "doesn't volunteer" to tame its Wilson coefficient, and that coef- 
ficient alone - that is, the perturbative series - must therefore be supplied with the 
usual "external" renormalization through a field-renormalization factor Z3. Thus the 
general scheme, in the presence of the nonperturbative OPE terms, is no more that of 
multiplication with one overall Z3 for the entire vertex function, but rather 



finite "'-^ finite 



I have been dwelling on this simple but important point because the OPE literature, 
vast as it is, does not seem to state it clearly, and because it means that in the presence 
of the nonperturbative terms, multiplicative renormalizability in the usual simple sense 
is no more a property of QFT amplitudes. We will encounter this peculiar "two-track" 
renormalization pattern again, in the context of a more dynamical approach, in ch. 3. 
below. 



In eq. ( 2.49 ) we have not yet separated the perturbative and strongly nonpertur- 
bative coupling dependences completely. To achieve such a separation, we need to note 
that in a massless theory, the renormalized condensates at e = must be of the form 

{02n,j)R;^ = (« (^))fi (A^)^" , (2.50) 

since Aq is then the only invariant mass scale available for this purpose. The coefficients 
Cnj, power series in a {^) and in general semiconvergent, 

oo 

Cn,j{a) = Y.^n\o^' . (2-51) 

g = 

are what gives the vacuum condensates their residual ^ dependence within the scheme 
R - "residual" because the terms without explicit powers of a, 

i^lr , (2.52) 

are RG invariants within R. It is these terms, with no perturbative power of a but with 
the strongly nonanalytic a dependence in Aq still present, that form the "incalculable" 
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cores of the condensates in the context of the OPE. One of the exhilarating experiences 
to which the formahsm of chapter 3 wih invite you is that these cores become calculable 
there as multiples of powers of Aq in a given (dimensional) scheme, with coefficients 
^^n \ given by a systematic sequence of approximations, which in turn can be determined 
self-consistently from the equations of motion of the basic vertex functions. 

Putting things together, and introducing modified Wilson coefficients 

l{n) 



(a) 



(2.53) 



that refer only to the mass dimension 2n but no more to individual condensates, we 
then have a more transparent form of the OPE, 



where 



(R) 



n = l 



In 



7 



a 



In —^-j a 



1 



A;2 



Zs a; - Vb,i In a; - 



1 



(2.54) 



(2.55) 



e J \ fJ' £ , 
denotes the renormalized perturbation series. Again, the terms without powers of a but 
with Ao kept nonzero. 



ia(o) 

R 



1 + 



n = l 



(0, 0) 



A2 

A;2 



(2.56) 



with 



l{n) 



0) 



.(0) (\-^ 

-n, j K-^^O 



2\n 



(2.57) 



^(^) (0, 0) {Kir = Vn,, (« 

i = i 

form an RG- invariant "core" for the two-point function, to which the Un power series 
add /i-dependent, perturbative corrections. 



As rearranged in the form of eq. ( p.54| ), the OPE tells us three things. First, 
the genuinely nonperturbative coupling dependence of a correlation function can be 
represented by a series in integer powers of Aq, each power being accompanied by its own 
series of perturbation corrections with the typical logarithmic momentum dependences. 
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This structure has been confirmed in ah cases where an OPE could be studied to some 



degree of exphcitness in renormahzable model theories with a lambda scale [ DAV82 1 



| BBK98 |, and may nowadays be regarded as well established, but it was not always so 
clear in the early years when attempts to calculate condensates by instanton methods 
seemed to indicate the necessity of terms with noninteger powers of A |TH07£]. Work 
by G. Miinster [MUN82| later established that these noninteger powers were really an 
artefact of the "dilute-instanton-gas" approximation used in most instanton calculations, 
and that by taking account of the denseness of the instanton gas one is led back to integer 
powers of A (see t 'Hooft 's introduction for the 1994 reissue of his 1980 lectures, p. 457 
ofref. |rH094| ]). 



Second, since Aq is dimensionful as opposed to the dimensionless a, such an expan- 
sion is automatically a large-momentum expansion in powers of Aq / fe^ . This automatic 
linkage between coupling and momentum dependences is much tighter and more direct 
than in the perturbative series, where rearrangement in terms of a momentum-dependent 
running coupling a (Q^) is a useful option for physics but not a mathematical necessity. 

Third, we have learned that each of the strongly nonperturbative {n > 1) terms 
must establish itself in an ultraviolet-finite manner - a requirement which obviously 
places nontrivial restrictions on any dynamical theory seeking to calculate those terms. 

The OPE, since it does not by itself produce a prescription to calculate its own 
vacuum condensates, remains a deep structural statement but does not represent a com- 
pletely known starting point for an ab initio dynamical theory. It is nevertheless - even 
apart from the successful applications already mentioned - of great value for clarifying 
the coupling dependence of amplitudes in the strictly renormahzable theories with a A 
scale. To see this dependence more explicitly, we may once more rewrite our result (|2^ ) 
by using the structure ( 2.25| ) of Aq. Introducing the perturbative expansions 



1 



/3? - Pop2 



4^) 



2n 



R 



P = o V / 



(2.58) 



we then have for our two-point function the form 

°° / ..2\ " , 



n = 



4tv \ n 

e 1^0°' 



/So 



a 



n OO 

p = 



(2.59) 
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a form that stands directly at the entrance door to the theory of resurgent functions. 



2.3 Resurgent Symbols and Functions 

You win remember that at the end of sect. 1.1 we made a point of the fact that 
the "escaping" functions of a semiconvergent expansion typically appear again as small- 
z expansions, eq. ( |r20D , whose power-series factors {1 + O {z)\ in general are again 



semiconvergent. This observation may be taken as a starting point for the following 
heuristic consideration. 

Assume we are dealing with a function / (z) having an asymptotic expansion 

oo 

/(z) ~ ao,pZ^ (2.60) 

p = 

but not sufficient analyticity to be reconstructed uniquely from it as in chapter 1. We 
may take some representative /o {z) from the class of functions having this expansion 
and write 

/ (z) = /o [z) + 50 {A , (2.61) 

with go (z) the "escaping" remainder. Assume further that / (z) has bounding constants 
of the form of eq. (1.13) so the bound of eq. ( |1.19D holds, and that go (z) actually has 



terms, as e. g. in ( |1.20| ), that exhaust the bound: 

oo 

go{z) = e-^F,{z) ; Fi{z) ~ z'^^ ^ ai^pZ^ {Resi > 0) . (2.62) 

p = 

Of course, may have other terms, of faster exponential or even super-exponential 
decrease, but the series in the -Fi factor, an object as formal as the original series ( 2.6C| ), 



has enough freedom of interpretation to accommodate all of these. In the simplest case, 
that series will turn out to have a finite radius of convergence and thus to determine 
uniquely an analytic function Fi (z), so that the ~ sign in ( 2.62| ) turns into an equals 



sign within that convergence radius. In this case, the combination of eqs. (2.61) and 
(2.62) is a complete and exact representation of our function / (z). 



But a much more interesting case, and the one we concentrate on in all that follows, 
is for the situation of eq. ( |2.60| ) to repeat itself: assume Fi (z) again has a small-z 
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expansion that is only semiconvergent, and "escaping" terms led by an exponentially 
suppressed one: 

Fi (z) = h {z) + gi {z) , gi (z) = e-"? F2 (z) {Rer2 > 0) , (2.63) 

where /i is some representative of the class of functions having the second member of 
(2.62) as asymptotic expansion. Again gi (z) may have other terms, of faster exponen- 
tial or even super-exponential decrease, but we are again free to stuff these into the 
"escaping" terms left undetermined by the formal series expansion 



p = 

With that understanding we may then write 

/ (z) = fo [z) + e-"^ /i [z) + e-"? F2 {z) 



is2 



si + r2j 



(2.64) 



(2.65) 



Imagine now that this pattern keeps repeating itself, or "resurging", so that at every 
stage the escaping functions turn out to be led by an exponentially suppressed term. 
Then upon putting conventionally 



So 



Bn 



, 



we arrive at a formal representation 

fiz) 



n = 



where 



and where 



/„ (z) : = a representative of the class of functions 
having asymptotic expansion ^ CLn^pZ^ , 

p = 







Re sq < Re Si < Re S2 < Re S3 



(2.66) 



(2.67) 



(2.68) 



(2.69) 



The representation may terminate if at some n the formal series for the function F„ (z) 
turns out to have a finite radius of convergence and determines uniquely an analytic 
function Fn (z), which is then equal to fn (z) without a nonanalytic remainder. (This 
case actually occurs [STS9C] in some of the classical transcendental functions obtained 
as solutions to linear, second-order differential equations.) 
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However if the decomposition process does not stop and new, nonzero, suppression 
constants r„ (or s„ = Sn-i + Tn) keep emerging indefinitely, then the representation 



( 2.67 ) remains highly formal in that it lacks proper definition at two fronts. First, 
we have said nothing as yet about the convergence or divergence of the n series, nor 
about its resummation in the case of divergence. Second, while it appears plausible 



that by going far enough in the expansion (2.67) we will eventually account for all the 



exponentially suppressed terms, we have completely lost track of the possible terms with 
super-exponential suppression, which at every step we swept under the rug of the next 
Fn + i{z). Thus at every finite stage of the decomposition, such as eq. ( |2.65| ), they 
were still accounted for by the simple device of the last term having an F„, instead of 
fn, that still had room for such terms ~ just as in an ordinary semiconvergent series 
all suppressed terms were still present in the exact remainder term of the finite Taylor 
expansion ( l.lOl) . We lost them the moment we discarded this subtle difference and 



jumped ahead to the infinite formal series 

(^M) " just as in earlier times we lost all 



terms with faster-than-powers suppression in jumping from ( l.lCj ) to the formal series 
( |1.1| ). The first and second deficiencies are of course related ~ by carrying the analogy 
one step further, we can guess that the super-exponentially vanishing terms, if any, must 
be hiding behind a divergence of the n series, and that accounting for such terms must 
be tantamount to giving a prescription for resumming that series. 

Making bad things worse, we now forget even about our list of specific choices of 
representatives /„ (z) for the sums of the various semiconvergent p series, and simply 
write 



= 



'^n,pZn , (2.70) 



a representation still more formal than ( ^.67 ) in that it now cries out for answers, not 



only to the above-mentioned and related two questions, but also to the question of 
the meaning of each of the p sums. This formal double expansion is called a general 



resurgent symbol in the mathematical literature [ECA81], and the successive terms of its 
n sum - each carrying a definite exponential-suppression factor exp (— s„ / z), except for 
the n = term - are referred to as the elementary resurgent symbols. Because of its 
highly formal character, it would at first seem like an exotic way of making one's life 
complicated; yet it is precisely such a double expansion that arises in several areas 
of both mathematics and physics - the perturbation theory of quantum-mechanical 
energy levels, the JWKB semi-classical expansion, the thermodynamic functions of the 
ideal Fermi gas, the perturbation expansions of celestial mechanics, or, to mention a 



purely mathematical problem that forms the subject of vol. II of ref. | ECA81 |, the 
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classification of germs of analytic functions around the origin through a study of their 
iterations. Most important for us, a glimpse at eq. ( 2.59| ) tells us that the OPE, in a 



strictly renormalizable and asymptotically free theory, is also of this form: in a QFT 
with divergent coupling renormalization ( p.6| ) and therefore a coupling-nonanalytic mass 
scale ( p. 20 ), the operator-product expansion is a resurgent- symbol formal representation 



with respect to the coupling of the correlation functions. This representation in terms of 
a formal double expansion is the new paradigm that replaces the formal perturbation 
series (p.4[) in the realm of the strictly renormalizable theories. 



You should realize that eq. ( 2.7C| ) is a highly peculiar kind of expansion: it looks as 



if one and the same functional dependence on a single variable z were to be expanded in 
two sets of basis functions simultaneously. If somebody would suggest this to you in the 
context of convergent expansions, you would surely protest against the uncontrollable 
double counting. What makes the resurgent-symbol expansion possible is the specific, 
and matching, incompleteness in both types of "basis functions": the expansion can 
never account for an exp (— const. / z), and exp (— const. / z), being an infinite Laurent 
series with only negative powers of z, can never account for positive powers z^. 

The possibilities of reconstructing analytic functions from such resurgent symbols 
form an important part of the theory of resurgent functions. As a preliminary definition 
- which will not satisfy the mathematician, but conveys enough of the agenda to be 
helpful to the consumer - we may thus formulate: 



A resurgent function is an analytic function that can be 
reconstructed from its formal small-z expansion if the latter 
is given at least "to exponential accuracy", i. e. in the form 
of a resurgent symbol (2.70). 



(2.71) 



In dealing with such functions, the physicist must revise his or her ingrained habit of 
neglecting exponentially suppressed terms in the small-z expansion "because they are 
smaller than any power". While this may be a useful approximation numerically at 
sufficiently small z and under favorable circumstances (namely, if the leading coefficients 
ctn,0) n > 1, are not too large), it is completely irrelevant when it comes to knowing 
the "true" function: in reconstructing the latter, one either keeps the exponentially 
suppressed terms or gets nonsense, tertium non datur. 

Our above heuristic argument emphasized the peculiar property that gave resurgent 
functions their name - the fact that semiconvergent expansions with the special proper- 
ties of "factorial" error bounding in N, and of the ensuing exponential error suppression 
in z, "resurge" at every level of refinement of the small-z expansion. Yet it does not lead 
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to the most general type of resurgent symbol. The more systematic mathematical ap- 



proach [ECA81] proceeds by firmly placing the problem from the outset in the framework 
of the most powerful and versatile reconstruction method, the Borel-Laplace process of 
sect. 1.3. The starting point there is a reformulation of the basic Laplace-transform 
representation (|1.45| ): one introduces a modified Borel-transform function B{s), called 
the major of / as opposed to the minor b (s), such that B (s) has as its sole singularity 
a branch cut along the positive real s axis with discontinuity equal to the minor b{s): 

B{s + ie) - B{s - ie) = e(s)6(s) . (2.72) 

For a h (s) sufficiently decreasing as s ^ +00, such a function may simply be constructed 
as a spectral integral, 

00 , 

B{s) = -^ f ds'i^, (2.73) 

ZTTl J S — S 



which by virtue of the well-known distribution relation 

1 1 



{s + ie) s' — {s — ie) 



2 in 6(3 - s) (2.74) 



fulfills condition ( [2.721 ). (Such functions have occasionally been used in physics work 
too [ DAV82| ].) If b{s') does not decrease fast enough, or even increases, so that the 



simple spectral integral fails to exist, one may use various "subtracted" or "divided" 
forms of the spectral representation: for example, if under Nevanlinna-Sokal conditions 
b (s) turns out to be unfriendly enough to exhaust its exponential-growth restriction 
\b{s)\ < const, exp at s > 0, one may set up a spectral representation for the 



function exp ■ B (s) with R < R, which retains all analytic properties of B 

since the exponential is an entire function, and obtain 

/ 

00 _ s , 

5e(.) = e(^) J- fds'^^^, (2.75) 

ZTTl J S — S 



where the integral now converges. Of course, this super convergence device may also be 
applied in cases where it is not strictly necessary; for example, when the pure-and-simple 
Cauchy integral (2.73) already exists, we obviously have 



B,{s) = B{s) + [ ds'—, -bis) , (2.76) 

ZTTl I S — S 
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where the second term is an entire function since it is now analytic both off and on 
the positive s axis. This example illustrates the general truth that an infinite family 
of functions with the prescribed right-hand discontinuity ( 2.72| ) can be constructed, 



which all differ from each other by the addition of entire functions. The notion of 
major therefore really stands for such a whole equivalence class of functions, known 
to mathematicians as a hyperfunction. In a QFT context we are usually dealing with 
minors h{s) sufficiently decreasing as s — > +00 - the toy example of eq. ( |2.191 ) below 



being typical in that respect - and therefore will usually be able to get along with 
the "canonical" representatives of eq. ( p. 731 ) without really using this generality of 
hyperfunctions. 

For an ordinary asymptotic series summable to a function / {z) by the Borel-Laplace 
procedure, the Laplace representation ( 1.45| ), by virtue of eq. (|2.72| ), can now be rewrit- 
ten in terms of the major B (s) as 

/(z) = i j dse-iB{s) , (2.77) 

the integral being taken along the contour on both sides of the positive s axis as 
shown in Fig. 2.1. In view of the analyticity of B (s), this contour can in turn be opened 
up, as indicated in Fig. 2.1, to turn into the contour Fq parallel to the imaginary s axis. 
Thus as an alternative to ( |1.45| ) we have an integral transform along the path Fq - that 
is, in essence a Fourier transform - of a major function B (s) whose sole singularity is a 
branch point at sq = 0. 
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Fig. 2.1 Integration contour r_(_ in major-based 

BL transform and its deformation into Fq 
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Fig. 2.2 Integration contours in BL integral for 

resurgent function with support points Sn 
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It is this alternative Laplace integral that forms a convenient point of departure for 
the generalization to resurgent functions. Assume we are now given a function B (s) 
which instead of a single branch point at sq = has a denumerable set of singularity 
points Sn, either branch points or poles, in the right half plane Re s > 0, possibly includ- 
ing one at So = 0, and ordered as in eq. ( |2.69| ). See Fig. 2.2. We label B (s), somewhat 
summarily, by this set {sq, si, S2, ■■■} of singularities, which in the mathematical liter- 
ature is referred to as the support of a resurgent function or symbol (or both). Assume 
that B (s) is such that the Laplace integral with contour Fq, 

f(^) = 1 I ^^e-tS{,,,,,,,,,,„}(s) , (2.78) 
To 

exists in at least some "anisotropic", right-hand z-plane domain. This is the basic 
representation of a general resurgent function, provided one further condition is fulfilled: 
in order to obtain a formula suitable for calculating / (z) at small, positive z, we must 
be able to bend the contour Fq back to the right as in Fig. 2.2 so as to turn it into 
a union of horizontal back-and- forth contours F„_|_, n = 0, 1, 2, ... , each slung around 
one of the singularity points Sn, so that 

/(^) = E- / ^^e-f i?|,„,,,,,,....|(s) . (2.79) 

n ^/ 

J- n + 

By using integration variables s' = s — s„ in the individual terms, this can be written 



as 




f{z) = ye-^{- / dse-'^bnis)} (2.80) 







with a set of discontinuity functions, or Borel minors, given by 

bn{s') = B {sn + s' + iO) - B {sn + s - iO) , n = 0, 1, 2, ... , (2.81) 
where s' is real and positive. The representation ( |2.80| ) is precisely of the form encoun- 



tered in eq. ( 2.67 ) above, but with the /„ (z) factors now being given a more concrete 
meaning through the bracketed Laplace transforms. The final step towards a small-z 
expansion is then to perform Taylor expansions at z = OH — in general semiconvergent 
- of these transforms. For the more common types of singularities of the major B (s), 
this step is straightforward: in particular, if B (s), locally around each of the s„, is of 
the logarithmically singular form 

1 



. In {sn - s) 

ZITT 



4>n (s - Sn) + {regular at s„} , (2.82) 
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with factors (pn analytic at zero argument and thus having expansions 



oo 
q = 



with finite convergence radii, then since 

ln(-s' - iO) - ln(-s' + iO) 
at real s' , the n-th Borel minor simply becomes 



bn (s ) = (pn {s) = ^ bn, 



-2i7rG(s ) 



{s > 0) , 



(2.83) 



(2.84) 



(2.85) 



g = 



and the small-z expansion is obtained in the form of eq. ( |2.7C| ) above, with all Bn 's 
equal to zero. A function f {z) with such an " all- logarithmic" Borel major is called a 
"resurgent function with simple singularities" , and is the kind of function most frequently 
considered in the mathematical literature on resurgence. In the QFT context, where in 
view of the form ( 2.59| ) of the Operator-Product Expansion we expect an expansion with 
non-integer Bn exponents, only a minor extension of this concept is necessary as long as 
< Bn < 1: we assume B (s) to be locally of the form ( 2.82| ) only near s = sq = 0, 
but of the form 



2 i sin {Bn vr) 



■ <Pn {s 



+ {regular at Sn} {n = 1, 2, 



(2.86) 



near the other s„, with Bn non-integer. (The numerator term of minus unity serves no 
purpose other than to make the limit Bn exist and agree with (|2^)). Then 



bn [s 



(0 < 5„ < 1) , 



(2.87) 



and the general form of eq. ( ^.70 ) is obtained. We see how a resurgent symbol arises 
as a small-z expansion "to exponential accuracy" of a resurgent function given by eq. 
(llTSl) . 

It has been essential, in deriving the resurgent symbol, that the original contour Tq 
could be deformed into the set of r„_|_ 's without changing the value of the integral. For 
this it is necessary, of course, that the growth (if any) of the major B(s), or at least 
of its discontinuities (2.81), as Res — > +oo be sufficiently limited to allow the Laplace 
transforms of eq. ( 2.80 ) to converge. But there is also a more subtle precondition, which 
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in the presence of an infinite set of singularities Sn as allowed by our definition is no 
more trivial: we tacitly assumed that the arcs-at-infinity that are invisible in Fig. 2.2, 
but necessary to turn the union of (minus) the old Fq and the new r„+ 's into a closed 
total contour and thus to establish the equivalence, can always be supplied even if the 
singularities occur at arbitrarily large distances |,s„|. In other words the set while 
possibly spreading out to infinity in the right half of the s plane, must not be allowed 
to accumulate anywhere in such a way as to completely obstruct, for some direction or 
sector of directions, analytic continuation of B (s) to infinity along paths between them. 
The mathematical literature on resurgence calls such functions endlessly continuable. 
We are thus in a position to state the mathematical definition of resurgent functions: 



A resurgent function is an analytic function / (z) permitting 
a Laplace-transform representation (2.78) / ( 2.79D with an 
endlessly continuable Borel major B (s). 



(2.88) 



It is clear that this definition is extremely general ~ it places almost no restrictions on 
the detailed nature of the singularities of the major at the points Sn, or on the detailed 
structure of the analytic factors cpn (s — Sn) associated with them, except that they 
should not destroy the possibility of a Laplace-transform representation, and that they 
should remain sufficiently isolated to permit endless continuation in between them. It 
is no more necessary, as our introductory heuristic argument seemed to suggest, that 
each Vn = Sn — Sn-i be identifiable as the exponential-suppression constant in the 
series of the preceding, (n — l)-th term. That argument, leading to eqs. (|2.67] ) and 



(2.70), really dealt with a special case, where one not only encountered the same type of 



singularity, namely (2.86), for all n, but also the same type of semiconvergent series Fn (z) 



(namely, with factorial error bounding), and therefore the same type of convergence in 
the analytic factors (pn, so that in fact the complete local behavior of the major near 
any s„ gets revived near all the others. It is this special but also especially interesting 
case which gave the entire class of resurgent functions their name. In Ecalle's words 



I ECA81 1 , these are functions " possessing a built-in mechanism that counteracts the usual 
"erosion" of local properties under analytic continuation and provides instead for their 
spotwise resurgence" . 

Even within this special class, you should still be prepared for all sorts of nasty little 
complications, when actually trying to reconstruct a resurgent function through, say, the 
series ( p. 801 ) of partial Laplace transforms. Our definitions required the expansions ( 2.83| ) 



of the analytic factors 0„ (v) to have finite convergence radii, but this does not prevent 
them from having the singularities that limit their convergence circles sit maliciously 
at real positive v. If these singularities are poles, then the corresponding partial Borel 
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minors, obtained by analytic continuation of (|2.87] ), have poles at real positive s' , and 
the associated Laplace integrals in (|2.80| ) have no unique meanings any more - they 
depend on the way one evades the pole when integrating, and therefore become two- 
valued functions. The complete resurgent function is then, in general, multi-valued. 

The mathematical definition ( 2.8^ ), apart from the weak restriction of endless con- 



tinuability, allows for very general forms of the support {s„}, the set of singularities 



of the Borel major, which via eq. (2.80) become the set of exponential-suppression 
constants. In many cases - particularly when considering nonlinear equations where 
products of a function with itself or its derivatives must be dealt with - one may want to 
move within sets of resurgent symbols that are closed under formal series multiplication 
(form a multiplicative algebra), and mathematical theory therefore prefers to deal with 
support sets that form additive semigroups. Of special importance among these is one 
particular kind of support, where all s„ sit equidistantly on the nonnegative real axis, 

Sn = nsi ; si real and > ; n = 0, 1, 2, ... , (2.89) 

and which not only plays a prominent role in resurgence theory (in particular, it is ubiq- 
uitous in Ecalle's treatment of resurgent solutions to nonlinear differential equations) 
but, as it happens, is also present in the resurgent symbols of OPE type relevant to 
QFT coupling dependence. These, as illustrated by our two-point example of eq. ( 2.591 ) 
above, feature the constants 

47r 

Sn = nsi si = — > ; n = 0, 1, 2, ... (2.90) 







Bn = nBi; Bi = ^ ; n = 0, 1, 2, ... (2.91) 

Po 



In this case, the resurgent symbol ( 2.70| ) assumes the special form of a double formal 
power series: 



n 



f{a) ^ Yl ^" [a-'^'e-^) ^n,poP . (2.92) 

ri = p = 

The occurrence of an n-th power of a dimensionless "counting parameter" u for the n 
summation, where in our case 

u = ^ , (2.93) 

is a typical feature of such symbols. Since all Sn are now on the real s axis, the repre- 
sentation of eq. (2.79) can alternatively be written as a single Laplace integral. 



f {a) = - I ds e-^b{s) , (2.94) 



a 
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with a combined Borel minor given, according to line 13 of Table 1, by 



6(s) = ^ G(s - nsi)bn{s - nsi) . (2.95) 



n = 



It is no more analytic along the positive s axis, but has finite discontinuities in the 
support points. For the individual minors 6„ (s), the class of functions ( 2.87] ) considered 



up to now must be enlarged here to include distribution- valued terms, since by eq. 
( 2.91[ ), becomes larger than one as n increases so the first line of Table 1 can no more be 



used. (Recall typical values such as i?i = || or i?i = ^ for SU (3) gauge theory with 3 
or 6 fermion species respectively, for which i?„ > 1 at n > 3 or n > 4.) For the terms 
m eq. (|]92|) where p — nBi = —k is a negative integer, the Borel transform is a delta 



distribution 5 (s) when k = 1 and a derivative of a delta distribution, 5^^ {s — Sn), 
when k > 2: 

a-^e"^ = £[<5('=-^)(s - s„); a) (A: = 1, 2, 3, .... ) . (2.96) 
This is verified immediately by using the definition 

ds5^'''Hs)9is) = H'~-'[^)^ ^ (2.97) 

in the Laplace integral ( |2.94| ). For the terms with noninteger p — nBi = —f3 < — 1, 
one even needs to introduce fractional derivatives of delta distributions, functionals that 
perform the minor miracle of the mapping 

oo 

^ j ds e-f (s - Sn) = z-^e"^ (/3 non-integer) . (2.98) 



I can hear your oophs - that is quite some machinery to deal with an object as simple 
as a""^^! Indeed the combined minor of eq. ( ^.95| ), having now not only the finite 
step-function jumps but also all kinds of distribution-valued terms at the support points 
s = Sn, becomes complicated in a way that appears somewhat out of proportion for 
the relatively simple structure of the symbol ( 2.92| ). While the Borel-Laplace framework 



remains essential for the definition and basic theory of resurgence, it may therefore well 
turn out to be a clumsy tool when it comes to the practical recovery of a resurgent / {z) 
from its symbol, and it may become legitimate, for that purpose, to look around for 
simpler alternatives - we emphasized earlier that Borelization, for all its power, has no 
monopoly. 
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On the other hand these observations show clearly that distribution- valued terms of 
delta and derivative-of-delta type are natural ingredients for Borel transforms of resur- 
gent functions - as indeed we might have guessed already from line 3 our Table 1. They 
do not create new problems, in the sense that they do not infringe on either the ex- 
istence or the univaluedness of the Laplace transform. This is the main reason why 
in the following we will not be concerned with the so-called "instanton singularities" 
of Borel transforms that arise in the semiclassical treatment of QFT's with instanton 
solutions. These, as discussed by t'Hooft [ TH079 |, are simply integer-order derivatives 
of 5 {s — Sn), which produce discrete additions to the Laplace transform by the pattern 
of eq. ( |2.96D above. 

When does a resurgent symbol ( 2.7C| ), or the more special double-power-series form 
(2.92), determine uniquely an analytic function f (a) whose small-a expansion it rep- 
resents? (In this case the symbol is again referred to as summable.) Here comes a 
big disappointment: to the best of my knowledge, no sufficient conditions are presently 
known that would come in as handy for a realistic, asymptotically free QFT as the 
Nevanlinna-Sokal criterion ( 1.24| ) does for the superrenormalizable theories. Such con- 
ditions, if they were to be of practical value, would as a minimum have to allow for the 
special analyticity situation in the complex a plane governed by the so-called t 'Hooft 
smg ularities |TH07g|] , at which we had a glance in Fig. 1.1 c) and whose origins we 
will explain in more detail in sect. 2.4 below - singularities that preclude even the min- 
imalistic assumption of a small-angle "sectorial" analyticity of / (a) that is routinely 
made in the mathematical treatments of resurgence. The modest contribution of sect. 
2.4 will be to prove, by explicit counterexample, that these singularities by themselves 
need not sabotage unique reconstruction. But a compact and convenient positive crite- 
rion accommodating them seems to be lacking. It is clear that if all the Borel minors 



obtained by analytic continuation of the series (|2.87| ) have sufficient smoothness and 
limited growth to make the partial Laplace transforms in ( |2.8[l| ) well-defined, and if the 
n series in ( 2.80 ), with these transforms as weights, converges on at least some segment 
of the positive a axis including the point a = 0, then at least on that segment the resur- 
gent function / (a) is uniquely recoverable. But this set of conditions is so restrictive 
as to be practically useless. For example it has been plausibly conjectured for realistic 
QFT amplitudes |DAV82], and will be strikingly apparent again in sect. 2.4, that the 
interesting case is for all b^'s individually to have positive-s poles, but for these poles 
to cancel in the combined minor b (s) of eq. ( |2.95| ) so the total / (a) is nevertheless uni- 
valued. (This scenario seems also to be known | |CNP93U in mathematical theory.) The 
resurgent symbol is then obviously an awkward way of decomposing such an amplitude, 
but we still need to rack our brains about it because physics makes us run straight into 
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such awkward decompositions. 

In this (for the time being) unsatisfactory situation, the best one can do perhaps 
is to ask what the "escaping" parts, the parts undetermined by the resurgent symbol, 
would look like if / (a) were not to be uniquely recoverable. We already noted in passing 
that just as the escaping terms for an ordinary semiconvergent expansion must vanish 
faster than any power as a ^ 0+, the escaping parts for the more comprehensive 
resurgent symbol must be what resurgence theory calls rapidly decreasing functions, 
functions vanishing faster than any exponential exp (— const, / a). To be more precise 
we must, as earlier, adopt more specific hypotheses: let us therefore focus on the special 
case of (|]9|) relevant for QFT, and let us assume a behavior of the CLn,p coefficients 
with respect to n such that the n sum at fixed p, 

oo 

an,pX'' , where X = ua'^^e'^ , (2.99) 

n = 



is no worse than semiconvergent with factorial error bounding, analogous to eq. ( 1.13| ) 



We may then repeat verbatim the error-minimization argument with respect to N of 



sect. 1.1, and conclude that the sum of ( 2.99 ) is determined up to escaping terms with 
the small-a expansion 

e""^ X-^p {1 + O {X)}p , (2.100) 

with constants Kp (> 0) and Lp (real), and with the bracketed power series in X possibly 
being again semiconvergent. This formula is interesting in that it shows that - under 
these more special but still entirely reasonable assumptions - the super-exponential 
suppression of the "rapidly decreasing" functions turns out to really be as fast as 
exp (—exp 1/a), or in other words, faster than exp(— c/ a'^) for any c and k. For the 
semiconvergent bracket, arguing heuristically, we may repeat our chain of argument at 
the beginning of this section, and conclude that the escaping terms form the nontrivial 
part of a resurgent symbol in X, 

oo 

X-^^-" {1 + O (X)}p,g . (2.101) 

g = 

By resupplying the p summation omitted up to now and reordering the entire expansion 
by increasing K's (here of course, is another tacit assumption - that there is a small- 



est K), we find that the terms left undetermined by the symbol ( 2.92 ), if any, should 
themselves be representable as a small-a expansion of the form 



Y e"^™ ' " ^ {a^' e^)^- (a) , (2.102) 

m = 
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where Rm (a) is some resurgent symbol in a. 

In the physics context, clarity may be helped by rewriting the symbol ( p. 92 ) as 



a double expansion in powers of Aq and a as in ( p.54D , i. e. by undoing the series 
multiplication of ( 2.58| ). Call dr the mass dimension, and some typical, Euclidean 
squared momentum of the vertex function T under consideration; all other invariant 
momentum arguments of T may then be replaced by their dimensionless ratios, Xi, with 
Q^. The heuristic argument just given may then be summed up in the small-a expansion 
(we mostly omit dependences on the uninteresting set {xi} for leaner writing): 




n, p 



91 

/.2 



(2.103) 



The new escaping terms are characterized by a Gaussian decay at large momentum, with 
Aq providing the decay length. In view of their exotic coupling dependence displayed in 
eq. ( 2.102| ), they are probably outside the range of what can usefully be treated by a 
Borel-Laplace process with respect to a. 

We should now pause and reflect on the fact that here, at a point where we have 
hardly done more about resurgence theory than quote a few general concepts and basic 
formulas, we have already been touched by its great organizing power: where in the 
matter of coupling dependence we originally expected a horrible mess, we now see an 
orderly hierarchy of small-a suppressions emerging, with progressively stronger suppres- 
sion at small coupling accompanied by progressively faster decrease at large momenta. 
This insight remains valuable even while we have no complete certainty about the pres- 
ence or absence of the third line (and possibly fourth and further lines) in eq. ( 2.103| ). 
If we find ways of treating the dynamical effects leading to the second line of ( 2.103| ), 
i. e. of constructing vertex functions whose small-a expansions are resurgent symbols 
rather than perturbation series, we will have made definite progress and omitted at most 
terms with a definite, rather exotic large-momentum behavior. Moreover, we will have 
gained a framework for adequately treating an impressive array of nontrivial dynamical 
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problems. Take nonperturbative mass generation and dynamical symmetry breaking: 
the conversion of a massless perturbative into a massive nonperturbative propagator, 
described schematically by 

^ fe2 ^ cKl + ... = 

involves power, not exponential, corrections in Ag / A;^. The breaking of dilatation sym- 
metries by renormalization involves trace anomalies that are powers, not exponentials, 
ofAg. 

In my opinion, the view is even brighter, since everything in present knowledge favors 
the assumption that for the strictly renormalizable and asymptotically free theories, eq. 
( 2.103| ) is complete after its second line, just as for the superrenormalizable theories it 
was complete after the first line. (Here "complete" refers to unique recoverability from 
a small-a expansion.) Let me list a few arguments (there may be more) in favor of this 
assumption. 



(2.104) 



A purely qualitative but nevertheless strong argument is that the completeness 
postulate underlying the OPE is very probably correct. In a theory based on a 
finite number of local fields as basic degrees of freedom, it is difficult to see what 
basic operators should be available for expanding products near a point other than 
the set of all local composites at that point. Therefore the OPE, very probably, 
does not miss further escaping terms as in the third line of eq. ( 2.103| ), and 
represents a summable resurgent-symbol expansion. 

Although the question of convergence for the OPE is not settled, the n summation 
at fixed p in the OPE, eq. ( p. 991 ), in contrast to the perturbative p sum at fixed n, 
is very probably not semiconvergent but has a finite radius of convergence, so that 
in fact it leaves no room for escaping terms of type ( 2.1001) . This is because the 



growth of the coefficients a^^p with n is controlled by the growth in the number I (n) 
of local operators with mass dimension n and vacuum quantum numbers, and this 
number, in contrast to the factorial growth of the number of Feynman diagrams 
in the p direction, increases only exponentially. Consider a theory with B bosonic 
fields (pb {x) (6 = 1 ... B) of unit mass dimension in D spacetime dimensions, and 
regard a local operator product of mass dimension n as an array of n slots. Each 
slot can be filled, either with one of the or (apart from the rightmost slot) by 
one of the D spacetime derivatives d'^. (We adopt the convention of having each 
d'^ acting only on the next 4> field to the right, since the action of a on two or 
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more fields gives terms like 



{(l^a <Pb) = {d^ <Pa) h + 4>a {0^ A) (2.105) 

and therefore nothing new) . Thus the number of local composites of mass dimen- 
sion n is 

B ■ {B + DY'-^ (2.106) 

The number we called / (n) above is distinctly lower, since only the operators 
contracted to vacuum quantum numbers will give nonzero vacuum condensates, 
but in any case we have the exponential estimate 

l{n) < {B + . (2.107) 

It is true that coincidence of the n 0-field factors, where n < n, at the same 
spacetime point x leads in momentum space to n loop integrations in a vacuum 
condensate, but this increase in the number of loops is already present in the 
perturbation sum in the p direction, where it is known to contribute only to the 
power-and-exponential growth in eq. (p.5|), not to the factorial one (the p-loop 
integrals can usually be bounded from above hy a p — th power of some one-loop 
integral). Thus 



V „.n 



0'n,p — > (const.)p n ii' 



1 + O 



{n oo) (2.108) 



is the likely behavior. (Well, yes, as I hear some initiates among you saying, there 
are still the "renormalon" type sequences of p-bubble diagrams that we will also 
meet in sect. 2.4 and that generate their own p ! growth in the p direction in 
the manner exemplified by eq. (2.139) below - but remember we are dealing here 



with behavior at fixed p. In the n direction, what we have diagrammatically is 
growth in the number of lines entering and leaving a single composite vertex, at 
constant number of interaction vertices in the remainder diagram, and this does 



not generate bubble chains). Thus the series of ( 2.99 ) is likely to have the finite 
convergence circle (in X) of a generalized geometric series, such as ( |1.39D , and the 
question of escaping terms ( 2.10C| ) does not arise. 



Dimensionful observables in lattice-QFT calculations, in the removal-of-regulator 
limit, scale like a power of Aq as functions of the coupling, but not like 



(power of Ao) x <^ cq + ci e + ... \ (2.109) 
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(I am, of course, simplifying a little: lattice calculations are parametrized in terms 
of a bare coupling gQ (a) at lattice spacing a rather than a renormalized one, but 
the dependence of (Ao)iattice ■ a on go is still the same as for our A / /i in eq. ( 2.25| ) 
on g). Indeed this behavior is routinely used by lattice theorists as an indicator 
that they are sufficiently close to the continuum limit. This is a special case of the 
general statement that momentum or coupling dependences as in the third line of 
( |2.103| ) have nowhere been encountered in any of the nonperturbative approaches 
applied so far to QFT. 

There are therefore sound reasons for believing that strictly renormalizable 
and asymptotically free field theories will be about summable resurgent symbols, 
just as the superrenormalizable theories are about summable perturbation series. 
The methods described in Chapter 3 are tuned to this exceptation. 



2.4 A Resurgent Model Amplitude 



Several features of amplitudes having resurgent coupling dependence through the 
lambda scale, such as their peculiar singularity structure with respect to the coupling, 
and the tricky things that can happen when such an amplitude is approached through 
decomposition into its semiconvergent perturbation series and a nonperturbative remain- 
der, are well illustrated by a simple model function that can be analyzed completely in 
both the a and Borel planes. This model depends, other than on the coupling a, on 
a single scalar Euclidean momentum variable k'^, as does a two-point function in QFT, 
and explicitly on the sliding scale fi, but since it is taken to be dimensionless, it can 
depend on the latter two arguments only through the dimensionless ratio k'^ / It is 
given by the spectral representation 



r ( ^; a 1 = 1 



^ k^ + 



(2.110) 



with the deceptively simple spectral function. 



p =e{q^ - (p, a)) 



a 



2 Va? 



./^^ / 1 + a2 In 

where as a caricature of a QFT lambda scale we adopt 



(fi, a) = fj.^ e a <^ a (fi) 



In 



A2 



-1 



(2.111) 



(2.112) 
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The first, "nonperturbative" factor in ( 2.111 ), which simply limits the spectral inte- 



gration to > A^, schematically represents the assumption that the intermediate 
states contributing to T have through nonperturbative mass generation acquired invari- 
ant masses of at least A. The second, "perturbative" factor mimicks a typical resum- 
mation of a perturbative subseries; it displays the logarithmic dependence typical of 
perturbative amplitudes, and contains a positive constant r] characteristic of its (imag- 
ined) building-block diagram. We will always restrict 77 to 

< 7/ < 2 , (2.113) 

since at = 2 the amplitude would cease to exist, even in its renormalized form of eq. 
( 2.115| ) below, due to a divergence of the integral at the lower limit q^ = h? . (This is 



merely a consequence of the fact that in order to minimize kinematical complications, 
we have avoided building correct threshold behavior into the spectral weig ht (|1T|). If 



we were more fussy about threshold behavior, we would include in the p of eq. ( |2.111] ) 
an extra factor 

/q2 iq2 _ A2) 

that would preserve integrability at q'^ = A^ even for r] = 2, but this would also burden 
us with an extra parameter, and would deprive the model of the useful generating- 
function property mentioned after eq. ( |2.222 ) below, thus rendering the whole exercise 
much less transparent and instructive). 



Since ( |2.111 ) decreases too slowly at large q'^ to make the spectral integral con- 



vergence at the upper limit, we apply a caricature of renormalization, rewriting the 
amplitude as 

r «) = r(l; a) + a / a) (2.115) 
in terms of the derivative function 



00 / ij2 







while absorbing the ultraviolet divergence in the constant 



r(l; a) = 1 - / d{q') / . (2.117) 

2 + p- 
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It is the finite amplitude F of eq. ( p.ll6| ) , rather than the original quantity of eq. ( 2.115| ) 
which serves more or less as a pretext, that will be the main object of our study. 



We first take a look at the singularity structure in a of (2.110) or (2.116) with the 
spectral weight ( p. 111 ). From the " perturbative" factor of the latter, the amplitude 
receives a ^^-independent singularity at 



1 + a — in 

2 VA* 







while from the spectral denominator we have a A;^-dependent singularity at 

+ A;2 = , 



(2.118) 



(2.119) 



which when varying k"^ produces a branch cut in the k'^ plane at " Minkowskian" values 
k"^ = — ... — oo, the standard pattern in QFT two-point functions. Here, however, 
we are after the a plane: in both conditions ranges from to -|-oo so we may 
parameterize 

g2 = A^e* = /i2e(*"^) , t = Q...oo. (2.120) 

Then since a 

In(gVAi^) = - 1) condition (|2ll8D gives a "perturbative" coupling- 
plane singularity line at 



a 



1 



... oo , 



i. e., a branch cut on the negative real axis at 



a 



-OO...0 . 



(2.121) 



(2.122) 



If this were the only singularity with respect to the coupling, we would again be in 
the enviable situation of having Phragmen-Lindelof analyticity in a sector of infinite 
radius and opening angle 2 (vr — 5) with infinitesimal 5. The physics behind this kind of 
branch cut - which, incidentally, is also inherent in Dyson's QED argument [ pYS52 | — 
is generally that the QFT 's underlying an amplitude such as ( ^.llOj ) become unstable at 
real negative a and produce complex, two-valued Euclidean correlation functions, whose 
imaginary parts are related to the decay width of the unstable vacuum state. 



The spoiler of this simple picture is the singularity arising from ( 2.119 ), which upon 

(2.123) 



using ( 2.120 ) reads 



el' = 
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or, since k'^ is Euclidean and therefore positive, 

1 f k'^\ 

t - - -hi \ -^ \ = In(-l) = ±z7r(2n + 1) ; n = 0, 1, 2, ... (2.124) 

At a given Euclidean k'^ , we therefore have a whole infinite sequence of singularity curves 
in the a plane, 

a{t) = j-^ ^; i = 0...oo; n = 0, 1, 2, ... (2.125) 

(t - In (gj T Z7r(2n + l)j 

By taking the real and imaginary parts of this equation and eliminating the parameter 
t between the two, one finds 

{Reaf + {ImaT Vnf = rl , (2.126) 

i. e., the equations of two sequences of nested circular arcs in the a plane, with radii 

' ' ' , (2.127) 



2(2n + l)7r 2-k Gvr lOvr ' 

and tangent to the real axis from above and below at a = 0+ (see Fig. 2.3). The k"^ 
dependence is in the lengths of these arcs, which are more or less than those of semicircles 
according to whether A;^ is more or less than /x^. (The situation depicted in Fig. 2.3 
corresponds to the first, that of Fig. 1.1 (c) to the second case). 
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Im a 



k -dependent 




Re a 



k -independent 



Fig. 2. 3 Singularity lines of the model amplitude in the complex a plane 
at Euclidean k'^. In addition to the /c^-independent branch line 
on the real negative a axis, only the first two pairs from the 
infinite sequence of nested, /j^-dependent branch lines are 
shown. 



Obviously, this "nonperturbative" kind of singularity is very general: it arises solely 
from the fact that the scale of the spectral integration, or of the k'^ branch cut in the 
Minkowskian, is set by a threshold proportional to the strongly nonperturbative A^. This 



is why t 'Hooft, in discovering this peculiar family of singularities [THOTE], did not even 
need an explicit amplitude function. It is clear that this set of singularities precludes 
not only the use of theorem ( |1.24| ) to ascertain uniqueness of reconstruction from the 
perturbative expansion, but also any form of sectorial analyticity - the "horn-shaped" 
region of analyticity defined by the two largest arcs allows for no finite sectorial opening 
angle along the positive real axis. Although we have emphasized that the conditions 



of theorem (1.24) are sufficient but not neccessary conditions, this has generally been 
taken as a plausibility argument against the possibility of unique reconstruction, and the 
t 'Hooft singularities ever since have hovered as a kind of vague threat over the entire 
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subject of coupling dependence. This still should not deter us from proceeding with 
the analysis of ( 2.116 ), since after all at physical a, which are real and positive, that 
amplitude is manifestly well-defined and single- valued, and since we know by now that 
there is more to physics than the perturbative expansion. 

To examine the perturbation series of the amplitude ( |2.116 ), it is convenient to use 



In 



In 



/.2 



(2.128) 



as integration and external variables respectively. Then 



dx ■ 



1 

" 2 Q 



cosh {x - X) I + arjx 



(2.129) 



Look first at the Taylor coefficients of F with respect to a at a = 0-|-: you quickly 
realize that the nonperturbative a dependence through the lower integration limit does 
not contribute here. For example, 



dF 
da 



1 
2 

1 
2 



dx ■ 



1 



1 

" 2 a 



_d_ 

da 



cosh (x — A) 



1 

2a 



'_d_ 

da 

1 



1 



1 -|- arjx 
1 



cosh {x — X) I + ar]x 



(2.130) 



1 

' 2 a 



As a — > 0-|-, the second term vanishes exponentially, and so do all terms in the higher 
derivatives (9"F/5a" except those in which all n derivations act only on the "pertur- 
bative" factor of the integrand: the lower limit in ( 2.129| ) only creates the exponentially 
suppressed, and therefore "escaping", a dependence we described in ( 1.20| ). Thus the 
perturbation series is identical with the formal expansion around a = of the integral 



ppert 



a 



+00 



dx 



1 



cosh {x - X) I + a'qx 



(2.131) 



Accordingly, we would call the remainder integral 



/x2 



1 

2 a 



dx ■ 



cosh^ {x - X) I + arjx 



(2.132) 
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the nonperturbative part, so that 



F 



a 



jppert _|_ jpnonp 



(2.133) 



But something weird has happened here: the prescription to let a 0+ in the Taylor 
coefficients - that is, the very definition of perturbation theory - has forced us to take 



the lower integration hmit in ( 2.131 ) to minus infinity, so that now we not only integrate 
beyond the circle of convergence, |x| < of the perturbative integrand ~ a sin for 
which we expect punishment in the form of a semiconvergent perturbation series -, 
but inevitably also over its pole at x = — 1/r/a, thus creating a branch cut and the 
associated multi-valuedness in a on the entire real a axis in the F^^^^ and the positive 
a axis in the _F"°"p portion. (On the other hand it is obvious that in t'Hoofts 
"circular" a-plane singularities of Fig. 2.2 have disappeared: with no more setting 
the scale of the q'^ integration as in ( 2.12[)| ), condition ( 2.119| ) no more translates into 
a singularity condition in the a-plane.) Clearly this problem is entirely spurious, since 
in the full amplitude ( 2.129| ) at physical positive a, the integration starts only at x = 



— j^ and never touches that pole. Thus if we want to interpret ( 2.131 ) as the sum 
of perturbation theory, we will have to specify a prescription of how to approach the 
spurious branch cut, and be careful to use the same prescription in the nonperturbative 



piece (2.132) so the two nuisances can ultimately cancel. Looking back at (2.128), we 
see that x = — oo corresponds to having = 0, rather than the "true" = A^, as the 
lower limit in the spectral integral of (2.110) or ( 2.116| ): perturbation theory has created 
for us a spurious infrared problem by forcing us to integrate over massless intermediate 
states that are plainly absent in the full theory. Again it is immediately clear that this 
problem is very general, and of significance far beyond the present simple model: the 
full theory has the spontaneous A scale in almost all of its thresholds (the only exception 
being electroweak amplitudes involving the massless photon), but perturbation theory, 
by "losing" that scale, gets us involved in false infrared problems that will manifest 
themselves in false nonsummability. 

For the moment, we may ignore this problem, since we are only interested in the 



formal power-series expansion of ( |2.131 ), 

^ p=0 

whose individual Taylor coefficients are well defined: 

1 



jp{p)pert^^-^ = -{-r]Y 



dx 



+ A)f 



cosh 



X 



(2.134) 



(2.135) 
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Here we used x = x — X. The evaluation is facilitated by the observation that i^'P^^* 
possesses a scaling property of the renormalization-group type: by writing in ( 2.131 ) 

1 



1 + arj {x + X) 



1 + Q A 



a 



1+ T+^h^' 



(2.136) 



one sees that apart from the square-bracketed prefactor, F^^^^ at general k'^ looks like 
the one at k'^ = jj?, i. e. at A = 0, but with a /c^-dependent effective coupling: 



ppert 



-: a 



^(A;2) 



a 



FP^"*(1; ^(A;2)) 



a(r 



a 



(2.137) 



(2.138) 



It is therefore sufficient to evaluate the coefficients (^.135| ) at A = 0. This can be done 
with the help of the formula []GRY94 ] 



dx ■ 



cosh X 



2P-1 



1 - ^) C(p)r(p + 1) (p / 1) 



(2.139) 



and of the symmetry properties of the integrand, giving 

F(P)P-^^{X = 0) = [(-r + 1] (l - ^) C(P) (|)%! (2.140) 



Sine C ip) ~^ 1 as p oo, this indeed displays the factorial- growth pattern of (|0^), 
with constants cq = I, Bq = 0, and Aq = ^. Due to the symmetric range of the x 
integration, "alternating" and "monotonic" contributions are present in equal mixture 
so only even-p coefficients are nonzero ~ an oversimplified feature of the model. This 
feature will disappear anyway when calculating the full expansion ( p.l34| ) at k"^ ^ /x^ 
from the scaling relation ( 2.137 ): by expanding the last factor in each term of 

)oo 
= FiP)P^ri{X = {))aP[l + q(7?A)]-(^' + i) (2.141) 



j-ipert ( ^ 



p = 



in powers of q, and collecting powers of a, one finds the series (|l3|) with 



E 4'^ (2A)^ 

q = 



(2.142) 
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fip) = [1 + {-)P-1] 



1 



1 



1 



C{p - q) ■ 



(2.143) 



These coefficients, for all the simplicity of the model, again display typical properties of 
QFT perturbation series: the k'^ dependence through logarithms, which in this case are 
infrared singular at A;^ = due to the spurious infrared problem we produced, and the 
fact that the p - th term contains a p - th degree polynomial in these logarithms. The 
" leading- logarithms" (LL) subseries, comprising the terms with q = p, has fp^^ ■ p\ = 1 
(use C (0) = — ^), and is therefore, as usual, a purely geometric series. 



F: 



ipert 



(LL) 



^ p = 



1 



1 + a (r? A) 



(2.144) 



which is identical with the leading term in a (the prefactor) of relation (2.137) and, 
moreover, with the complete perturbative factor of our spectral function ( p.lll| ) taken 
at fc^ instead of - in LL approximation, that factor seems to sail unchanged through 
the folding operation of ( |2.131 ). It is clear from a look at ( 2.142 ), however, that such an 
approximation barely scratches the surface of the resummation problem: as mentioned 
already in chapter 0, the summation of subseries with finite convergence radius usually 
evades the difficulty of the p\ growth that resides here in the g = (i. e. A = 0) terms. 

For later use we note a different, still formal representation of F^^'^^ that follows 
upon using the property 



1 



1 



2P- 



Cip) 



E 

n = l 



(-) 



n-1 



nP 



(2.145) 



of the zeta function |GRY94]. In the ensuing, double formal series 



* (1; a) 



E(-r<lE[H^ + i]^'' 

n=l I p=0 



2n. 



(2.146) 



we would be tempted to identify the inner series as the asymptotic expansion of a 



combination of Ei functions according to (1.63): 

^ oo 

- EH"2n 



i^r* (1; a) 



r]a 



n = l 



-liL ^ / 2n 
e Ei 



2 n 

ei" Ei 



In 
T] a 



(2.147) 



but even this partial resummation would confront us with the necessity of supplying a 
prescription for approaching the spurious branch cut, which at real and positive values 
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of a would reside in the second Ei function of ( 2.147 ) at each n. (The function at 
k"^ / jj? 7^ 1 would then again follow from the scaling relation ( p.l37] ), i. e. by making 
the substitution 



1 

a 



1 



a 



(2.148) 



and supplying the prefactor). Exploiting a representation of type ( 2.146 ) will therefore 
be much easier for the Borel transform. 

Borel-Laplace analysis of the F?"^^* amplitude begins with the observation that the 
right-hand side of the scaling relation ( p.l37| ) is precisely of the form appearing in the 
9th line of our Table 1 (page 30) on the left-hand side, with w = — ryA. Thus 



If- 



ert 



r; s 



-{rj X) s inert 



V^'' (1; s) 



(2.149) 



is the Borel image of the scaling relation ( |2.137| ). The "Borelized" perturbation expan- 
sion aX k'^ / ^j? = 1 has coefficients ( 2.14C1| ) without the p\ factor: 



,pert 



(1; s) 



E{[(- 

p = 



+ 1] 1 



(2.150) 



Its circle of convergence is \s\ < 2 / rj. Analytic continuation beyond the circle is, in 
this case, possible by summing the series in closed form. One may, for example, apply 
relation ( p.l39| ) in reverse after having divided it by p\ : 



hP^""' (1; s) 



dx 



cosh X 



-OO 
+ 00 



OO ^ 



p = 



dx ■ 



^rj sx 



cosh X 



dx 



cosh (rj s x) 
cosh^ X 



(2.151) 



where we used the symmetry of the integrand. Use of the integral formula [ GRY94 ] 



dx 



cosh {y x) 
cosh^ X 







r (1 + I) r (1 - I) 
r(2) 



(2.152) 



and of a well-known property of the gamma function | MOS66 |, 

vr z 



r(i + z)r(i - z) 



sm Ivr z 



(2.153) 
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gives finally 



and, therefore, 



(1; s) 



sm 



pert 



; s 



sm 



(2.154) 



(2.155) 



The consequences of our having spuriously integrated over the perturbative pole in 

( 2.131| ) are now manifestly on us: this Borel transform has poles on the real s axis 

at 2 

s = ±-n , n = 1, 2, 3, ... ; (2.156) 
V 



and while those at negative s do not obstruct unique Laplace reconstruction ( |1.45| ) of 
jppert physical values of a, those at positive s do. The perturbative expansion of F is 
not Borel summable. We already pointed out that this property has at best indirectly 
to do with the presence of t 'Hooft 's set of circular-arc singularities in Fig. 2.3, since 
in i^P^*^* these are altogether absent; instead it arises from our having created, through 
the very definition of the perturbative expansion, a spurious branch cut in F^^^^ on the 
real positive a axis - that is, from our having lost the threshold that gave rise to the 
t 'Hooft arcs in the first place. I am emphasizing this in order to bring more precision to 
the vague notion that the perturbative nonsummability in, for example, QCD amplitudes 
is somehow connected to the t 'Hooft singularities. 

It is useful for later purposes to make the poles of bP^^^ more explicit, either by 



looking up the Mittag-Leffler expansion of the inverse sine in ( |2.155 ), or by using the 
representation analogous to ( 2.146| ), with a replaced by s and p\ omitted, and with its 
inner series therefore summing to 



1 + 



.2nJ 



+ 



Either route leads to the representation 



rj ^ 



n = \ 



It consists of a meromorphic part. 
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poles n = l 
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.2nJ 



-n 



+ 



+n 



s + [^]n 



- I n 
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(2.157) 
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and an entire remainder function. (Remember e^^ 
meaning to the Laplace integral 



1? I [i}). In trying to give a 













a J 



(2.160) 



we are therefore forced |ECA81| to specify a path of integration around each of the 
infinite sequence of positive-s poles in the second term of (2.159): at each of these we 
may deform the path of integration along the positive real s axis to run slightly above 
(s s + ie) or slightly below [s — > s — ie) the pole. We thus face an infinitely 
multivalued function of a, whose branches may be indexed as 



ppert 
{(Tl , CT2, 0-3, ...} 



(T„ = +1 or — 1 (all n) 



(2.161) 



and are computed by replacing the second term in (2.159) by 



E 

n = l 



/f2 



? I n 



s + iear, 



- I n 



(2.162) 



Among these, the two branches where all ct„, 's are +1 or -1 may be characterized by a 
single overall shift, s — > s + ie or s — > s — ie, of the Borel variable, which by virtue 
of ( |1.47] ) may in turn be transferred to the variable a: 



b{{t ± ie)a) 



b{ta lb ie 



b{t{a ± ie )) . 



(2.163) 



So these two branches meet at the branch cut on the positive real a axis of what by 
definition is the physical sheet with respect to a of the function i^P*^*"*: 



ppert 



jppert 



11^ 



;] a + iO 
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iO 



-ipert 
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(2.164) 
(2.165) 



With this clarification of branch-point structure, both the perturbative amplitude 
ppert main branches ( p.l64| ) and ( |2.165| ) and its Borel transform may now 

be given in well-defined closed forms involving known transcendental functions. The 
function ||M0S6q] 



Giz) = i; 



z + 1 



HI 



E ^ 

n = " + ^ 



(2.166) 
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and its property 



may be used to express 



G{z + 1) 



-G{z) + - 
z 



(2.167) 



G 



s] + G [I 



(2.168) 



where the first term in the square bracket contains the contribution from "harmless", 
negative-s poles, while the second term brings in the undesirable poles at positive s. 
Analogously, the integral representation ( 2.131| ), now given a precise meaning as a rep- 
resentation of the two main branches through the a — > a it zO prescription, may be 
written (at fc^ = /x^) 



dx 



cosh X 



1 



+ 



1 



1 + ar]x 1 — (a±iO)77x 



(2.169) 



and compared to an integral representation of the function ip' [z) = In F (z) / dz'^ 
| MOS66 | along its branch line Rez = ^, 



i I f 1 i 
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Try \2 Try 
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dx 



cosh x 
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1 + yx 1 — {y X ± iO) 



(2.170) 



to give 



1 i 

2 TTija 



1 . 



(2.171) 



These are evidently unphysical constructs, since an Euclidean amplitude at real momenta 
and physical couplings must be real and single- valued. 

Enough of the perturbative portion - there is more to life! We turn to the nonper- 
turbative remainder ( 2.132 ), at first treated formally since again we wish to derive only 
a formal resurgent-symbol expansion. In the variable 



y 
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2x + - 

a 



1 , 

h In ^ 



(2.172) 



we have 
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dy 
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{y + ^) + A) (1 
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(2.173) 
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We expand 



i-)"-^ re-(^ + 2A)l" g-„y (2.174) 



4cosh2(iy+(^ + A)) 
and obtain in the form ( 2.70| ) of a resurgent symbol, 

(1,2 \ ^ / 2 \ ^ 

^'«) = E (f) /„,/(«), (2.175) 



with support at the points s„ = n, n = 1, 2, 3, and with the formal series resulting 
from the expansion of 

oo 

/„(») = (2.170) 
as coefficients. By using t = ny as an integration variable, we obtain 



a 



fn{a) = {-r-'h [-) ; n = 1, 2, 3, ... , (2.177) 



where 



2 - r] 



oo 

e-* 



/i («) = - n Z dt ^— ^ (2.178) 



V2-»? 



a t 



has the formal expansion 



/.,/(") = -(j4^)E(^)''!"»'- 

The series fn,f therefore all exhibit the factorial-divergence pattern of eq. ( p.5| ), with 
constants 

Cn = i-TiT^ , Bo = , An = n (^—^) , (2.180) 
2 - r/ \ V J 

and, thererfore, also a pattern similar to the perturbative expansion, whose coefficients 
(at k"^ = fi'^) were given by ( 2.14C| ). For higher n's they all follow from the series fij 
by simple dilatations, a — > ^. This is a special, simple form of the phenomenon of 
resurgence described in sect. 2.3. With a look back at eq. (|1.59| ), fn may alternatively 
be expressed as 

/„(„).(_rl!ie-^i;i(P^) . (2.181) 
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It is one of the oversimplified and unrealistic features of our model that at k'^ ^ /i^, 
the fnj series for n > 1 do not develop the perturbative logarithms, 2 A = In (k'^ / /i^), 



(2.182) 



that were present in eq. ( 2.142|) for the perturbative series 

/o,/(a) := Ff^* ") ■ 

This is essentially a consequence of the simple factorizing structure of eq. ( |2.111 ). The 
only fe^ dependence in ( |2.175| ) is therefore the purely nonperturbative one of "power 
corrections" through the appearance of n - th powers of the dimensionless parameter, 



u 



-2 A 



(2.183) 



just as in the generic form of eq. ( 2.92| ). In view of the lambda-scale definition ( |2.112| ), 
this gives to the complete resurgent-symbol representation of F, 



(2.184) 



the form of an operator-product expansion as in ( 2.54 ) 

It is easy to find the "Borelized" /„ series, 

n \ °° / \ p 



K.f(s) = (-)' 



rij ^ \2 
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(n > 1) 



and their analytic continuations through resummation, 

'2\ 1 



bn{s) = (-) 



n-1 



n 



^/ s + n 



(2.185) 



(2.186) 



- n 

■n 



Each has a pole on the positive real s axis, but the interesting observation is where these 
poles are when we form the combined Borel transform ( |1.19| ) of the nonperturbative part 
(|2.175 ) by using the last line of Table 1: 



^nonp 
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J2 0(s - n) 



(I) n i-e-^r 



n = l 



- I n 



(2.187) 



The poles, because of ( 2.113| ), are at points where the step functions are unity, so in 
order to isolate them we write 



G (s - n) = 1 - e (n 



(2.188) 
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and have 

j^nonp 



fc2 
/i2' 



£ e(n - s) 



n 



n = l 



n 



1 (i) 



n 



(2.189) 



What we suspected all along is now manifest: the second term, which sums up the 
poles on the positive Borel axis, exactly cancels the corresponding second term in the 
meromorphic part (2.159) of If"^'^^. The full Borel transform 



b { s 



yiert _|_ j^nonp 



(2.190) 



is only piecewise smooth on the positive s axis due to the first term in ( 2.189] ), but 
otherwise perfectly well-behaved in the Laplace integral, which therefore yields a well- 
defined function having the resurgent symbol (2.184) as its formal small-a expansion. 

It only remains for us to check whether the Laplace reconstruction indeed gives back 
the original function ( [21291) . We use (|2l68|) and (|2l89| ) to write 
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where R collects the pieces with smoothed-out, positive-s poles: 
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(2.192) 



For the first term in (2.191), use of relation (2.167) and of the integral representation 

oo 

'z + r 



| GRY94|] 
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dx ■ 



cosh X 



(2.193) 







gives, after a few lines of calculation, the contribution 

oo 



dx 



cosh X 1 + ?? a (x + A) 



(2.194) 



to the Laplace integral. Note the lower limit of zero, rather than of minus infinity as in 
ppert^ In the second, R part of (2.191), the only minor trick you may need is to observe 
that in the term coming from the step functions, which may be written as 

1 



1 
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dy 
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n = l 
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-(f+2A) 
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with integration variables y = s / the convergent series in square brackets can be 
resummed in the manner of eq. (|2.174|) above. The rest is straightforward and the 
contribution from the R term becomes 



1 

2 



dx 



2 a 



cosh^ X I + ria{x + \) 



(2.196) 



The two contributions together indeed reproduce the function ( 2.129| ). The Borel- 
Laplace process at physical values of the coupling a has worked perfectly, unobstructed 
by either the t 'Hooft set of singularities, or by the spurious poles on the positive Borel 
axis introduced by the artificial perturbative-plus-nonperturbative decomposition. It 
has worked because we based it on the full resurgent symbol as the adequate asymptotic 
representation of a resurgent amplitude, rather than on the inadequate perturbation 
series. 

The spurious right-hand poles in the perturbative Borel transform (second term of 
eq. ( p.l59|) ) are of course, as some of you already noted, a not-so- unrealistic model of 
the so-called infrared renormalons, just as those on the left-hand Borel axis (first term 



of eq. (2.159)), which are not spurious at all but remain a true and important feature 



of the full resurgent amplitude, represent the so-called ultraviolet renormalons [GNV| 
I WEI96II ||ZNJ89|] . (The completely symmetric manner in which both arise here from eq. 
( 2.154| ) is presumably among the oversimplified features of our model.) The renormalons 
- another specialty on which we will be unable to dwell in any detail [|BEN99|| - generally 
arise through one of the perturbative partial resummations with finite 
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Fig. 2. 4 Example of a family of diagrams whose summation in perturbation 

theory gives rise to renormalons. The diagrams represent loop contri- 
butions to the two-point function of an elementary or composite opera- 
tor denoted by a cross. 



convergence radius mentioned in the introduction: in this case, the geometric-series 
summation of chains of two-point or four-point bubble graphs inserted into a larger loop 
for some two-point correlation (of elementary or composite fields), as illustrated in Fig. 
2.4. Although in such a diagram family there is only one chain diagram per loop order 
p - namely, the one with p — 2 bubbles -, their insertion into the larger loop generates 
p\ growth by the mathematics expressed in formula ( 2.13S| ), so that here one has a 



mechanism producing factorial divergence independent of the " standard" mechanism of 
growth-in-number-of-diagrams. That these infrared renormalons must get cancelled, on 
the level of the combined Borel minor ( p. 95 ), by opposite contributions from what we 



would now call the higher elementary resurgent symbols, seems to have been inferred 
first, in the context of an OPE analysis of the nonlinear sigma model, by David [DAV82|, 



who also conjectured the cancellation would persist in realistic theories. What the 
present model study can add to his (exact but somewhat indirect) argument is to show, 
with the greatest simplicity and explicitness, how the problem arises directly from the 
definition of perturbation theory, and how it disappears when the resurgent coupling 
dependence of an amplitude is viewed and treated as a whole. 

What, then, are the morals to be drawn from this model? It seems to me there 
are at least two. First, there is absolutely no reason to be afraid of t'Hooft's a-plane 
singularities: in the perturbative portion, they are missing anyway; in the full resurgent 
function, they do not impede BL reconstruction. In my opinion they are actually a sign 
of health - they indicate that the nonperturbative thresholds proportional to Aq have 
been left intact, and have not been perturbatively mutilated. 

Second, as for the recoverability of the amplitude from its resurgent symbol, it is 
clearly imperative to keep together the perturbative and strongly nonperturbative por- 
tions as much as possible, since only the sum of both leads (at positive q) to an univalent 
function. As much as possible? Since we were dealing here with a simple model ampli- 
tude having an explicit integral representation, we were able to "sum the whole thing", 
i. e. to perform Laplace reconstruction from the combined Borel minor in the form of 
eqs. ( ^T94|) and (pl96|) - but that route is clearly not viable in the real world where 



all we may know about an amplitude is that it should be (part of) a coupling-resurgent 
solution to some complicated Dyson-Schwinger hierarchy. But here a simple idea almost 
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suggests itself: if we cannot hope to sum everything in one blow, can't we at least es- 
tablish what one might call a quasi-perturbative series, i. e. an expansion in which the 
n summation in the "nonperturbative direction" has been performed, either completely 
or in some systematic approximation, at each fixed p7 In the physically relevant case of 
the real-and-equidistant support, eq. (p.92|), we would then have something like 



/ (a) ~ X] 



(2.197) 



a semiconvergent series in a in which the coefficients Cp carry the strongly nonpertur- 
bative a dependence through Aq oc exp (— si / a). Since each Cp keeps the perturbative 
(n = 0) and nonperturbative (n > 1) terms together, and since we know the full 
function is ultimately univalued at real positive a, we may be bold and conjecture that 
this series should be Borel summable as a formal series in a, i. e. if we treat the Aq 
dependence in the Cp as a dependence on an extra parameter unrelated to a. It will 
be interesting, and yield a little harvest of extra insight into the resurgent structure, to 
check this out for our model amplitude. 



We write the complete resurgent symbol Fj of eq. ( |2.184 ) in the form 



oo oo 



7 = E Y.^-^v[^) «^ (2.198) 

n = p = ^ ^ 

where the coefficients be assembled from eqs. (P.142| ) and (|2.179D : 



E 4'^ (2A)^ 

.q = 



+ (1 - 5, 



n,0) 



(l _ nP j 

Performing the n sum at fixed p then gives the form of eq. ( |2.197 ), with coefficients 



(2.199) 



Cp(A ) = ao,p + 2^ ffln.p ( ) 



(2.200) 



n = l 



Here, the sum of the n > 1 terms. 



(2.201) 
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converges in the domain 



A2 



A;2 



< 1 , 



(2.202) 



which illustrates our general expectation that the n summation should not be semicon- 
vergent but have a finite radius of convergence. Within that radius, it can be performed 
in terms of a known class of transcendental functions, the polylogarithms, 



Lip (z) 



E- 



n = l 



(kl < 1) 



(2.203) 



whose p = 2 member you may have met in calculations of two-loop (or three-point 
one-loop) Feynman graphs; the p = and p = 1 members are actually elementary 
functions: 

Lio (z) = ■ Lii (z) = -ln(l - z) . (2.204) 

I — z 

The analytic continuation beyond ( 2.202| ) is provided by the integral representation 



Lip {z) 



p\ 



ze 



(1 



z e" 



(2.205) 



which again encompasses the two elementary cases ( 2.204| ). (We have rewritten the usual 
integral definition of Lip found in tables [DEV84] through a change of variable and a 
couple of partial integrations.) Thus, 



Cp(A2 



: P\ 



q=0 



In 



+ 




(2.206) 



Now for Borel summability we do not seem to have made much progress - in particular, 
the telltale alternating signs that would indicate banishment of Borel singularities to the 
left-hand, real Borel axis are not visible. But wait! Up to now we failed to express both 
parts of eq. (2.206) in the same variable. Maybe in order to really fuse them into one 
quantity, what we should do is to reexpress the perturbative portion of Ff, 



7 (pert) 
/ 



E "o.p 0^ ^) 



a'- 



(2.207) 



in terms of the variable k"^ / h? , rather than k"^ / /x^ , by using 

In 



1 

a 



(2.208) 
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(remember eq. ( p.ll2| )) and reshuffling the perturbative a series. This requires a page 
of algebra but leads to a simple result, expressible as another RG-type scaling relation: 



F 



(pert) 



fc2 

^2' " 



2 



a 



A2' 1 



(2.209) 



(Note that the ^ dependence hidden in a = a (/x) is still present on the r. h. s.; we 
merely changed the scaling of the momentum variable k'^.) In other words we now have 
for a rearranged expansion 



(pert) 



F 



p = 



In 



fc2 

A2 



1 



2^P + 



In 



A2 



(2.210) 



(2.211) 



21 r = 

which is no more the perturbation expansion in the strict sense - the Taylor series of 
F aX a = Oh — but has the formal advantage of being expressed in the same variables 
as the nonperturbative terms. (Let me once more emphasize it is an oversimplified 
feature of our model that this reshuffling of logarithms should be necessary only in the 
perturbative part. More realistic amplitudes always have dependence on In (A;2 / ^2^ ^Iso 
in the modified Wilson coefficients re > 1, of eq. ( p. 54 ), and will need regrouping 

of logarithms in these quantities as well). It is then plausible that we should redefine 
our quasi-perturbative expansion to mean 



E 



p = 



(2.212) 



with coefficients now given by 



(2.213) 



rather than by eq. ( |2.200| ). Since the structure ( 2.143| ) of the fr^^ coefficients enforces 
{—Y = {—Y, we have 



'A2 



{r]/2Y 



[1 - (7?/2)]f + i 



(p) 




(2.214) 



This redefinition now pays off in a surprising way: what we have in the curly bracket is 



the r. h. s. of a known "duality relation" for the polylogarithms [ DEV84 ], a relation 
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connecting Lip's at large and small arguments: 



{-r+'Lip{z) = ^ /, 



r = 



In 



+ LL 



(2.215) 



Under this new definition we can therefore indeed merge the perturbative and nonper- 
turbative parts at a fixed p into a new quasi-perturbative object given, up to factors, by 
a single poly logarithm: 



A2 



-LL 



'A2 



(2.216) 



•[l-(r?/2)P + i 

From eq. ( |2.205| ), it is clear that —Lip {—z) at real, positive z is positive, so the al- 
ternating signs (— )^ we were looking for are now present. That the series ( p.212| ) with 
coefficients ( 2.216| ), and with / fc^ treated as a parameter unrelated to a, is indeed 
Borel summable with respect to a can be seen by yet another rewriting of the integral 
representation ( p.l2g| ), using as integration variable 

1 

t = 2x + - = -y , (2.217) 



a 



which is the negative of the one employed in (2.173) for the nonperturbative part alone. 
We get 



F l^; a 



1 



2(2-7?) 



dt- 



1 



cosh^ [i^] 



where 



2A + 



1 



a 



In 



fc2 
A2 



(2.218) 



(2.219) 



Due to the strong convergence provided by the cosh"^ factor, the function F (e**; a) at 
fixed K is obviously analytic with respect to a in Re a > 0. Note that the t'Hooft 
singularities do not now disturb this conclusion; they would arise from the "hidden" 
a dependence in e'^ = k'^ / A^ which for the present purpose we have eliminated by 
treating fixed parameter. Moreover, one easily evaluates 



dPF _ {-rj)Pp\ 
daP ~ 2(2 - 77)P + i 

and by estimating 

1 + 



tP 



cosh^ (t^) 



1 + 



2~ri 



a t 



p + i 



(2.220) 



7] 



a t 



-(p + i) 



< 1 



{Re a > 0) , 



(2.221) 
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< 4 6" *, one establishes a bound of type (|1.14| ), similar to 



as well as cosh 
what we had for superrenormalizable amplitudes in ( 1.27 ): 



QPF 



daP 



< 



V 



V 



(p!) 



(2.222) 



Therefore we have Nevanlinna-Sokal conditions satisfied in the half-plane Re a > 0, and 
F is Borel summable with respect to a. (As a purely mathematical aside, we have shown 
the function F of eq. (|2.218D at r/ = 1 to be a generating function for the functions 
2 (— )P"*"^p! Lip {—e'^), i. e. to generate them as coefficients when expanded with respect 
to a.) 



The quasi-perturbative expansion ( |2.212| ), in this respect, is better behaved than 
the perturbative one, but we should not lose sight of the fact that there is now no easy 
way back from the former to the latter. The beautiful duality formula ( 2.215| ) became 
applicable only after we replaced the perturbation expansion in the strict sense by the 
rearranged series ( 2.210| ), which is no more the Taylor expansion at a = 0+ of F. 
The rearrangement is infinite, in the sense that one quasi-perturbative term collects 
contributions from infinitely many perturbative terms and vice versa. As a result, in 
writing ( ^.212| ) we have "lost the perturbative limit": the formal limit of A ^ with 
a^'s remaining finite, which previously could be obtained by simply dropping the last 
term in each of the coefficients Cp of ( 2.206| ), cannot now be recovered any more from one 
or a finite number of quasi-perturbative terms. In the zeroth quasi-perturbative order, 



the analog of the OPE object ( 2.56 ) which sums up the terms without powers of a but 
with all powers of A^, we have 



CO [J, 



1 



(2.223) 



1 - (f) fc2 + A2 ' 

and although the limit A ^ here exists, it differs from the perturbative _p(o)pe''t = i 
by the factor (l — 2) ' '^hich sums a geometric series coming from all perturbative 
terms. In the higher quasi-perturbative orders, the A^ — > limit does not even exist 
termwise; we could retrieve it only when undoing the rearrangement of eq. ( p.21C| ) 
by applying relation ( 2.208| ) backwards - i.e. by acknowledging that A^ and a are, 
in fact, not unrelated. In a calculation where we can determine only a finite number 
of quasi-perturbative terms, we therefore have no easy way of making contact with the 
perturbative series, and of checking whether it comes out as a limiting case. One might 
take the view, why bother at all about a perturbative limit in this sense? After all, the 
very point of an RG-invariant scale is that it cannot be switched off in any way, and since 
the interesting regime we would like to explore for nonperturbative phenomena, in both 
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Quantum Chromodynamics and electroweak theory, is the regime we could 

dismiss the A — > hmit as uninteresting. Still, the difficulty of staying in touch with the 
perturbative expansion remains a nuisance in some respects. Take, for example, the RG 
results on large-momentum behavior: they stipulate that the perturbative zeroth-order 
vertex functions should be approached, with logarithmic corrections otainable through 
perturbative leading- log summations, and since these are exact limiting results one would 
like to preserve, one feels unhappy about losing the easy way back to their verification. 

The loss is offset, to be sure, by considerable gain in the description of the small- 



k regime. By eq. (2.203), we now have the correct behavior F — > as A; ^ in 



each quasi-perturbative term. What is more, due to the analytic properties of the Lip 
functions - for p > 1, they all have branch cuts at real argument larger than unity each 
term now also reflects the correct "spectral" singularities of F at Minkowskian values 
/c^ < — A'^, including the pole at k"^ = — carried by the p = term ( p.223| ) and the 
physical cut at k"^ < — A^ in the p > 1 quasi-perturbative corrections. Thanks to the 
tight linkage between the A^ and k'^ dependences, the quasi-perturbative resummation 
of (A^)" powers in the resurgent symbol has automatically furnished a very satisfactory 
analytic continuation in the k'^ plane, which resummations of perturbation theory alone, 
even disregarding its false infrared singularities, could not provide. It would then seem 
that the quasi-perturbative series deserves a try in the "real" problem, and chap. 3 will 
sketch some steps in that direction. 
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3 Construction of Resurgent Amplitudes 



3.1 The Approximating Sequence 

We have seen that the operator-product expansion, in a strictly renormaHzable and 
asymptotically free QFT, is of the form of a resurgent symbol in the coupling constant 
a. As such, it asserts that amplitudes in those QFT 's are resurgent functions of a. We 
have also seen that there are good reasons for assuming that this symbol is summable, 
that is, the small-a expansion of a unique analytic function of a. In physics terms, there 
are good reasons for expecting the OPE to be complete, i. e. to leave no room for further 
terms that would escape the resurgent-symbol representation by having stronger-than- 
exponential suppression as a — > 0+. 

The implications of these two insights are quite strong. If we find methods of con- 
structing systematically, or at least in a systematic approximation, the unique analytic 
coupling functions waiting behind the resurgent symbols, we will have accounted not just 
for one special type (as an easy misunderstanding might suggest) but for the totality of 
the strongly nonanalytic coupling dependence. It therefore seems worthwhile to look for 
and try out such methods. The method described in this, our last chapter, builds on 
the plausible hypothesis of summability of the resurgent symbols, and attempts to sys- 
tematically reconstruct the strongly nonperturbative coupling dependence of correlation 
functions - the dependence manifesting itself in the resurgent symbol through powers of 
the RG-invariant scale A - from their equations of motion. 

In developing such a method we will find ourselves compelled to deviate somewhat 
from standard strategy. Standard procedure relies on the function being a solution to 
some differential or differential-cum-functional equation. Its first step is then to obtain 
the function's resurgent-symbol expansion by purely formal manipulation, inserting the 
symbol into the governing equation and comparing terms. (If the equation is nonlin- 
ear, it becomes important at this point that the set of elementary resurgent symbols 
form an algebra under multiplication.) In mathematical theory, completion of this step 
is usually taken for granted; it belongs, so to speak, to a realm of pre-mathematical 
heuristics. The procedure then continues in the manner we illustrated in sect. 2.4 with 
our model amplitude: one "borelizes" the individual power series, separately for each 
n, and performs analytic continuation to obtain the Borel minors 6„ (s). If these are 
sufficiently nonsingular (integrable singularities, as pointed out by Sokal | SOK80(| , can 



be tolerated) and restricted in growth on the positive Borel axis, one may proceed to 
evaluating their Laplace transforms that appear in eq. ( 2.80| ). Finally, one confronts the 
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n series of that equation: if it has nonzero convergence radius, it may require another 
analytic-continuation step; if it is summably divergent, it will call for another BL resum- 
mation. If all goes well, one finally reaches the holy grail, a globally valid representation 
of the full resurgent function. 

For physicists dealing with QFT, while this standard procedure should be in the 
back of their minds, there are cogent reasons for proceeding in different order. As we 
noted in studying our model, it is a consequence of the very definition of perturbation 
theory that individual bn's are in general afflicted with spurious Borel nonsummability, 
and thus force us to cope with spurious multivaluedness at all intermediate stages of 
the reconstruction, while we are convinced that the full function must ultimately be 
univalued at physical, positive couplings. We saw that the easiest way to avoid this is to 
keep perturbative and nonperturbative contributions tightly together by using a quasi- 
perturbative expansion, that is, by first performing the n summation in the symbol's 
nonperturbative direction at fixed p. You might choose to call this a "resummation of 
the OPE in the n direction" , but you know by now why I am reluctant to use that term - 
it sounds as if the OPE, for realistic amplitudes, were something fully calculable ab initio 
that one could simply take and subject to some process of continuation. In truth, as we 
have emphasized, the OPE is not known ab initio, since is does not by itself determine its 
own strongly nonperturbative ingredients, the vacuum condensates. We will therefore 
be forced to recover our strongly nonperturbative contributions from exact equations 
of motion, which supply the dynamical principle the OPE is lacking. In the process, 
the vacuum condensates will lose some of their central role in describing the strongly 
nonperturbative effects, that role being taken over by dimensionless nonperturbative 
parameters more easily accessible through the equations of motion (and from which one 
may later recover the condensates if one so desires). 

Of course, with realistic amplitudes there is again no way of performing the entire n 
sum at fixed p in one blow, as we did for our model in writing eqs. ( 2.206| ) or ( 2.214 ). But 



we may attempt to do it in terms of a suitable systematic sequence of approximations, 
amounting to resummations of the (A^)" power series that grow progressively more 
comprehensive with increasing approximation order. These resummations, as again we 
learned in our model, will automatically provide a sequence of analytic continuations, 
increasingly accurate, to the momentum region of ^ A^. 

As we stressed already at the end of sect. 1.2, it is not compulsory for us to use 
" Borelization" in doing this - all the more as we trust the n summation at fixed p to 
already have nonzero radius of convergence. Indeed we shall find other approximat- 
ing sequences better adapted to the task. Enumerating successive approximants (we'll 
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discuss their detailed form in a moment) by an index r = 0, 1, 2, we denote 

p[np) _ approximant to c„.p ( — | (3-1) 

n = o ^'^ ^ 

at a fixed quasi-perturbative order p. 

There is, to be sure, a subtle drawback to this necessity of employing approximants 
to the n sum. At finite level r, we can no more expect a quasi-perturbative expansion 
based on such approximants to achieve the complete cancellation of spurious poles in the 
combined Borel minor that we were entitled to expect from the full n summation. This 
is not much of a problem as long as we calculate amplitudes at an energy scale /i where 
the coupling a (//) is small enough to ensure good semi-convergence for the sequence of 
partial p sums, 

rhp] ^ ^ ^[r,,) (A2)a(/i)^ (r, p = 0, 1, 2, ...) , (3.2) 

q = 

as we do in most of present-day electroweak theory: then these sums may be taken 
directly as approximants to the amplitude, without ever confronting the Borel plane. 
But when a (/i) is not small, as in the hadronic region of QCD, we will have to rely 
on a BL process for the p summation, and at low approximation levels r we will then 
be able to avoid only a small number, increasing with r, of the "nearest", i.e. next- 
to-the-origin, spurious Borel poles. Since nearest poles in the s plane dominate large-a 
behavior in the a-plane, and since in an asymptotically free theory the large-a regime 
is also the low-momentum regime, we may qualitatively expect our approximation to 
have its largest errors at low momenta and, on increasing r, to gradually work its way 
down from the large-momentum regime. Nevertheless, the crude analytic continuations 
it provides at low r to the low-momentum region will have finite error and be usable 
even down to = 0, where the error of the original resurgent symbol, an expansion in 
terms of / Q^, would become infinite. 

If we then start our reconstruction scheme by inserting such a continuation-through- 
resummation, rather than the original OPE-type expansion, into the equations of motion, 
we reap several benefits. First, it turns out it is really only for p = that we need to 
specify an explicit sequence of approximating functions: the nonperturbatively modified, 
zeroth quasi-perturbative orders T^'^'^\ whose exact counterparts we met in the form of 
eq. ( p. 56 ) for the OPE, and in the form of eq. ( 2.223| ) for our model amplitude. The 
higher quasi-perturbative corrections, T^'^'P^ with p > 1, will then be generated for us 
by the equations of motion themselves. Now choosing approximants for r'^'") is a lot 
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easier than for general p, since they do not need to be saddled with typical complications 
that only the higher-p corrections are supposed to bring - in particular, with logarith- 
mic corrections to large-momentum behavior, and with the associated dependence on /x 
(compare eq. ( p.53 ) in our discussion of the OPE). Drawing on the body of experience 
we alluded to in connection with formulas (1.28) and (1.53), we therefore again choose 
rational approximants, but this time with respect to A and not to a - i. e. expressions 
of the form 

'" = ° (3.3) 



E 5n({fc})(A2)" 
rt = 



(for bosonic amplitudes). We seem to need two indices, the polynomial degrees l' and 
I, to label such an approximant, but here we reap our next benefit, since asymptotic 
freedom greatly simplifies the matter. Recall that the RG leads to exact asymptotic for- 
mulas for vertex functions when the entire set {k} of momentum arguments are scaled up 
uniformly, and that these formulas simply feature the zeroth-order perturbative expres- 
sion, r(o)peri^ multiplied by certain logarithmic corrections. But the logarithmic factors 
arise from (resummation of) perturbative corrections with p > 1, and in our context 
are therefore for the higher T^^'P^ 's to bring; they have no business with r[^''*'\ So of the 
latter we must demand naive asymptotic freedom: 

r[;'°)({fc},A) T^^^'''\{k}) (allfc»A). (3.4) 

In exploiting this condition, it is best to proceed in steps and deal first with the simplest 
case of a two-point vertex function (A^ = 2), or negative-inverse propagator, depending 
on a single scalar Euclidean k'^; vertices with > 3 that require more momentum 
variables will be taken up later. Since A^ is dimensionful, and we assume no Lagrangian 
masses present, the approximant (3^) is then automatically rational with respect to fc^ 
too, with 

fm = Cmik^y-"" , gn = Vniky-'', (3.5) 

and with real constants Cm and ??„ since the complete approximant should be real at 
real Euclidean k"^. (At this point, things would seem to become substantially more 
complicated when Lagrangian masses rh are present, since the Cm and 77^ could then 
still depend on the dimensionless ratios / k'^ - but let us continue to concentrate on 
the massless case. ) Now the T^o)pert condition (0), which are none other than the 
ordinary Feynman rules, are polynomial in k, and in particular for N = 2 and a massless 
scalar field, 

-rf ^^'■*({A;}) = A;2 . (3.6) 



91 



Thus condition ( ^.4| ) fixes the relative degrees of our approximant in advance at I = 
I + 1. We do not face the prospect, as one usually does in working with rational 
approximants, of having to work through the entire Pade table at given total degree 
l' + /; we need to explore only the next-to-diagonal line — / = 1 of it. So in fact we 
need only one labeling index, and we will from now on adopt the convention of using 
the denominator degree I as the classifying index r for our approximants to irreducible 
vertex functions: 

r = / = /'- 1 iN = 2) (3.7) 
(For propagators, or connected two-point functions, 

D = G^™"") = -(Fa)-^ , (3.8) 

and the index r is then also the numerator degree of rational propagator approximants.) 
Since without any loss of generality we may choose 

% = 1 , go (k) = {k^Y , (3.9) 

the level-r approximant is now characterized by a set of 2r + 2 dimensionless, real 
parameters 

ct\cP,...^U; vt\vP,...vl^-^ (3.10) 

and assumes the explicit form 

_rr\k') = ^P^ (3,11) 

A2) 

for a scalar two-point vertex, where 

ivM (k^ A^) = V cM (fc2)^- + i-™(A2)- ^ (3.12) 



m = 



am ik', A') = y: vl:^ ik^y-'' i^'r = i) • (3.13) 

n = 

(Of course, if r2 has Lorentz-tensor or Dirac-matrix structure, as do the fields in a gauge 
theory, such a formula applies separately to each of the Lorentz-scalar "invariant func- 
tions" in a decomposition of the vertex into basis tensors or basis matrices. For example, 
an SU (N) gauge-field propagator Z?^^ (k), with Euclidean vector indices fi, u = 1 ...4 
and color indices a, b = 1 ... N'^ — 1, would decompose into transverse and longitudinal 
pieces with respect to k, 

Kbik) = 5atr''{k)DT{k'') + l>'''{k)DL{k^)] , (3.14) 



92 



with projectors t^'^ and l'^'^ given by 



6^" 



l^^'ik) 



(3.15) 



and the functions [Dt (^^)] ^, [Dl (^^)] ^ would each have an approximant of the form 
(|3llD). Theubiquit ous indices [r] are cumbersome but necessary at this stage to remind 
us that the coefficients of a rational approximant must be recalculated at every new order 
r. 



You may ask at this point why we did not exploit condition (| 
up to now, we only ensured 



1)1 



cfi k^ 



{k^ » A^) 



more fully 



(3.16) 



so why don't we impose Co = 1) i- e. /o {k) = [k'^Y'^^'? Here I have to remind you of 
the experience with our model amplitude, and ask you to cast a glance at eq. (p.223| ) 
above. If we insist on establishing a quasi-perturbative expansion in the sense originally 
envisaged, and exemplified by eqs. ( 2.197| ) and ( p.206| ), then indeed we should require 
(^0 = 1. (For the inverse of the function Di (fc^) defined by eq. ( p. 14 ), the correct analog 
would be 

[Diik^T^ Ik" (A:2 » a2) , (3.17) 



with ^ the gauge-fixing parameter, in the linear covariant gauge fixing usually employed.) 
But if we allow for the possibility that our approximation might turn out to be quasi- 
perturbative in the redefined sense of eq. (1121^ ) and ( ^^21^ ), we must allow for a 
constant factor Qq ^ 1, analogous to the (l — 2) ^ °^ (E!1H)' '^^ich amounts to 
a finite renormalization. As long as we have not begun to determine our approximant 
( |3.11 ) from a dynamical equation, we cannot know whether that equation will reveal a 
tendency to prefer one or the other of the two expansions, or something in between. So 
caution would advise us to keep the option of Co 7^ 1 (or 7^ 1 /^) at this time. 



It is trivial to recover from eq. ( 3.11 ) the expansion corresponding to what we met 
in the OPE context as eq. ( p. 56 ): by expanding ( ^.11 ) in powers of A^ / k"^ for k"^ ^ A^, 
we have 



A2 

A;2 



+ 



V2 



(3.18) 
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SO if we imposed Q = 1 we could identify, at level r, 
C/J^) (0, 0) = Cf 1 - r?f 1 , 



(0' 0) = Q - ^1 (^Cl - ??1 j - % 



(3.19) 



etc. Indeed it would be nothing but standard Pade procedure to determine the Cm, Vn 
coefficients from the first 2 r + 1 of these equations if the Un (0, 0) were known quantities. 
They aren't, and so eqs. ( 3.19| ) are useful mainly as a kind of sum-rule check on the 
vacuum-condensate combinations (0, 0) (Aq)" of eq. ( 2.571) that can in principle be 
computed once the C,'s and -q's have been determined from dynamics. 

More interesting, particularly with a view toward generalizations to higher vertex 
functions, is another simple rewriting of ( |3.1lD , namely the partial- fraction decomposi- 
tion of its inverse (the propagator approximant) with respect to the variable A:^. Decom- 
pose the numerator polynomial ( 3. 121 ) into its root factors. 



r + l 



(fc2,A2) = Con [k^ + '^V^') (3-20) 



then 



with residues, 



DK0)(fc2) ^ (3.21) 



= 
r + l 



Co n [o^-o? 



k = l 

(k^l) 



This may be viewed as arising from the Euclidean Kallen-Lehmann representation, 



oo 







by a discrete approximation of the spectral function, 

P2 is) pt''^ is) = X: S (s - aM) . (3.24) 



1 = 1 
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Since the spectral function is positive and has its support in the Minkowskian, we ex- 
pect a sensible approximant to have all its pi 's and o"/ 's non-negative; for the residues, 
comparison of the fe^ 3> limits of ( p.2lD and ( 3.11 ) gives the sum rule 



r + 1 



Co E p\ 



1 = 1 



The 2r + 2 dimensionless quantities 



... a, 



r + 1 ' 



[r] [r] 
Pi , P2 



■Pr + 1 



(3.25) 



(3.26) 



represent a parametrization equivalent to ( p. 10 ); when (q = 1, equation ( |3.25| ) turns 
into a relation among the pi 's that reduces the number of independent parameters 
by one. For later convenience, we also note the corresponding decomposition of the 
vertex function ( 3.11| ) itself, which proceeds from the root-factor decomposition of the 
denominator polynomial ( |3.13 ), 



(3.27) 



Since in this case the numerator polynomial is of higher degree, the decomposition now 
takes a slightly different form. 



\r,0) 



so 



r] 



2s- 



(3.28) 



with the sum now representing a discrete approximation to the branch cut of the self- 
energy that mirrors the Kallen-Lehmann propagator cut. By comparison with eq. ( 3.1^ ), 
its parameters obey 



Ci 



[r] 



(3.29) 



'■2S + 1 



Co 



Cl , 2 

Co 



etc. 



If the field whose two-point correlation D represents had exactly one particle as- 
sociated with it (as is usually assumed without discussion), then y^oi A and pi would 
approximate its mass and field-renormalization constant, and the higher pi's and ai's 
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would provide a discrete approximation to its dressing cut. But in a nonperturbative 
context, where the zeroth-order perturbative spectrum may suffer not only small shifts 
but also qualitative changes, it is quite possible for a quantized field to give rise to more 
than one particle, or none at all, and in such cases we encounter a practical disadvan- 
tage of our rational approximation: at low degrees r, there is no simple way of deciding 



whether a given pole term in (3.21) approximates a genuine particle pole, or is merely 
part of the approximation to the dressing cut. The answer can only be found indirectly 
by observing the approximation as r increases: all experience with rational approximants 
indicates that poles approximating "true" poles then display reasonable convergence of 
their positions and residues, whereas the poles approximating a cut keep "hopping", 
and at the same time spread to ever higher values of s, since at each successive r the 
approximation invents a new set of positions and strengths to achieve an improved scan 
of the continuous part of the spectral function. 

We should again pause, since some of you may have taken offense at this point 
with our apparent lightheartedness in talking about such treacherous objects as par- 
ticle masses and field renormalizations. Haven't we all learned that these things are 
ultraviolet divergent, and can at best be traded for something experimental but not 
really computed? Here you should recall that at this stage we are dealing with a zeroth 
quasi-perturbative order, and at zeroth order in powers of a, renormalization constants 
are not yet divergent. The ultraviolet divergences will catch up with us soon enough 
once we start computing quasi-perturbative corrections, but at the zeroth order we are 
only into reconstruction of the quantity that in the OPE context appeared as eq. ( 2.5£| ) 



a quantity which we noted is RG-invariant and finite. At this level, things like the 



nonperturbative mass shift of eq. ( 2.104 ) do indeed have finite values. Even when it 



comes to computing quasi-perturbative corrections, you should still remember the pat- 



tern we noted in eq. (2.49) - it is only for the perturbative portion that we expect to 
have to cancel net ultraviolet divergences. For the nonperturbative remainder, what we 
are aiming at is a direct continuation-through-resummation of the n sum ofeq. (|2lg| ) 



with its ultraviolet-finite ■ coefficients - we do not attempt to separately deal with 
vacuum condensates and Wilson coefficients, whose ultraviolet divergences then cancel 
by eq. ( p. 46 ). That the dynamical equations should provide a mechanism allowing for 



all these features is of course not trivial, but nevertheless true, as the next section will 
show. 



Let us return to eqs. ( |3.21| ) and ( 3.281) where we exploited the tight linkage, inherent 



in the dimensionful character of A^, between the A^ and k"^ dependences to reinterpret 
our original A^ approximation as an approximation to the expected /c^-plane analyticity 
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structure of our propagator or self-energy. It is this reinterpretation that provides the 
best guide for our next and shghtly more difficuh step: choosing zeroth-order approxi- 
mants for higher-than-two-point vertex functions, consistent with the two-point form of 
eq. ( 3.11| ). Here again we start with the simplest and most important case of a single 



scalar three-point vertex function with Bose symmetry, depending on the three squared 
Euclidean momenta 

kf, k\, k^ (3.30) 
of its three external legs, and being otherwise dimensionless with unit perturbative limit, 

Tf^P--\{k}) = I . (3.31) 

(We have tacitly adopted the widespread convention of defining the three-point vertex 
without a coupling factor ^q, so what actually enters into diagrams is g^V^ rather than 
Fs alone.) This again is not precisely what occurs in a realistic gauge theory - a vertex 
for three SU {N) gauge- vector-fields, for example, has mass dimension plus one rather 
than zero, and that dimension is carried by the Lorentz basis tensors in a decomposition 



of the full tensorial vertex. In the index notation adopted for eq. ( 3.14 ) above, the 
beginnings of such a decomposition would read like 



+ fc2 + fca = 



= /a6c{(5^"(fc2 - hYF^J^ (kl kl; kl) 

+ 5-f^{k3 - h)^Fy^ (kl kl; kl) 

+ 5^'^ikl - k2rF^^^ {kl kl kl) + ...} (3.32) 
+ dai,c {5^" {k2 - h) F^') {kl kl kl) 

+ 5-^(^3 _ fcO^F^'^) {klklkl 

+ 6^^^{ki - k2)''Fl,''\klklki) + ...] , 

with fabc and dabc denoting respectively the antisymmetric and symmetric SU {N) 
structure constants, which enter here in their capacity as Clebsch-Gordans for the two 
possible ways of coupling three adjoint-representation objects to a (global) gauge singlet. 
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What is dimensionless is then the Fq'^^ and Plf"^ invariant functions, but they wih of 
course display only the reduced Bose symmetry required to match the partial symmetry 
of their associated tensor structures, and thus be symmetric or antisymmetric only with 
respect to the first two of their arguments. The Feynman rule F^y^^^^ , as you will 
remember, has only the fabc terms of eq. ( |3.32| ) with Fq'^^ 's replaced by unity, so we 
have the perturbative limits 

(fr)'°'"" = l, (Fr)'"""'* = . (3.33) 

as A ^ 0. 



The dots in eq. ( 3.32| ) in each case stand for a finite but lengthy list of other indepen- 



dent Lorentz-tensor structures (actually, no less than eleven of them) , forming triplets or 
singlets under the group of permutations of the three vector-field legs, and each accom- 
panied by its own invariant function. It will emerge only much further down the road 
(and be commented upon in our closing section) that in the zeroth quasi-perturbative 
order we are considering here, this cumbersome (2 x 14)-term tensorial structure fortu- 
nately simplifies - most of the Lorentz tensors not made explicit in eq. ( 3.32| ) are tensors 



involving three powers of momentum rather than one, and being therefore accompanied 
by invariant functions of mass dimension minus two rather than zero, and in the zeroth 
quasi-perturbative order these turn out to be excluded by considerations of renormal- 
izability. Here, however, we should not let ourselves be sidetracked too much by these 
purely kinematical and enumerative aspects of tensor structure, so you will allow me 
to refer you rather summarily to the standard reference [|BCU81 ] on these aspects. For 



the issues we wish to address, consideration of a purely schematic, scalar three-variable 
function characterized by eq. ( ^.31 ), or equivalently restriction to the F^'^^ function of 



eq. ( |3.32| ) with limit as in the first of eqs. (|3.33| ), will be quite adequate. 

For such a function, a rational approximant with respect to A^, of the form of eq. 
( |3.3| ), must have l' = / by virtue of condition ( |3.4| ), so again we need only one labelling 
index. But otherwise, eq. ( |3.3D alone is now much too unspecific to represent a computa- 
tionally useful approximant. (All we can say at this stage about its fn and gn coefficients 
is that each must be of the form (Q^)'""", where is some dimension-two combination 



of the variables ( |3.30| ), times some dimensionless function of two independent dimen- 
sionless combinations.) It might be sufficient if the dynamical equations to be invoked 
were differential or integral equations with respect to the coupling, which would corre- 



spond to the situation called "equational resurgence" in mathematical theory |ECA81| 



| STS96|] . But here the equations at our disposal will be integral equations with respect 
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to momentum variables, and for an approximant to be capable of exploiting these equa- 
tions, it must then detail the momentum dependence in a manner comparable to eq. 



( 3.11 ) for the one- variable case. In addition to the resurgent coupling structure, which 
we are resumming at the p = level through functions ( p.3D , we need some additional 
input here. What can help us across this gap are considerations of analyticity. From 
the general properties of QFT amplitudes, we expect to be analytic (in the usual 
loose sense, i. e. except for a denumerable set of poles and cuts) in the complex plane 
of each of its variables ( |3.3[1| ) at fixed values of the other two, in a manner analogous 
to the single-fc^ analyticity embodied in the spectral representation ( p.23D . It therefore 
seems sensible, and not unduly restrictive, to postulate that our T^, at least at the zeroth 
quasi-perturbative order, fits into the form of a triple Cauchy-integral representation, 

Ci C2 C3 

For a Bose-symmetric function, the spectral weight p^, which may include distribution- 
valued terms, should be symmetric in its three variables, and the contours Ci in the 
three Si planes should be identical to one C. Also the full function must be real at real 
Euclidean momenta, so if C comprises any sections away from the real s axis (which for 
off-shell variables like our kf cannot be excluded), these, as well as the associated portions 
of p^^\ should only appear in complex-conjugate pairs. We will always understand the 
form ( |3.34| ) to allow for terms depending on less than three variables, by interpreting 
the compact notation p^^^ as 

pf=io n [^m^^)] + A' Y.^k^^{sm^^5{s,)]p,,,{s,) 

i=l i<j 
3 

+ J2 [klHsk)]p3,2{Si, Sj) + A6p3,3(si, S2, S3) , (3.35) 
k = l 

with a positive constant and with (z, j, k) cyclic in the second and third terms. This 
function then has large-momentum limit 

r^^"°^ ^ ^0 (all kf > A^) , (3.36) 

and all our earlier comments on the C,q^ of eq. ( 3.16| ) are again applicable. 



The main consequences of postulate ( |3.34 ), for our purposes, are that if we now 



generate an approximant to by adopting a discrete approximation analogous to 
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dm, 

(0)^ \ [r,0)/ ^ A M 

(Sl, S2, S3) > (^1' ^2, S3) = E ?rii 



n2 ns 



X [<5„,o kl + {l- 6n,o) A'] + (1 - <5„,o) A'] K30 A^I + (1 - -^nao) A2] (3-37) 

X 5 (^si - ri'^i'^ 5 (^S2 - ri^J) (s3 - ^3) , 

then automatically that approximant will (i) be rational not only with respect to 
but also to the kf, and (ii) possess a denominator factorizing in the three kf which 
mirrors the product structure of the Cauchy-integral denominator. (To accommodate 
the structure (|3.35| ) we add the convention that 

rfl = (all r) (3.38) 

in eq. ( p. 37] ).) For full Bose symmetry we should have 

^n}n2n3 Completely symmetric in (ni, ?i2, 713) . (3.39) 

Therefore if we put the entire approximant over a common denominator we arrive at the 
form 

AM(fcf;A2)Al"l(t|;A2)AW(*:iiA2) 

with a numerator polynomial of degree 3 r in and of degree r in each of the squared 
momenta, 

3 r r 

+ 7712 + n^a , 3 r ^'mi rn,2 ma 
A = mi , 7712 , 7713 = 

X {klr' (A;i)™2 (A;!)™^ (A^)^ , (3.41) 
^mim2m3 Completely Symmetric in (mi, m2, ms) , (3.42) 

■^irr ~ ?0 ^ (3.43) 

and with a denominator factorizing into three degree-r polynomials similar to ( |3.13| ), 

aM {k^; A2) = f[{e + tM a2) . (3.44) 



iVfl ({fc2}; A^) = ^ ^ 5A + mi+m2+ma,3rxM 
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In the notation of our original eq. ( |3.3D for the rational approximants, eq. ( ^.401) therefore 
has 

l' = I = 3r (iV = 3) , (3.45) 



rather than / = r as in eq. ( p.ll| ). Indeed it will turn out that on this point analyticity 
has given us good guidance, since it will be just this combination - of ? = r for r2 with 
I = 3r for Fs, and more matching / 's for the r4 vertex ~ that will prove itself capable 
of self-consistency in the equations of motion. 



When all kf > A^, eq. (l3lo| ) gives an expansion 



1 + 



X, 



r, r, 1 



?0 



r] 



E 



A2 A2 A2 

1.2 ^ h.2 ^ h.2 
'^1 1^2 '^'i 



(3.46) 



that may be regarded as a three-point counterpart of eq. ( 2.56 ) - a three-variable OPE, 
restricted here to terms of zeroth perturbation order with respect to a. 

The foregoing discussion of and Fg approximants has, I hope, provided 

sufficient illustration of the guidelines and technicalities that go into a construction of 
zeroth quasi-perturbative orders as a systematic sequence. You will therefore allow me 
to skip the corresponding construction for a four-boson vertex - the highest function for 
which such a construction will be found necessary. With no less than six independent 
kinematic invariants, and with a plethora of possible tensor structures, this function 
requires a much lengthier but not much more illuminating construction. Only much 
further along the way - that is, after having learned some lessons from trying to get 
these constructions self-consistent in the equations of motion - will it emerge that for the 
zeroth quasi-perturbative order., F^' , tensorial structures and functional dependences 
are strongly restricted, so that the final approximants are still manageable. For the 
moment, it may be more important to address your feelings that my remark about 
this being the highest construct necessary sounds completely obscure. The equations 
for correlation functions, in whatever variant, are always hierarchically coupled - so in 
order to establish a four-point amplitude, aren't we going to need a five-point one, and 
so on forever? In other words, mustn't we live in fear that unlike their perturbative 
zeroth-order counterparts - the Feynman rules y^^^^"^ - our new "rules" will endlessly 
proliferate'^ It seems nontrivial (at least to me) that such fears in fact should turn out 
to be unfounded. 
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3.2 Self-Consistency of the Generalized Feynman Rules 



The best of the benefits from our n-resummed, quasi-perturbative reconstruction 
strategy are yet to be reaped. They emerge once we insert the matched set of zeroth- 
order approximants, as constructed in the last section, into the interaction (loop) terms 
of exact equations of motion. As announced, and as expected, these loop terms then start 
producing the quasi-radiative corrections T^^'P\ p > 1, in what is basically the first step 
in an iteration process around a nonperturbatively improved starting point. But more 
importantly, we find four pleasant surprises. First, despite the strongly nonperturbative 
character of the T^'''^^ 's with their ubiquitous A-scale terms, the loops turn out to be 
readily calculable by the known methods of loop computation and regularization, simply 
due to their rational structure in momenta. Second, in addition to producing p > 1 
corrections, the loop terms prove themselves capable of reproducing the nonperturbative 
parts of our zeroth-order input, thus conferring true self- consistency on the improved 
starting point. Third, while the p > 1 corrections continue to carry the ultraviolet loop 
divergences, the mechanism operative in the self-reproduction sees to it that the non- 
perturbative zeroth-order terms, the quantities r''"''^^ - r(o)p^'"* at all levels r, establish 
themselves in a finite manner, as one would expect and demand from our consideration 
of the OPE. Fourth, the same mechanism relieves us from our nightmare of proliferation 
of Feynman rules, by strictly limiting the formation of zeroth-order proper vertices to 
a small finite set of these functions. It is these features together that seem to justify 
designation of this finite and self-consistent set of T^"^'^^ 's as generalized Feynman rules. 

The system of exact equations of motion we are going to exploit is the oldest known 
such system, the integral equations first established by Dyson [pYS49 |, who derived them 



by Feynman-graph summation, and by Schwinger SCH51 ], whose functional derivation 
made it plain that they are in fact independent of any perturbative framework. They are 
referred to as the DS equations. Although they still were fundamental field-theory tools 
for such researchers as Schwinger, Bogolyubov, and Symanzik, they are rarely discussed 
in textbooks nowadays (two 1980 's texts that give them more than passing mention 



are [tTZ80| and [ |RIV87 |), but I hope you may have had at least a superficial encounter 



with them and will recognize a DS system (from its typical, hierarchical structure) when 
you see one. For the DS equations of non-abelian gauge theories, the standard original 
reference is [EIF74]; condensed statements of their derivation and main properties can 



also be found in the review articles |ROW94] and AVS01 |, which concentrate on exact 



numerical solutions - another line of research which our focus on coupling dependence 
will not allow us to pursue in any detail. 
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To collect the minimum of formal relations we will need to refer to, we temporarily 
ignore the fact that we actually have to deal with non-abelian gauge-theories (these 
being the only ones for which our basic premise of asymptotic freedom is correct in 
the physical dimension D = 4) and write simplified equations for some generic field 
ip{x), using notations as in eq. (|1.25| ) above. The generating functional of Euclidean 
correlation functions, 

Ze[J] = Jl J I^Me-^^[^l + (-^'^) , (3.47) 

with normalization M such that Ze [0] = 1, is postulated to be invariant under simul- 
taneous infinitesimal translations of all its integration variables. 

If (x) — > (p (x) = If (x) + 6 If (x) (all x) , (3.48) 

in the sense that 2? [92 ] = ^'[v']- (For the mathematically minded, it is not unprob- 
lematic to use translation invariance for an infinite-dimensional measure, but since the 
procedure leads quickly and elegantly to equations of motion that can also be establish- 
ed independently by operator techniques, we will accept it here as a formal generating 
tool.) This gives the statement that 

which may be written as 

5 Se [v>] 

6 ip{x) ' ^ V J 
using the expectation-value notation 

{A)j = ^ J P^AMe-^-I^l + W'^) (3.51) 

for functionals A of Further rewriting along familiar lines puts this into the form 

] f SSeI^pY 



+ J{x)) = , (3.50) 



6 if (x) 



J{x) } Ze[J] = , (3.52) 



a relation sometimes called the master DS equation. The equations for the 1-, 2-, 3- 
... point correlation functions then follow from this relation the same way that the 
functions themselves follow from ( p. 47 ): by functionally differentiating 0, 1, 2, ... times 
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with respect to J and putting J = 0. For example, after one differentiation with respect 
to J {y) one arrives at 

Uy)'-^)=5-{y--), (3.53) 

\ Oip[x) / 

and if 5^, as usual, consists of a bilinear kinetic term 

if {x) ip (x) (3.54) 
with symmetric, positive-definite differential operator K^^ plus tri- and quadrilinear 



interaction terms collected in a "potential" y((/3(x)), then ( p. 531) will be an integro- 
differential equation for the two-point function G2 (y, x^ of eq. ( |0.lD of the form 



G2 (y, x) = 6\y -x) - {if (y) V ((^ (x)))o , (3.55) 

coupling that function to three- and four-point functions of the special forms 

((^(y)((/p(x))2)o , {ip{y){ip{x)f)^ . (3.56) 

It is this characteristic coupling to the next higher functions that is referred to as hier- 
archical. Upon converting the equations from x space to momentum space, the iden- 
tifications of two or three points in the functions ( p. 56 ) will turn into single or double 



internal- momentum integrations (loops) . Finally, one transcribes the equations from full 
correlation functions Gn to connected, amputated, and 1 PI functions F^v to give them 
their most transparent form. 

In Fig. 3.1 I have attempted to write diagrammatically the three lowest levels - the 
equations for two-, three- and four-point vertices - of what might be called a generic 
DS system. (These three equations, for our purpose, will actually turn out to bear a 
crucial intrinsic distinction, being the equations for the superficially ultraviolet-divergent 
functions). For each equation. Fig. 3.1 selects a set of typical and recurrent features: in 
the first line of each, the r. h. s. starts out with an " inhomogeneous" term representing 
a zeroth-order perturbative vertex pWpe'"* (^^gj-j^^^g 2 A, 3 A, 4 A). Next on each r. h. s is 
a loop term containing no vertex function higher than the one on the 1. h. s., so that 
with such an interaction term alone the equation would close upon itself (terms 2B, 3B, 
4B). Spoiling this closure is, on the second line of each equation, a loop term featuring 
the next higher vertex function, which brings in the hierarchical coupling (terms 2C, 
3C, 4C). 
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Fig. 3.1 Beginnings of a schematic Dyson-Schwinger hierarchy in 
diagrammatic form. Small circles prepresent bare vertices; 
heavy lines, exact propagators; large circles, exact proper 
vertices as specified. Dashed line in r4 equation indicates 
limit beyond which loop terms do no more contribute to 
generation of zeroth quasi-perturbative order. 



Also present in most cases, but not made explicit, are a finite number of other loop 
terms, which in the case of non-abelian gauge-boson vertices will include two-loop terms 
arising from the second type of function in (|3.56D . By suitably interpreting the dressed 
propagators (heavy lines) and dressed vertices (large circles) in their explicit terms, and 
suitably filling in the non-explicit ones, these three graphical equations may be viewed 
as encompassing the lower levels of practically every known DS hierarchy. To guard 
against misunderstanding, I must interject a remark here. Because of the hierarchical 
coupling, DS systems, like their non-relativistic counterparts in statistical mechanics, 
are usually solved by "decoupling" approximations whereby the next-higher vertices in 
some equation are expressed in terms of lower ones, thus closing the system. Although 
these decoupling approximations are efficient devices for generating partial resumma- 
tions, and in my opinion do not deserve the low esteem in which they seem to be held by 
most contemporary particle theorists, it must be emphasized that what we are going to 
study in this section - the self-consistency of nonperturbatively modified zeroth orders 
in a quasi-perturbative resurgence reconstruction - has nothing to do with decoupling 
approximations. The one "decoupling" that will occur here will be both automatic and 
exact, and will arise from the basic property of renormalizable theories of having only 
a limited number of superficially divergent correlation functions. Decoupling approxi- 
mations may be useful at later stages, particularly when dealing with the equations for 
superficially convergent higher amplitudes of which our Fig. 3.1 stops short, but for the 
purpose at hand, they play no role whatsoever. 

The common form under which all the equations of Fig. 3.1 may be subsumed is 

= tP^'"-' + [T] {N = 2, 3, 4) , (3.57) 

which makes explicit the fact that all interaction terms carry at least one power of the 
squared bare coupling, g^. The DS interaction functional, consists of loop inte- 
grals over integrands involving all vertices F up to Fat + i (or F7V + 2 in the presence of 
quartic interactions), and is nonlinear, as is directly visible from the graphs of Fig. 3.1. 
Straightforward iteration of eqs. (3.57) around the zeroth perturbative orders is one way 
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of generating the perturbation series ( p.4[) , but in such an iteration, in which at every 
stage one evaluates integrals over functions known from previous steps, one never fully 



faces the true challenge of ( 3.57 ), which is a self- consistency problem: find a set of Tn 's 
such that, when inserted into the functionals, they are capable of reproducing them- 
selves. Nonperturbative treatments, along whatever lines, must address that challenge: 
to determine, in mathematical terms, a fixed point of the complicated mapping defined 
by the r. h. s. of eqs. (|3.57] ). In our specific context, not the least part of that challenge 
is to establish from eqs. ( |3.57] ) the nonperturbative parts of the zeroth-order input - the 
terms 

_ p(0)pe.t (, = 0,1,2,...) (3.58) 

which carry powers A" with n > 1 but no powers {g'^)^ with p > 1 - in spite of the 
fact that the interaction terms always have at least one power of Qq in front. For this 
to be possible, the interaction functionals must obviously be capable of producing 
terms with factors of 1 / (7^ from sufficiently nonperturbative input, and indeed that is 
what they will be seen to be doing below. (One might say that the basic "mistake" of 
perturbation theory is to miss precisely this possibility.) 

Besides a multitude of aspects of DS systems that will have to go unmentioned in the 
restricted time and space of these lectures (a few will be taken up in our closing section), 
there is one that is so rarely discussed in the literature, and so intricately connected 
to renormalizable ultraviolet-divergence structure, that I feel I should comment on it 
before plunging into the technicalities of self-consistency. In a theory with up to quartic 
terms in its action functional 5^; [(/?], the "master DS equation" ( 3.52| ), of which our 



individual equations (3.57) are functional moments, is a functional-differential equation 
of the third order for the generating functional Ze [J]- Prom the analogy with differential 
equations, one would therefore expect that the first three functional Taylor coefficients 
in the expansion 

^e[J] = ^ [ d^xi... I d^XNGN{xi...XN)J{xi)...J{xN) (3.59) 

will be left undetermined by the equation and will have to be specified as "initial con- 
ditions" (for the evolution of the functional away from J = 0). The = coefficient, 
the vacuum Norm Go = Ze [0] = (0|0), is fixed at unity by normalization, and the 
= 1 coefficient, the single-field expectation Gi = (0|<^(x)|0), vanishes in theories 
without elementary scalar fields (or indeed remains an empirical parameter in Higgs- 
field theories, but then that's what we have been accustomed to). But the idea that 
the N = 2 coefficient, the propagator G2 {xi — X2), should have to be put in as an 
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initial datum sounds positively bewildering. I must refer you to the excellent discus- 
sion in sect. 1.4 of ref. | RIV87 |, which uses the so-called static ultra-local model of 
Caianiello and Scarpetta | CAI74 |, for insight into how the perturbation series seems to 
evade this problem - basically, it steers you tacitly into a particular solution in which 
the initial datum G2 is assigned a special value, while displacing the freedom of choosing 
initial conditions into its undetermined "escaping" terms. (Again, analogs of this have 
been known in Ecalle's theory of resurgent solutions to nonlinear differential equations, 
where the perturbation series is referred to as the "formal solution"; see the "epilogue" 
of [ CNP93U ). In any truly nonperturbative treatment, however, the initial-conditions 
problem is bound to reappear, and indeed the numerical studies of propagator equations 
described in [AVS01| have met with indeterminacies presumably related to this question. 
In our context, the problem will resurface in a very clear-cut and interesting way: while 
the standard DS system will produce self-consistency conditions determining the poly- 
nomials (A;^, A'^) of eq. ( ^.11 ) above, which at each level r encode the particle-pole 
positions of the propagator and the mesh points for the discretization of its continuum 

integral through their root-factor decompositions ( 3.20| ), conditions that would deter- 

r 1 ^ ' 

mine the coefficients of the polynomials (^^, A^) of eqs. ( 3.13| ) and ( p. 27 ) - that 
is, the weights to be assigned to those poles and mesh points - will be found lacking. 
When reflecting on the origins of this indeterminacy, both in DS systems in general 
and in our more special problem of zeroth quasi-perturbative orders, we are led back 
to looking at the OPE in e. g. the form of eq. ( ^.54 ). We noted that the coefficients 
appearing there contain information about renormalized vacuum expectations of 
composite operators of mass dimension n - singular products of up to n elementary 
fields at the same space-time point. In the ordinary DS system, the loops do represent 
some of these composites - the quadratic and cubic composites of eq. ( 3.56| ) - but not 
more, while the OPE, when following it to higher n's, obviously carries more informa- 
tion about higher composites. The natural conjecture then is that the indeterminacy 
revealed by the initial-conditions argument is related to insufficient composite-operator 
information in the ordinary DS equations. To supply the missing information in a way 
that is compatible with, and forms a natural extension of, the ordinary DS system, one 
may therefore think of using DS equations with increasing numbers of coalescing points 
in X space, which in momentum space would correspond to integrating out increasing 
numbers of the independent momenta. On the simplest level, for example, the y ^ x 
limit of eq. ( 3.53 ) would give 



(p{x) 



6 Se [(p] \ 







(3.60) 



(remembering that 6^ (0) 



6ip{x) / 

in dimensional regularization) , which in momentum space 



108 



will turn into an equation for the integrated quantity 



27r 



^(O)pert(^) (3.61) 



rather than for D {k) itself. (This equation is a single-number condition, and at the order 
r = 1 of rational approximation, where only a single coefficient Ug^' = rj^^ needs to be 
determined in the polynomial Ag^', it would already suffice for completing the "initial" 
conditions). On the next level, coalescence of two of the three spacetime points in a 
three-point equation would lead to an equation for the quantity 

^3 f-p, \v + (i,\p - q] , (3.62) 



2TrW \ " 2 



which is still an equation for a function of p, and coalescing all three points would give 
the integral of this equation over p, which like the equation for ( 3.61| ) is a single-number 



condition. Equations of this kind have been termed " equation-of-motion condensates" 
m some OPE studies [ |AHL92|| , and have b een used as checks on Wilson-coefficient cal- 



culations, but in our context they stand to assume a more fundamental role. In ordinary 
integral-equations theory with well-behaved functions and convergent integrals, once you 
have solved an equation, the integration of that equation over its domain of definition 
will not, of course, generate new information. In a theory with renormalizable ultra- 
violet divergences, the situation is different - composite-operator divergences require 
new renormalization steps beyond the usual field and parameter renormalizations, and 
the above-described momentum-integrated DS equations therefore do bring in new in- 
formation. In these lectures we will take the view that they are an adequate means of 
completing the initial conditions. 

When taking up self-consistency in earnest, it is again best to first concentrate on 
a partial task in which the basic mechanisms can be seen at work in an exemplary way. 
Let us therefore focus on the first, inverse-propagator equation of Fig. 3.1, and let us 
single out, among the interaction terms on its r. h. s., the one-DS-loop diagram (2C), 
which is the simplest one to bring in the crucial hierarchical coupling. We take the 
D^^ on the 1. h. s. to be a scalar function of /c^, k being the external momentum 
flowing through the equation, with perturbative limit (dC^)?*^^*)"! given by ( |3.6D , and 
with zeroth quasi-perturbative order given, at level r of rational continuation, by ( 3.11] ) 
or ( |3.28|) . The equation reads, in the notation of ( |3.57|) , 

D-i (A;2) = k'^ -gl {$2c [r; k^) + [other loop terms]} , (3.63) 
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where, remembering our parenthetical remark after eq. ( 3.31 ) 

^ ^^""^ (3.64) 
X D{q^)D{{k - qf)T^{q - k, k, -q) . 

In order to get mass dimensions (and the associated degrees of divergence) right without 
having to delve again into the paraphernalia of gauge-theory tensor structure, we will 
patch up things a little and supply both and p^'^)*"^''* with an extra unit of momentum 
dimension in an ad hoc fashion, such that their product develops an extra contraction 
factor 

[{k -q) - qr i-q - (q- k)r = {2q - kf . (3.65) 

(Never mind this little dose of patchwork - it will do all the essential things right. If 
you feel uncomfortable with it, you may do an instructive piece of homework by setting 
up, say, the gluon-loop contribution to the gluon's transverse D^^ as defined by ( 3.14| ), 



with full tensor structure for r3_g;„o„ as in eq. ( 3.32| ), working in Landau gauge to avoid 



mixing in Z?l's, and grind out the tensor contractions. You will find that our ad hoc 
factor does things qualitatively right for the ultraviolet-divergent terms, which turn out 
to be what matters here). 

In a first iteration, we now substitute D and by their zeroth quasi-perturbative 
orders, using eq. ( 3.4C| ) at level r of rational continuation, with perturbative limits given 



by ( |3.6D and ( |3.3lD . That is, we decide - in order to keep all inessential complications out 
of the following discussion of the self-consistency procedure - to only look for solutions 
with 

Ct^ = 1 , it^ = 1, (3.66) 

which represents a kind of boundary condition. Indeed we remember from sect. 2.4 that 
there may be at least two different meanings, related by finite rescaling, to the concept 
of quasi-perturbative zeroth-order, and we expect that the DS functionals in principle 
can produce both (and perhaps more), so a supplemental condition of some kind will be 
necessary. Eq. ( |3.63| ) assumes the form 

(3.67) 

= k'^ - {$2C [r[^'°); k^) + [other loop terms]!^'^) + O {g"^)} , 



110 



where 




(3.68) 

iVf^ (<? - g2. A2) I 1 



At^J {{q - fc)2, A2) A![' A2) J A^^J (A;2, A2) 

At this point, to shorten the discussion, we anticipate a simphfying fact that will soon 
be obvious: the rational-approximant poles at k'^ = — Ug'^j, A^ (s = l...r), which ac- 



cording to eq. (3.11) we are postulating for the 1. h. s. of eq. ( |3.67| ), can at zeroth 



quasi-perturbative order arise only from the last factor 1 / A3 of ( |3.68| ) but not from 
the bracketed integral. Therefore we immediately have a first set of self-consistency 



conditions, which in the notation of eq. (3.44) read 

rM=41 is = l,...r), (3.69) 

and which imply that quite simply 

aM (A;2) = aM (A;2) (all r) . (3.70) 

We remark, without elaborating, that this situation will repeat itself when considering 
the three-point equation, to the effect that 

A^ (A;2) = aM (A;2) (all r) , (3.71) 

where A^ is the polynomial appearing in denominators of level-r approximants to the 
invariant functions of the 4-point vertex. In other words, the vertex-denominator (or 
propagator-numerator) polynomials Ag , and the vertex-function branch cuts which they 
represent, turn out to be a common element in all zeroth-order, basic vertex functions. 
Depending on viewpoint, this element is "handed up" or "handed down" the ladder of 
DS equations by the hierarchical coupling, which is seen here playing a prominent role. 

For our immediate task, this brings substantial simplification since A's now cancel 
in the integrand of ( p. 68 ). For the remaining fraction, it is straightforward to write a 



partial-fraction decomposition, first with respect to the loop variables q^ and (k — q)"^ 
and then also with respect to the external k"^ . The former step reduces ( |3.68| ) to a linear 
combination of the basis integrals 

(3.72) 
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These may be evaluated in D = 4 — 2e by the standard techniques, the result being 



1 



(4vr) 



vlM „ A2; s)-k' + Bt\nik', A'; e) ■ A^l . (3.73) 



The functions ^m,n and Bm\n are given by the Feynman-parameter integrals 



AM,(fe^ A^;^) = / dz 




(3.74) 



X <^ [12z(l - z) - 1] 



1 



7 + In 4 vr + 4 z {I - z) + O (e) 
SW„(A;2, A2;z) 



X ^8[zctM + (1 - .)aM] 



1 1 

7 + ln47r + - 



(3.75) 



+ 0(6) 



and near D = A have weak, logarithmic dependences on k'^ and A^ through the (— e)-th 
power of the quantity 



(fc2, A2; ^) = ^(1 _ ^)A;2 + L^;] + (1 _ 



a2 . 



(3.76) 



(To maintain correct mass dimensions, we have stopped short of expanding the non- 
integer power). After collecting the debris of the partial-fraction decompositions, eq. 



( 3.67 ) now assumes the form 

[D-^ (P)][r,0) ^ g2 [^-1 (^2)][r,l) ^ q (^4) 



9o 



^ + { ) { ^0 ' -k' + rl'J (A:^ A^; e) • A^ 



(3.77) 



'2s + l 



(A:2, A2; e)A4 



A;2 + n2 s A2 



+ [other loop terms] + O {g 



\r] 

where the quantities have structure 

r + 1 



5^ [a!:L,„ • A^Anik', A'; e) + b^^.n " ^Hn (fc', A^; 6)] , (3.78) 



m,n = l 
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and therefore weak k"^ and dependence near D = A. About the dimensionless coef- 

[rl \r\ 

ficients a\ ^ ^ and 6- |^ we need to know no more than that they depend hnearly on 
the level-r, three-point vertex coefficients x^, and in general nonhnearly on the set of 
parameters of the inverse propagator itself (the of eq. ( 3.28| ), or any equivalent set). 



With eq. ( |3.77| ), we are ready for posing the self-consistency question. It is, by the way, 
no later than at this stage that one realizes that condition ( p. 70 ) is a necessity. 



This is, once more, a moment to pause - two observations are worth mentioning. 
First, while our integrand, with its A^ terms all over the place, is manifestly a strongly 
nonperturbative object which has no power-series expansion in the coupling at all, noth- 
ing has kept us from computing the loop by the standard techniques. These techniques 
are sometimes, misleadingly, referred to as '"perturbative" , which sounds as if they were 
applicable only in pure perturbation-theory calculations, but in fact they depend, of 
course, only on the rational structure of integrands with respect to momenta. 

Second, it is clear that from a purely technical standpoint, nothing would have stood 
in the way of our using directly the resurgent symbols - that is, the operator-product 
expansions - of D and of in the loop integral ( 3.64 ). Instead of ( 3.72| ) we would have 
encountered momentum integrals of type 



oc , (3.79) 



(2 7r)^ [{q - kYY[q^]' \k^ 

which are increasingly singular as k"^ 0, and we might have proceeded to matching 
the series involving these integrals to the series of eq. ( 2.56| ) at increasing orders 



n, 



( R) 

thus determining in principle the resurgent-symbol coefficients Un (0, 0). But it would 
nevertheless have been foolish to do so. On the mathematical side, the coefficients of 
proportionality in ( p. 791 ) would be Feynman-parameter integrals increasingly divergent 
for s or t > 2, and would force us to throw in a new set of rather obscure regularizations 
that might (hopefully) go away in the symbols's final resummation. Even if we were will- 
ing to adopt such an ugly formalism, after all the effort we would still be around with 
our resurgent symbol, a large-A:^ asymptotic representation, and would have made no 
progress at all toward knowing the function (D^^) at intermediate and small k"^ ~ just as 
in Borel-Laplace reconstruction when we forget the crucial continuation step ( |1.51| ). On 
the physical side, we would have decomposed the well-defined physical effect described 
by our diagram Fig. 3.1 (2C) - the self-dressing of a field quantum by virtual pairs 
of quanta of its own kind - in a totally unreasonable way: into a sum of increasingly 
infrared-singular terms with increasingly violent small-/c^ cancellations between them. 
Finally, this procedure would have broken down completely for the "tadpole" diagram 
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of Fig. 3.1 (2B): by using the resurgent symbol for its single-D integrand, we would 
have decomposed it into a sum of scale-free integrals which vanish in dimensional reg- 
ularization, whereas with any reasonable analytic continuation this term contributes a 
nonzero multiple of A^, in marked contrast to perturbation theory. Thus it is abun- 
dantly clear that (as mentioned already in the introduction) we cannot simply plug our 
resurgent symbols into the dynamical equations in order to determine their parameters 
- some amount of analytic continuation is necessary from the outset to deal with objects 
globally usable in momentum space, and it is such objects that we are trying to establish. 



Returning to eq. (3.77), we note that although our task here is more intricate than 



the standard task of matching a rational approximant to an explicitly known function, 
the two problems seem to share the property that there is in principle a continuum of 
matching prescriptions for determining the approximant 's coefficients - Pade I or II, 



as mentioned after our eq. ( 1.53 ) above, are just two among many (and Pade II in 
itself already represents a continuum of possibilities). But a moment's reflection will 
show that in the present problem there is a distinctly preferred matching strategy: we 
should clearly see to it that we get the poles (at the zeroes k'^ = — i -^^ ^2 
denominator) right, from the r. h. s. of ( |3.77| ), for the postulated zeroth-order {D~^)^'^'^h 



If we failed to ensure this, the remainder in eq. (|3?r7| ), which is to give (the (2 C)-type 



contribution to) the first-order-quasi-perturbative correction (D^"*^ )''"'*'•*, would have its 
typical one-loop analyticity properties contaminated by splinters from the zeroth-order 
poles, and thus would be plagued by approximation errors that become locally infinite in 
the fc^ plane. Having already ensured correct pole positions through conditions (|3.6S| ), 
we therefore need to match residues at those positions. Looking at eq. ( 3.2 j ) (with 
Q = 1) above, we are led to imposing the self-consistency conditions 



41 + 1 = (|^)^rl + l(-2.A^ A2;e) {s = l...r) . (3.80) 

For the ri A^ term, there is no such obviously preferred matching prescription, so here 
we may freely choose a matching point k'^ = —uji A^: 

= (^)^rn-^lA^A^^) (3.81) 

(This freedom of choice is actually part of the freedom of performing finite renormal- 
izations, since the ui term can be shifted around by a local Lagrangian counterterm) . 
In all these conditions, the anxious question of how's-this-going-to-give-a-zeroth-order, 
with the explicit factor in front, is still hanging over us. 
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But only apparently! With k"^, in these conditions, now nailed to fixed multiples of 
A^, the Am. n and Bm_n functions entering in eq. ( 3.78 ) now appear with values 



m, n 

^m,n ("cA'^, A^; e) 



(A2 



1 



--{l + [0(e)](c)} 



bM„(-cA2, A2;e) = (A2)-1 {4(cT^ + + [0{e)]{c)} 



(3.82) 



(3.83) 



so all the rj*^^ at matching points — c A^ in conditions ( p. 80 ) and ( p. 81 ) are found to have 
factors of (A^)""^ j up front. Now take a closer look at the quantity 



- \ 2 
in 



(A2 



1 



47r 



A2 



(3.84) 



where relation (p. 61) has been used. Proceed the way you are accustomed to, Laurent- 
expanding the last two factors around e = 0. You have 



A2 

/i2 



In 



A2\ 1 
^ 1 + - + 0(e) 



(3.85) 



The logarithmic term survives the e — > limit in the same, familiar way as the well- 
known In (/c^ / /x^) 's in a perturbative calculation, but it differs from them in one impor- 
tant respect: it leads them by one order in the quasi-perturbative {g'^Y classification. 
From eq. ( p. 25 ), 

,2\ / . _ \ 2 



In 



47r 



1 



g (^) J Po 



[l + Oig' In 



(3.86) 



at the one- loop level. This is precisely the 1 / g'^ factor needed to "eat" the overall g^ 
factor in ( |3.84 ) and make zeroth-order-in-gr^ terms possible. There's more: we now have 



- \ 2 

in 



(A^)- 



1 

Wo 



1 + Po 



9 

in 



2 1 



+ 



(3.87) 



but Za, to the same order, is given by eq. ( ^.13|) , so that to the order at which we 
calculate here, the ultraviolet divergence too gets "eaten" in the special quantities that 
enter our matching conditions. From eq. ( ^.80 ), this in particular means that no nonlocal 
counterterms will be needed for the pole terms in eq. ( 3.77|) : as we expected from our 
discussion of the OPE, they establish themselves in a finite manner. 



Our simple one-loop-level argument has established the values 

^M„(-cA2,A2;0) = 



(3.88) 
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4 



[r] 



(3.89) 



at = 4 only up to terms O {g^ In g'^), which stih seem to be a potential nuisance 
in conditions ( ^.80 ) through ( ^.811) . It is therefore useful to know that a more careful 
argument [ 3TI96| establishes these values exactly to all orders. This argument avoids 
the use of the unrelated expansions ( |3.85 ) and ( ^.13 ) completely, and instead employs 
the exact integral representations for the two middle factors on the r. h. s. of eq. (|3.84| ). 
For Za, we already noted t'Hooft's representation, eq. ( 2.12| ) above. For (A^)~^ - be 
careful here to consider the e > quantity formed with the full e > beta function 
( |2.8| ) - it is a matter of a few lines of calculation to deduce from definition ( |2.20| ) the 
exact representation 



(9/4^)2 



9' (^) 



exp 



(si/ 4 77)2 



(3.90) 



involving again the function X ('^) defined by ( 2.11 ). Upon combining the two represen- 
tations one sees that in the quantity ( 3.84 ), the dependence on g^ (/x), and therefore on ^u. 



cancels exactly (remember gi is a trivially RG-invariant integration constant). Further 
analysis of the integral in 



9l 



(91 / 4^)2 



(47r)2£ 



exp 



dK 



(3.91) 



by exactly the same technique as applied in eqs. ( ^l5l) and (|2l|) gives, finally, 

2 



(A^)- i = 1 [1 + O (e, e In e)] , 

e po 



(3.92) 



rigorously independent of g^ {ji). Despite the superficial similarity, this should not be 
confused with eq. ( 2.16| ); it says something conceptually quite different. While ( 2.16 ) is 



a representation for with a g^ (^) dependence of the r. h. s. that is difficult to control 



as long as one does not know the exact RG beta function, the result (3.92) states the 
presence of an exact RG invariant, which, moreover, is finite and scheme-independent 
at e = 0. Its point is the exact contragredience of the g"^ dependences in the quantities 
(2.12) and ( p. 901 ) already ate^Q. Both a g^ factor and an ^ divergence get "eaten" in 
this remarkable combination. 



116 



It is not necessary, in these seminars, for us to pursue the consequence of eqs. (13. 8q ) 
and ( 3.89| ) through the rf^ functions to the self-consistency conditions ( |3.80|) / (|3.81|) in 



every detail. It is sufficient to state that the latter now assume, at £ = 0, the general 
form 



2S + 1 (WM, {x}M) {s = 0...r) , (3.93) 

with functions i?2s + i depending linearly on the set {x}''"' of three-point vertex coeffi- 
cients, and in general nonlinearly on the set {u}''^! of the inverse-propagator parameters 
Ui themselves. These conditions allow one, remarkably, to establish functions (D^^)^"^'^^ 
of zeroth quasi-perturbative order from DS loops with a prefactor, and to do so in an 
ultraviolet-finite manner not requiring nonlocal counterterms. 



As we have anticipated, apart from equations such as ( p.70| ) and ( p. 71 ) that establish 



[rl 

a common set of denominator coefficients Ugl, for all vertex functions, we have not found 
conditions that would really determine these coefficients. We have already sketched how 
one may go about filling this gap by invoking higher degrees of compositeness; setting 
up equations for the quantities ( 3.61 ) and ( p. 621) one obtains the r additional conditions 



needed. Apart from the technical complication that such equations now have at least 
two-loop terms on their right-hand sides, and involve the trickery of extracting two-loop 
divergences, their treatment requires no truly new ideas; in particular, relation ( 3.92| ) 



again plays a central role. We therefore skip details of these "integrated" conditions 
entirely. For the same reasons, we do not enter into details of the next higher ver- 
tex equations, which involve application of the same technology to more variables and 
lengthier loop integrals. For example, the middle (three-point) equation of Fig. 3.1, 
when analyzed for terms of zeroth quasi-perturbative order, will produce conditions for 
the three-point parameters {x}'^' of the general form 

•^mim2m3 ~ ^ 'S'mi ?n2 ma (^{"u}' \ {x}^ \ {t}^ , (3.94) 

with right-hand sides now involving the set {t}^^^ of four-point numerator coefficients. 
The structure of these conditions, which after multiplication with enough denominators 
become a system of nonlinear algebraic equations, mirrors quite closely the structure of 
the underlying DS system with its nonlinearities and hierarchical couplings. 



As for the remainder terms in eq. ( 3.77] ) that we expect to yield the first quasi- 
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perturbative correction, we are now left with 



A' 



^ rM^,(A:^A2;e)-rM^,(-^2.A^A^;e) 



'2s- 



(3.95) 



s = l 



+ [contributions from other loop terms] 

The unsubtracted term r^^ k'^ carries the A^ — > " 'perturbative limit" as it would 
come out from the usual perturbative calculation of the loop, including the perturbative 
i ultraviolet divergence. The other terms, which the self-consistency process has left 
subtracted, are both ultraviolet finite and (for the terms with A^) regular at k"^ = 
— n2sA^, since by eqs. ( 3.7^ ), ( p. 74 ) and ( |3.75| ) they are linear combinations of the 
quantities 



m, n 



SM (-cA2 a2- 



-£{k^ + cA2)z(1 



1 

wzil 



(3.96) 



+ 



z a. 



M + (l-z)aM-(l-t«)cz(l 



A^ 



With respect to renormalization, we are therefore led into precisely the situation we 
expected from our OPE discussion around eq. ( p. 49 ): we have to follow the "two-track" 
procedure of supplying an explicit counterterm, proportional to (7^ ^, for the perturbative 
portion, while the terms with powers of A^ establish themselves in an ultraviolet-finite 
manner. As for the in expression ( |3.95D , we note that its rigorously (^^-independent 
combination ( |3.92| ) with the quantity (A^)""^ ^ no more occurs, and since we are working 
within a quasi-perturbative scheme, we must therefore revert to the usual pattern of 
treating it as ^i'^^ {pi) + O {g'^), the serving as usual to render the (SM)~^ terms 
in the quantities (|3.74D and ( |3.75| ) dimensionless. 
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I is clear that our condition ( ^.66 ) of maintaining strict perturbative limits bars 
us from " oversubtracting" our result ( ^.77 ), in the sense of subtracting also the r^^ 



term at some point k'^ = — wqA^ and invoking the mechanism of eq. ( 3.92 ). Such a 
procedure would inevitably lead us away from the perturbative limit = 1 and into 
a different version of the quasiperturbative expansion, analogous to the one we noted 



for our model amplitude in eqs. (2.212) and (2.214). We are not going to explore this 
route here, interesting as it may be, since it would also force us to reorganize the entire 
perturbative-renormalization procedure, and therefore would lead us too far away from 
our stated focus. 

Several features of this self-consistency mechanism are worth commenting on. First, 
the final self-consistency conditions at I? = 4 (i. e. e = 0) are independent of the chosen 
subtraction points —U2s and —loi A^, since the latter only enter into the terms O (e) 
in eqs. ( ^.82|) and ( 3.83| ). The subtraction points do, however, enter into the first-order. 



quasi-perturbative correction [D ^]['''^) of eq. ( 3.95 ), and are essential for keeping it 



uncontaminated by the zeroth-order rational-approximation poles. 

Second, the conditions at D = 4, while derived here in a dimensional-regularization 
scheme, make no reference any more to a special renormalization scheme, since Pq is 



scheme-independent. Now different schemes, as a noted in eq. ( 2.23 ) above, have A 
scales differing by constant numerical factors. This leads us to suspect that upon scaling 
A by a factor, the conditions will remain invariant when scaling all the nonperturbative 
X and u coefficients inversely to the A^ powers with which they are associated. This 
scaling property indeed turns out to be present in the system, and has the effect that 
the conditions in fact determine only the ratios of all coefficients to one of them. The 
undetermined coefficient, with its associated A^ power, effectively defines a modified scale 
A, which can be fixed only by calculating and matching some dimensionful observable. 

Third, it is manifest that the entire mechanism is intimately tied to the presence 



of the ultraviolet loop divergences: the crucial term ( 0.86 ), and the ensuing terms ( 3.8^ ) 



and ( |3.89| ), while themselves finite, have survived the removal-of-regulator limit e ^ 
only due to their association with a ^ divergence. We have all known such an effect from 
our experience with perturbative logarithms, but here it assumes a new and much more 
prominent role. The self-reproduction of zeroth quasi-perturbative orders is therefore not 
only a genuine quantum effect (as it comes from loops) but also a genuine quantum-field 
effect (as it comes from divergent loops). 

This has the immediate consequence that, fourth, the formation of zeroth quasi- 
perturbative orders becomes rigorously restricted to the small finite number of superfi- 
cially divergent vertices - by and large, the same functions that already possess a zeroth 
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perturbative order. You need no big calculation to see this. In a superficially conver- 
gent amplitude, the one-loop contributions have factors of Qq but are, by definition, 
genuinely convergent, so they cannot trigger the mechanism of eq. (|3.92| ) and remain 
quasi-perturbative corrections of the first order. The two-loop contributions, with factors 
of gQ, may exhibit ultraviolet divergences, but only in the form of one-loop subdiver- 
gences producing a single j factor, and generally in the higher /-loop terms the number 
of i factors produced lags behind the number of overall factors by at least one unit. 
So the mechanism of eq. ( p. 92 ) can "eat" at most I — 1 oi the I factors of g^, and all 



contributions remain p > 1 quasi-perturbative corrections. 



The same consideration applies to terms such as (4C) of Fig. 3.1 that couple the 
highest superficially divergent to the lowest superficially convergent vertices. Since there 
is no Y^^^'^^^ in a renormalizable action, the lowest contributions to term (4C) have 
replaced by its convergent one-loop terms and two factors of g^ in front, of which at 
most one can be eaten. Thus as far as self-reproduction of zeroth orders is concerned, 
the hierarchical coupling automatically stops acting after the four-boson vertex equation. 
Once more we have a fine line (depicted in the third equation of Fig. 3.1) marking a 
fundamental divide, this time within the strictly renormalizable theory itself: between 
the equations for superficially divergent and for superficially convergent functions. 

The self-consistency conditions, fifth, do nevertheless retain an infinite-expansion 
aspect of a different kind, at least in the somewhat conservative treatment we have af- 
forded them here. We have been classifying the contributions to the right-hand sides of 
the integral equations in the time-honored way by their number of loops, or indepen- 
dent momentum integrations, that arise upon iterating the equations around the newly 
established, nonperturbative starting points r'^''^^. It is clear that this will translate 
into an analogous loop classification of contributions to the right-hand sides of the self- 
consistency conditions. In a superficially divergent amplitude, there is no lagging of the 
number of g^ prefactors; an /-loop diagram here produces terms of type 

^..^.i} ,3.97) 

and through the mechanism of eq. ( |3.92D therefore gives contributions to the quasi- 
perturbative series of the forms 

Po Pq Po 

all the way down to the zeroth order: a given loop number / produces not only terms of 
order (g^)' but also corrects all previous orders (g^)^ with p < I. That the tight coupling 
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between I and p orders characteristic of a perturbative single-coupling theory should get 
loosened is presumably something that happens in any truly nonperturbative treatment. 
(It complicates diagrammatic argument, in that a quasi-perturbative correction T^'^'P^ 
with p > 1 can no more simply be written as a set of p-loop diagrams, but must include 
finite counterterms to account for the pieces that go into orders p < p). But here it 
has the specific consequence that the self-consistency condition for a generic coefficient 
of a zeroth-order nonperturbative approximant in principle gets corrections from all 
higher loop orders that form a series in powers of 1 / Pq: 



This type of series will not create new summability problems of its own, since from ( 3.97] ) 



you see that it arises from the "most strongly divergent terms", terms proportional to 
(1 / eY , in all loop orders I, and these terms, like the "leading logarithms" to which they 
are closely related, generally do not exhibit factorial divergence but form a series with 
nonzero convergence radius, behaving roughly like a geometric series. But I have not yet 
had enough experience with this formalism to have a ready suggestion for its practical 
resummation, which at this time remains a (somewhat technical but not unimportant) 
research problem. In pure SU (N) Yang-Mills theory, where po = ^N, eq. (|3^ ) 
clearly is of the nature of an expansion (though only formally, since may well 
contain contributions comparable to (/9o)' of its own), and since general experience with 
expansions indicates they tend to work much better than anyone has a right to expect, 
it seems not unreasonable to start out by drawing conclusions from the lowest orders of 



(I3.99D only. 



Sixth and last for these comments, while we have stopped short of technical detail 
for three-point and four-point vertex equations (on that subject I must refer you to refs. 
| DRI99| ]), it is already clear from the inverse-propagator equation alone, in e. g. the 



form ( p. 77 ), that the mechanism allows us, in the well-defined approximation of a zeroth 



quasi-perturbative order, to watch the spectacle of the spontaneous generation of mass 
for a field that was massless on the Lagrangian level - the process we referred to in eq. 
( 2.104| ). If the field decides to form a single massive particle, its squared mass at level r 



will be given by the first root of the propagator denominator ( |3.2Cl| ), 

(m2)W=aflA2, (3.100) 

Ir] 

with a[ calculable from the self-consistent Ui coefficients. Because of the scaling prop- 
erty, solution of the [r, 0) self-consistency problem alone does not, of course, make this 
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an absolute prediction, since at this stage you do not yet know which A scale you are 
working with. But ratios of particle masses do in principle become predictable. It is 
interesting to note that in some sense this is already true on the r = level (which is 
not to be confused with the perturbative limit). The r = approximation is so crude 



that it can be accused of internal inconsistency ~ as you see from (3.41), it keeps the 
three-point, and in fact also the four-point, vertex functions at their perturbative forms 
-p^o^peri^ at least under our "boundary condition" (3.66) above, but it does allow the 



inverse propagator a nonperturbative mass term proportional to A , 

so that With such a crudely restricted scheme, you are not surprised when 

you work out the self-consistency condition and find the form uf^ = const. • u^^^ which 
allows only the trivial solution u^^^ = 0, since for realistic values of /?o the constant 
refuses to be unity. But it is interesting that for an asymptotically free theory with 
Pq > 0, the constant is positive, so that even at this oversimplified level the system 
sends a muted signal of its tendency to produce a spontaneous mass. 



3.3 Things done and not done 

At the end of these seminars, let us take a stroll along the periphery of the landscape 
that we tried to begin to explore, casting glances at the few things than have been done 
or attempted in the matter of construction of coupling-resurgent amplitudes, and at the 
larger number of things that have not been done but could be worthwhile taking up if 
you find these questions interesting. 

The method described in the last two sections for constructing resurgent vertex 
functions from quasi-perturbative series, based on a set of generalized Feynman rules, 
has not at this time been tested to a degree that one could call comprehensive, but there 
has been a reasonably complete application to massless QCD at r = 1, the lowest non- 



trivial level of rational continuation, for which I must refer you to the articles [DRI99|. 
This study employs Landau gauge fixing and plausibly motivated restrictions on the 
oppressively complicated tensor structure of the four-gluon vertex, and concentrates on 
the case of at most two flavors of massless quarks. If finds that the two superficially 
divergent vertices for Faddeev-Popov ghosts, at the zeroth quasi-perturbative order, re- 
main perturbative in Landau gauge fixing, a conclusion I would like to see corroborated 
by independent studies. For the zeroth-order self-consistency problem of the remaining 
five superficially divergent functions (three for the gluon and two for the fermion sector), 
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it finds acceptable solutions with all-real approximants, which after throwing out a few 
on physical grounds are even remarkably unique. Unlike the pure Yang-Mills case (no 
quarks) which is also considered, the r = 1 transverse-gluon propagator obtained is not 
entirely satisfactory, since its zero at k'^ = — M2^^ slips a bit into the Euclidean region, 
thus violating reflection positivity on a small Euclidean interval - an effect coming from 
a strong coupling to virtual quark-antiquark pairs, which the approximation obviously 
exaggerates. You will tolerate my not entering into further discussion of the propagators 
obtained in those articles - they are interesting, but their interpretation would drag us 
into issues of confinement far removed from the focus of these talks. Here I am quoting 
refs. I DRI99[ | mainly as a demonstration of technical feasibility. 



One problem for which refs. [ DRI99| do offer a reasonably complete solution, even 



for arbitrary r, looks at first like a purely technical one. But it turns out to have a bearing 
on the fundamental issue of unitarity (total-probability conservation) , and therefore let 
me briefly stop at it on our stroll. When we found we had to impose eq. ( 3.70| ) as a self- 



consistency condition, you may have wondered in passing what would have happened 



in the diagram of Fig. 3.1 (2C), or its analytic expression ( 3.6§| ) at level r, without 
that condition. Then we would have had on each of the two internal lines, in addition 
to the r-level propagator poles of eq. ( |3.21 ) that provide a discretized Kallen-Lehmann 



description of the propagation of the virtual particle, r additional poles from an 1 / A3 
factor that would plainly contradict the postulated propagator approximant. Condition 



( 3.70 ) saved us from that embarrassment, but only because in that diagram both internal 
lines end in a bare vertex on the left-hand side. But now take a look at the second of 
the DS equations of Fig. 3.1, the three-point vertex equation, and its loop term (3B). 
In that loop, the right-hand part involving the vertical internal line represents a dressed 
one-particle exchange of the form 

r3(..., ...,Q)D{Q^)T^{Q, ...,...) , (3.102) 

with Q the momentum of the vertical line. Since that line now extends between two 
dressed vertices, only one of the two 1 / A3 's from these vertices gets cancelled by the 
A2 numerator of D {Q"^) by virtue of eq. ( p. 701) , and now the embarrassment of the 
"superfluous poles" seems to be really on us: such poles will, by the diagram-cutting 
rules for unitarity relations, produce contributions to these relations that do not corre- 
spond to any accessible external states, so unitarity would be endangered. It is a minor 
surprise (at least to me) to see that the formalism solves that problem automatically by 
its own devices: one finds that the four-point vertex r4 in diagram (3 C) unavoidably 
develops poles in its crossed channels that are the exact opposites of the "superfluous" 
poles in (|3.102 ). (That conclusion, remarkably, follows from residue-taking operations 
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in the three-point and two-point equations alone, that is, from lower equations in the 
hierarchy). Note that since the denominators in these "compensating poles" are not 
propagators of any of the elementary fields of the theory, their presence does not con- 
tradict the one-particle-irreducible nature of Their effect in the present context is 
to turn the vertical internal line described by (3.102) into what refs. [ DRI99|| denote by 



a dotted internal line, that is, a partial-fraction decomposition of expression ( |3.102| ) in 
which only the "true" propagator poles have been kept. Further analysis quickly reveals 
that these dotted lines really are almost everywhere, being present on every internal line 
between two dressed vertices. They are an important ingredient for keeping the formal- 
ism in line with the unitarity requirement. One amusing consequence of their presence 
that has not even begun to be exploited is that they provide a natural discretization of 
Bethe-Salpeter bound-state equations that refines itself with increasing r. 

What is entirely lacking at present - and that brings us to the part of the countryside 
where the things not done are lying around - is tests of the formalism at levels r > 1 
of rational continuation. (They are necessary not only for numerical refinement but also 
because rational approximants have a reputation for sometimes behaving erratically at 
low orders). The reasons are, to a large extent, purely technical and have little to do 
with physics: the nonlinear algebraic systems of self-consistency equations turn out to 



be a surprisingly tough assignment. Refs. [DRI99] were still able to get away with ad 



hoc methods of solution, exploiting accidental simplifying features of their Landau-gauge, 
r = 1 system, but higher r 's call for more systematic methods. From present experience, 
it seems fair to say that no presently known algorithm or available routine is capable of 
locating the rapidly growing number of solutions of realistic r > 2 systems completely in 
tolerable amounts of computing time, and most methods even have a hard time coping 
with r = 1 systems. So here the formalism also creates a new challenge for applied 
mathematics. The provisional solution currently available is to rely on minimization 
methods that do not guarantee to find all solutions, to let them search a reasonable 
portion of the space of starting values as finely as one can afford, and to rely on the 
experience that a fully real and physically acceptable solution, as a rule, stands out quite 
conspicuously and essentially uniquely. 

The correlation functions of QCD, by the way, may not represent the optimal prov- 
ing ground for first tests of the method described here ~ or, for that matter, of any 
method you may invent for coupling-resurgent reconstruction. The fact that the ele- 
mentary QCD fields have no asymptotically detectable particles associated with them 
tends to intermix and obscure all questions of nonperturbative methodology, in which 
we are interested here, with questions of the physical interpretation of various features 
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of the functions obtained - particularly of the propagators. A much cleaner environ- 
ment conceptually is the SU (2) gauge system of electroweak theory, with coupling to 
fermions but without strong interactions. Here, quasi-perturbative corrections are al- 
most certainly small at energies around the relevant A scale, and resurgent amplitudes 
of the zeroth quasi perturbative order should allow one in principle to explore with some 
amount of systematics the question of precisely how much can and cannot be done with- 
out the help of a Higgs field, particularly with respect to dynamical mass generation. 
All elementary fields here have asymptotically detectable particles, so the requirements 
on and interpretation of propagators are clean and simple. 

The same environment may be well suited for studying the question of when-and- 
how-and-how-much the coupling-resurgent amplitudes break symmetries, particularly 
the local gauge symmetries. Here I may remind you of the fact that quantum local sym- 
metries place much more stringent invariance requirements on the functional measure 
of the generating functional ( ^.47 ) than the simple translation invariance we used for 



deriving the DS equation ( |3.52D . The DS equations, consequently, are more fundamen- 
tal than the Slavnov- Taylor identities that express local gauge invariance on the level of 
quantum correlation functions; they reserve the right of producing symmetry-breaking 
solutions. This possibility becomes particularly acute in connection with nonperturba- 
tive solutions containing A terms - since the self-consistency of the latter, as we saw in 
sect. 3.2, is a genuine quantum effect arising from ultraviolet-divergent loops, they have 
in principle the power to break classical symmetries, much like anomalies in that respect. 
On the other hand the perturbative limits of amplitudes, denoted r^^)^''^''* in eq. (|o!^), 
do preserve the quantum-symmetry identities, and that is all that is both necessary and 
sufficient for maintaining the most important practical consequence of gauge invariance 
- the perturbative renormalizability |WEI96 |. I could therefore imagine (why not spec- 



ulate a little on a casual stroll like this) that a picture of " subasymptotic symmetry 
breaking" might emerge in general, where perturbative r('')^"^''*'s as large-momentum 
limits, with their proven Slavnov- Taylor compliance, uphold renormalizability, while the 
A terms coming into view one by one on moving down the momentum scale increasingly 
deviate from Slavnov- Taylor. In QCD, things are probably different - the simple fact 
that gluons obviously do not get saddled with b- or even t-quark masses (through virtual 
quark-antiquark pairs) would seem to require much stricter adherence to Slavnov- Taylor, 



and the fact that the solutions of refs. | DRI99 ] are not doing very well in that respect 



is another indicator that the approximation level there is still too crude. 

A purely methodological but definitely challenging cluster of things-not-done cen- 
ters around the question of whether the DS equations, as employed here, really offer the 
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best framework for the reconstruction of resurgent amplitudes. Clever rearrangements 
of the pure-and-simple DS system, as it emerges from eq. ( 3.52 ) and is sketched in Fig. 
3.1, have long been known. Already in his first paper on the subject | DYS4S| ], Dyson 
used what may be called a Bethe-Salpeter resummation of the QED three-point vertex 
equation, which allows one to replace the bare three-point vertices at the left ends of 
diagrams (3B) and (3 C) in Fig. 3.1 by dressed ones, while replacing the right-hand por- 
tions of both diagrams by a Bethe-Salpeter kernel K4 one- and two-particle irreducible 
in the horizontal channel. Such a rearrangement presumably shifts a substantial amount 
of physics from higher into lower loop orders, and may make a crucial difference when 
studying problems of symmetry breaking. It also provides partial, though not complete, 
relief from a perennial nuisance in most DS equations with > 3 - the lack of manifest 
Bose or Fermi symmetry due to the unsymmetric distinction of the leftmost leg that cor- 
responds to the primary functional differentiation in eq. (|3lg| ). In a non-abelian gauge 
theory, however, one pays for the progress with equally substantial new complication: 
the information on the four-point amplitude, already by far the most complicated of 
the superficially divergent ones, now enters the system in two different ways, as r4 and 
as Ar4, whose relation to each other must be stated as an additional dynamical equa- 
tion. This roughly doubles the effort required for dealing with the superficially divergent 
subsystem. Still more sophisticated rearrangements, originating in work by Symanzik 
|SYM61] and carried out for non-abelian gauge theory by Baker and Lee |BAL77|, allow 
the elimination of bare vertices in favor of dressed ones to be continued down into the 
propagator equations (thus creating a system with the remarkable property that its per- 
turbative iteration does not produce diagrams with overlapping ultraviolet divergences). 
But it once more brings substantial additional complication in that internal lines now 
carry momentum derivatives of propagators in addition to the propagators themselves. 
In the search for improvements along these lines, just about everything seems to come 
at a heavy price. None of these improvements touches the one feature of all DS-type 
dynamical equations that the mathematician studying resurgence would find most disap- 
pointing: that they make no direct reference, either differential or integral or otherwise, 
to coupling dependence, which enters them only in a somewhat indirect, parametric way. 
Equations like the t 'Hooft- Weinberg RG equations, which directly react to coupling de- 
pendence through their trademark P (g) terms, would be much more after the heart 
of the resurgence theorist, but while formulating important restrictions on correlation 
functions, they are not statements of their full dynamics. Equations combining the 
coupling sensitivity of the t 'Hooft- Weinberg with the full dynamical content of the DS 
equations could perhaps be found among various versions of the exact RG differential 
equations of Polchinski [POL84]; one example that I happen to know of is ref. |WIE97|. 
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Equations of this kind could be well suited for studying nonperturbative structure of the 
resurgent type. Once more the subject has not even begun to be explored. 

Still strolling on, we dimly perceive from a distance a set of problems that lie at 
the borderline to mathematics, and may well be the deepest and most long-term ones. 
I emphasized at the outset that I would be appealing only to a few simple notions from 
what may be called the surface of resurgence theory, not to its deeper constructs and 
methods. But to bring these into play in the elucidation of nonperturbative coupling 
structure may be fascinating in the short and unavoidable in the long run. To mention 
one example, the ultraviolet-renormalon poles on the left-hand, real Borel axis that we 
found to persist in our resurgent model amplitude arc presumably a common feature 
of most if not all correlations in asymptotically free theories. In trying to character- 
ize quantitatively the features of coupling dependence controlled by these singularities, 
and in bringing to the fore the common traits of different correlations arising from them, 
Ecalle 's alien-derivative calculus, with its beautiful constructive tool of resurgence mono- 
mials, may have a nontrivial contribution to make. Another conjecture brings us, once 
more, to QCD. There, the quasi-perturbative corrections are not small enough at the 
scale A to be treated by "pragmatic" semiconvergence (as one would not hesitate to 
do in electroweak theory). Accordingly, the ultimate methods of solution may well be 
those that abandon completely the n-versus-p-direction dichotomy still prevailing in our 
quasi-perturbative scheme, and perform partial resummations that cut across the en- 
tire resurgent symbol. The prominent role, aready noted, that resurgent symbols with 
support on the non- negative integers play in Ecalle 's treatment of differential equations 
may have heuristic significance here: perhaps the resummations needed may be differ- 
ential approximants in the coupling, resurgent solutions to a well-conceived sequence of 
differential equations, characterizable by holomorphic invariants of presently unknown 
physical meaning ... 

Our little stroll has led us on to flights of imagination. For here and now, in my 
opinion, the greater value of resurgence theory lies not in specific techniques, intriguing 
as they may be. It lies in providing a framework for thinking, and in connecting dis- 
parate pieces of physical knowledge. The study of QFT, to the extent that it has used 
continuum methods, has been dominated conceptually by the perturbation expansion to 
such an extent that "summing the series" was all but synonymous with establishing the 
full solution for generations of physicists. Resurgent functions provide a much wider yet 
precise framework, wide enough presumably to cover the strongly nonperturbative cou- 
pling dependence of asymptotically free theories. They quite literally open up a second 
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dimension. Prom this, I would expect, will arise their longer-term impact. 
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